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Abstract: Non-human primates (NHPs) are phylogenetically close to humans, with many similarities in terms of physiology, 
anatomy, immunology, as well as neurology, all of which make them excellent experimental models for biomedical research. 
Compared with developed countries in America and Europe, China has relatively rich primate resources and has continually aimed to 
develop NHPs resources. Currently, China is a leading producer and a major supplier of NHPs on the international market. However, 
there are some deficiencies in feeding and management that have hampered China’s growth in NHP research and materials. 
Nonetheless, China has recently established a number of primate animal models for human diseases and achieved marked scientific 
progress on infectious diseases, cardiovascular diseases, endocrine diseases, reproductive diseases, neurological diseases, and 
ophthalmic diseases, etc. Advances in these fields via NHP models will undoubtedly further promote the development of China’s life 
sciences and pharmaceutical industry, and enhance China’s position as a leader in NHP research. This review covers the current 
status of NHPs in China and other areas, highlighting the latest developments in disease models using NHPs, as well as outlining 
basic problems and proposing effective countermeasures to better utilize NHP resources and further foster NHP research in China. 
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Employing animal models for biomedical research 
has been done to speed along research by bypassing 
ethical restrictions on using humans as experimental 
subjects, which averts risks associated with direct 
research on humans, especially in clinical trials. There 
are likewise numerous advantages to using animal 
models such as strong controllability of experimental 
conditions, high repeatability, ease of scale, comparabi- 
lity of results. However, the effectiveness of animal mod- 
els is hampered by intrinsic species differences between 
animal models and humans. For example, experimental 
data obtained from rodents, who are phylogenetically 
distant from humans, does not always appropriately 
predict the efficacy of a drug treatment or its potential 
toxicity in humans (Van Der Worp et al 2010; Xu, 2011). 
To bridge the gaps, non-human primates (NHPs) have 
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been widely used in lab settings since the 1950s. There 
has been a corresponding increase in the prominence and 
importance of using NHP models that have been widely 
used as subjects in infectious diseases (Kaushal et al, 
2012; Liu & Zhang, 2010), mental and neurological 
disorders (Perretta, 2009), cardio-cerebrovascular diseases 
(Sy et al, 2014), and endocrine diseases (Kamath et al, 





Received: 11 June 2014; Accepted: 15 August 2014 

Foundation items: This work was supported by the National Natural 
Science Foundation of China (81172876, 81273251, U1202228, 
81471620); the National Special Science Research Program of China 
(2012CBA01305); the National Science and Technology Major 
Project (2013ZX10001-002, 2012ZX10001-007) and the Knowledge 
Innovation Program of CAS (KSCX2-EW-R-13, KJZD-EW-L10-02) 


“Corresponding author, E-mail: zhengyt@mail.kiz.ac.cn 


Volume 35 Issue 6 


448 ZHANG, et al. 


2011; Pound et al, 2014), to name just a few. 

The key advantage of using NHPs is that NHPs are 
phylogenetically the most proximate models to humans, 
bearing many distinctive similarities not found in other 
animal models. Previous reports found that DNA 
sequence similarities between NHPs and humans can 
reach up to 98.77% (Fujiyama et al, 2002) and the 
concordances rates of gene transcription levels of spleen, 
peripheral blood monocyte cells (PBMCs) and liver 
between NHPs and humans can reach as high 91.41%, 
84.36%, and 74.29%, respectively (Lu et al, 2008). As a 
consequence, NHPs remain the first or the only choice 
for investigating major human diseases including 
HIV/AIDS, measles, malaria, hepatitis, but also in the 
study of human cognition and brain diseases (Xu et al, 
2013). In addition, NHPs are also used for human cancer 
research (Xia & Chen, 2011). The data garnered from 
such NHP experiments is invaluable as a starting point, 
because it ensures greater efficiency and surety of future 
clinical application, making NHPs the “gold standard” 
for preclinical studies (He et al, 2013). 

It was half a century ago that the NHP model was 
gained traction in Europe and America where most of the 
They 
began to develop the NHP resources by importing NHPs 


scientific research were conducted at the time. 


from other countries, developing husbandry programs 
and exploring the nature of NHPs. As more recent 
advances in medical research and biotechnology invol- 
ving gene modification, gene knockout, epigenetics, 
precision gene editing and fine regulation of local gene 
expression in experimental primates become more cost- 
effective and refined, we expect that there will be greater 
need for NHP resources to meet the demands of many 
new applications, which developed countries with 
advanced knowledge of NHPs may struggle to meet. At 
this point, countries with indigenous NHPs that can also 
advance effective protection, sustainable exploitation and 
utilization of their NHP resources may stand to benefit. 
Indeed, as its growing life sciences and pharmaceutical 
industries, China has a great opportunity to become a 
leader in NHP resources and even research (Cyranoski, 
2004; Hao, 2007). 

This review evaluates the ongoing situation of 
NHPs in China and other areas as well as the latest 
development of human disease models using NHPs in 
China. We analyze the existing problems for further 
developing NHP models and propose targeted counte- 
measures that may allow for better exploitation and 
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utilization of NHP resources in China. 


Status of experimental primates abroad 

There is a disagreement on the total number of 
primate species in existence. According to the 2005 
Mammal Species of the World (http: //vertebrates.si.- 
edu/msw/mswCFApp/msw/index.cfm), there are 376 pri- 
mate species. However, the latest statistics from IUCN 
Red List of Threatened Species (http: //www.iucnredlist.org) 
declare 420 primate species. Occasionally, new primate 
species are described, such as Rhinopithecus strykeri that 
was discovered recently (Geissmann et al, 2011; Long et 
al, 2012), make the numbers subject to change. More 
important to the total number of primates are threats 
from humans, whose development often interferes with 
primate habitat. As a result, most primate species are 
considered to be endangered. To overcome the threats to 
NHPs, numerous countries have adopted strict protective 
measures on NHP resources, and some have began to pay 
greater detail to utilizing these resources more rationally. 

While America has no wild NHPs, it is home to 
primate research institutes including 8 national primate 
research centers (NPRCs) funded by the National 
Institutes of Health (NIH) that develop NHP models for 
basic and applied studies of human health: California 
National Primate Research Center, New England Primate 
Research Center, Southwest National Primate Research 
Center, Tulane National Primate Research Center, 
Washington National Primate Research Center, 
Wisconsin National Primate Research Center, Yerkes 
National Primate Research Center, and Oregon National 
Primate Research Center (Kaiser, 2013). These 8 NPRCs 
have a total number of 26, 000 NHPs (Hayden, 2008), 
representing more than 20 species of NHPs, mostly 
macaques (largely for practical purposes and ease of 
research). Annually the US federal government invests 
about $100 million each year into these NPRCs to 
financially support ongoing studies of NHPs (http:// 
www.humanesociety.org/animals in _research/general_in 
formation_on_animal_research/an_introduction_to_ prim 
ate_issues.html). As advances in research techniques and 
methodology using primate models have been created, 
there has been a consummate increase in investment. 
From 1999 to 2006, funding through the US federal 
government for NHPs research increased 7.9-fold, ($15.4 
million to $122 million), while the total number NHPs at 
the eight NPRCs increased from 24 182 to 27 914 (http: 


//www.all-creatures.org/saen/res-nprcs-8yr.html). Other 
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institutions with strong focuses on NHP research also 
receive federal and state funding, including the New 
Iberia Research Center at the University of Louisiana at 
Lafayette, the University of Texas MD Anderson Cancer 
Center, Alamogordo Primate Facility and Primate 
Foundation of Arizona, the Caribbean Primate Research 
Center, a squirrel monkey colony of the University of 
South Alabama and a baboon breeding facility of the 
University of Oklahoma. Even a few private facilities in 
NHP breeding receive government grants to maintain the 
US’s overall capacity for conducting NHP research (http: 
//www.sourcewatch.org/index.php/National_Primate_Re 
search Center System). One further aspect of NHP 
research in the US worth noting is that the US is home 
to the largest colony of Chimpanzees used for research, 
with around 1 000 currently in government sponsored 
labs (Wadman, 2011la). Chimpanzees are the only 
great apes known to be captured for biomedical 
research, even though they are also an endangered 
species, and given that they are the primates most 
closely related to humans, are used in invasive rese- 
arch, and defined as “inoculation with an infectious 
agent and/or drug testing.” 

Compared to the US, Europe is the next largest 
center of NHP research (Carlsson et al, 2004), though 
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with some different trends. While the number of NHPs 
used in experiments in America has grown annually from 
2001 to 2010 (Figure 1, data from http: //www.aphis.- 
usda.gov), the total amount of NHPs used for research in 
the EU remains lower than in America, without stable 
growth. Within EU, total NHPs used for research in 2003, 
2005 and 2008 were 10 360, 10 450 and 9 600, 
respectively (https: //www.gov.uk). The most popular 
NHPs used for research were rhesus macaques (RMs) 
and cynomolgus macaques (CMs), with pig-tailed 
macaques, squirrel monkeys, African green monkeys, 
common marmosets, cotton-topped tamarins, baboons, 
chimpanzee, among others taking a minor share. No apes 
were used in experiments anywhere in the EU from 
2002-2008 according to a survey by the European Com- 
mission (http: //vertebrates.si.edu/msw/mswCFA pp/msw/- 
index.cfm). There is much debate regarding whether 
National Institutes of Health (NIH) research chimp- 
anzees be retired (Wadman, 2011b) and US primate 
centre faces scrutiny with stricter criteria (http: //www.- 
nature.com/news/us-primate-centre-faces-scrutiny- 1.10317). 
In June 2013, NIH announced that it will retire to 
sanctuary nearly all of its research chimpanzees (http: 
//blogs.nature.com/news/2013/06/nih-retires-most-rese- 
arch-chimpanzees.html). 
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Figure 1 Number of NHPs used in experiments in America per year from 2001 to 2010 


While the US and Europe are the key centers of 
NHP research, Southeast Asian countries that possess 
indigenous NHPs have recently begun to pay more 
attention to these resources. Southeast Asia has made 
rapid advances in its NHP capacity: several NHP 
research institutions and breeding facilities have been set 
up in the area, and a few NHP breeding facilities have 
passed the American Association for Accreditation of 
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Laboratory Animal Care (AAALAC) certification. For 
example, Thailand has NHP facilities and primate 
research centers, among which Primate Research Insti- 
tute of Thailand (PRIT; http: //www.primate.or.jp/se- 
mi/Abst_f2nd/suchinda.pdf) possesses 14 species of 
primates, and focuses on primate breeding and producing 
animal models suitable for diseases research, especially 
for tropical diseases research. These advances have made 
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the PRIT a key biomedical research center on primates 
and an information center for NHPs among Southeast 
Asia. and their 
corresponding resources greatly enhances the entire 


The presence of such facilities 


region’s NHP research capacity, and complements the 
growth in NHP resources and research being conducted 
in China. 


Status of experimental primates in China 

Compared to US and Europe, NHP resources are 
rich in China. More than 24 species were found. It 
accounts for about 6% of primate species around the 
world. Yunnan Province in particular possess the most 
NHP resources, and undertakes a bulk of the earliest and 
present work on NHPs (Shen & Yang, 2003; Pan & Ben, 
1984; Li & Lin, 1983). The anchor institution of this 
research, Kunming Institute of Zoology (KIZ), Chinese 
Academy of Sciences (CAS), was established in the late 
1950s. KIZ is among the first institutions in China to 
start a primate breeding center to conduct basic research 
and clinical applications on NHPs in taxonomy, ecology, 
anatomy, genetics, behavior, evolution, reproductive 
biology, neurobiology, immunology, and with a series of 
laboratories dedicated to animal microbial detection, 
animal nutrition, clinic pathology of animal; KIZ has 
also produced a series of monographs on NHPs including 
The Anatomy of the Golden Monkey, The Anatomy of 
Macaca mulatta. The primate centre in KIZ hosts 
currently about 2 800 NHPs, including M. mulatta, M. 
leonine, M. fascicularis, M. assamensis, M. arctoides, 
Rhinopithecus roxellana, and other species. The facility 
is well-equipped to deal with these NHP resources, as it 
has received AAALAC certification, and the centre’s 
NHP resources and research has received a worldwide 
attention (Cyranoski, 2004; Hao, 2007). Alongside KIZ, 
Institute of Medical Biology (IMB), Chinese Academy of 
Medicine Science, which was set up in 1958 and before 
1980s hosted the greatest number of RMs in China (Shen 
& Yang, 2003), focuses on the research and development 
of medical products using experimental primates, 
particularly for vaccines. IMB is currently the largest 
base of the research and production of both attenuated 
poliovirus oral live vaccine and hepatitis A vaccine in 
China. 

Prior to China’s development of domestic NHP 
research into exploitation and utilization of NHP reso- 
urces, China was—and remains—a leading producer and 
a major supplier country for NHPs to the United States, 
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Europe, Japan, and now Korea. In 2013, the number of 
breeding experimental monkeys in China, mainly rhesus 
macaques and cynomolgus macaques, in China there 
were about 250 000 breeding cynomolgus macaques and 
40 000 breeding rhesus macaques (data from The China 
Laboratory Primates Breeding and Development Asso- 
ciation), with approximately 50 000 to 60 000 experi- 
mental monkeys born in breeding facilities each year. 
Though many NHPs are used for domestic research 
projects, a large number are exported abroad. For 
example more than 65% of the experimental monkey 
imports in the US from Asia are from China, and in 
Europe, around 70% of research primate imports are 
from China. Essentially, though NHP research remains 
predominately a Euro-American venture, China is the 
leading supplier of resources necessary to this research. 
To support China’s shift from an NHP resource 
producer to an NHP research driver, the Chinese 
government has increased investment for experimental 
primates and new drug development. Over the last 
decade, there has been marked progress on research and 
innovation into experimental primates. Specifically, the 
Ministry of Science and Technology (MOST) of the 
People's Republic of China has set up an experimental 
primate germplasm resources center in Hainan and 
Suzhou, and in 2011 MOST launched a number of large 
science projects under the platform construction of 
promoting the development of primate animal models of 
human major diseases. The General Administration of 
Quality Supervision, Inspection and Quarantine of China 
has likewise set up a state key laboratory for experi- 
mental primate animal quarantine in Guangxi Province, 
and as we mentioned previously the CAS established the 
Kunming Primate Research Center of Chinese Academy 
of Sciences in Yunnan. In addition, the governments are 
providing more support and financial resources to NHP 
research. In short, these actions will undoubtedly 
accelerate the innovation and capacity for NHP research 
in China. Prior to 2005, most high-quality experimental 
NHPs were exported to European, America and other 
“developed” countries, while the domestic consumption 
was relatively limited, but thanks to the additional 
support and development of the Chinese scientific 
community, outsourcing of services related to 
experimental NHPs have continued to spring up and 
domestic services have exploded in Beijing, Kunming, 
Shanghai, Suzhou, Chengdu, Guangzhou, among others, 
to improve the development of experimental NHPs in 
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China. 


Current status and progress on the NHP models of 
human diseases in China 

Chinese researchers have established numerous 
NHP models to model a variety of human diseases, 
which have made large contributions to the study, preve- 
ntion and treatment of human diseases. In China, most 
NHP models for human diseases focus on exploring 
infectious diseases, cardio-cerebrovascular diseases, 
endocrine diseases, reproductive diseases, mental and 
neurological disorders, and ophthalmic diseases. Simi- 
larly, some basic research into NHPs by Chinese 
researchers have greatly extended the usability of NHPs 
as animal models by exploring transgenic strategies in 
experimental primates and in 2010 researchers of 
Kunming Institute of Zoology, CAS created the first 


transgenic monkey in China (Niu et al, 2010). 


Primate animal models for infectious diseases 

Though advances in medicine and medical 
technology have effectively controlled or eradicated 
infectious diseases, a variety of new infectious diseases 
have continuously emerged, largely as a result of human 
activities’ adverse impact on the surrounding environ- 
ment and bringing individuals into closer contact with 
these situations. As a developing country with the worl- 
d's largest population, China is actually threatened by 
many infectious diseases that are exacerbated by 
population density, intense land-usage, and industrial- 
ization, such as HIV/AIDS tuberculosis (TB), viral 
hepatitis, foot-and-mouth disease, avian influenza. This 
problem is only magnified on the global scale, especially 
among developing countries that experience the bulk of 
population growth and the increase in population 
densities as a consequence of urbanization and 
industrialization. Given the nature of infectious diseases 
and the larger global trends, developing more ideal and 
targeted animal models is a crucial step for opening up 
new strategies against them. 


HIV/AIDS 

To our knowledge, no animal models for HIV/AIDS 
research are more suitable than NHPs. HIV-1, for 
instance, infects only humans and a handful of primate 
animals such as the chimpanzee, gibbon and pigtailed 
macaque (Kuang et al, 2009). While HIV-1 can infect 
these NHPs, for several reasons the disease's progression 
is quite different from that in humans, though AIDS 
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models of RMs infected with SIV closely resemble the 
human disease progression of AIDS. Consequently, 
NHPs are the most widely used animal models for HIV/ 
AIDS research (Lei et al, 2013). Though NHP AIDS 
models were established using simian retrovirus (SRV) 
in the early years, NHP AIDS models created with SIV 
better simulate human AIDS. Currently, most NHP AIDS 
models in China for developing vaccines and 
therapeutics are established via infection of SIV or SHIV, 
a simian/human immunodeficiency chimeric virus. 

An initial Chinese RM/SIV model established by Li 
et al (2007) was accomplished by infecting Chinese RMs 
with SIVmac239. Furthermore, Xia et al (2010, 2011) 
and Li et al (2012) investigated the dynamic and 
functional changes of dendritic cell subsets and Treg 
cells in SIVmac239-infected Chinese RM. The results 
confirmed a previous study and showed that disease 
symptoms of SIV-infected Chinese RMs were more 
similar to human responses to HIV-1. Feng et al (2007) 
also set up RM models of SIVmac251 infection using 
several different routes (intravenous, intravaginal and 
intrarectal routes) and observed differences in disease 
course after SIVmac251 infection resulting from the 
different routes; notably, this was the first case of a 
Chinese NHP AIDS animal model infected with SIV 
through a mucosal route. Later, Cong et al (2011b) 
established a Chinese RM model of SHIV-sf162p3; an 
HIV-1 subtype B, infected through repeated low-dose 
intravaginal inoculation, which more closely resembles 
the natural transmission of HIV infection in humans. 
Though HIV-1 subtype C comprises the predominant 
strains in China (Wu et al, 2006), Cong et al’s model 
provides technical reference for establishing further 
AIDS animal model for evaluating vaccines and 
treatment strategies suitable for China's predominant 
epidemic strains, which may change given growing 
infection rates and variances, especially in Yunnan. 
Additionally, researchers have also successfully 
established different HIV domestic Chinese RM models 
(Wang et al, 2011b; Cong et al, 2011a; Yao et al, 2011). 
The characteristics of these models remains to be 
elucidated in future, however, as HIV/AIDS infections 
continue to increase globally in general and in China 
(and Yunnan) in particular, these models may become 
increasingly important research models. 

While the existing SIV/SHIV-infected macaque 
models are useful for application in AIDS research, 
genetic differences between SIV/SHIV and HIV-1 
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significantly limit the use of SIV/SHIV (Ambrose et al, 
2007; Van Rompay, 2012). Accordingly, NHP AIDS 
models established with HIV/HSIV infection should be a 
better model of human HIV/AIDS infections and 
pathology; fortunately, researchers at the KIZ are 
currently establishing northern pig-tailed macaque AIDS 
model with HIV/HSIV. If successful, such a novel model 
may provide an ideal platform for the study, drug 
screening and the vaccine evaluation of AIDS. 


Tuberculosis (TB) 

Recent studies reveal that mouse models are quite 
questionable, but NHP models avoid many of the pitfalls 
that make rodent models less popular; in particular, NHP 
models of TB for immunological research and vaccine 
test show promising results (Kaushal et al, 2012; Guo & 
Ho, 2014). The study of NHP models of TB was later 
than other countries in China, despite the ongoing 
prevalence of the disease. For TB, there is a certain gap 
in the types and construction methods of NHP models, 
for example, the NHP models of TB latency and 
reactivation have not been done in China. 

Currently, there are several reports on the ongoing 
construction of NHP models for TB study in China. For 
instance, Wuhan University began to establish of NHP 
models of TB using fiberoptic bronchoscopes since 2006, 
and later delivered Mycobacterium tuberculosis (MTB) 
into the lungs of Chinese RMs (Xian et al, 2010; Zhang 
et al, 201la). This laboratory has become the fifth 
standard laboratory for TB vaccine study in the world 
(http: //news.sohu.com/20060715/n244270131.shtml), and 
models have been used in preclinical testing of candidate 
vaccines and drugs. Wang et al (2012) later established 
Chinese RM models of TB via bronchoscopic and 
intratracheal instillation into the lungs and found the 
clinical manifestations, disease progression and 
pathological changes of these models closely resembled 
human primary TB and hematogenous disseminated TB. 

NHP models of TB have marked advantages in 
immunological research and test of drugs as compared 
with non-primate models. Unfortunately, the lack of 
currently available inbred strains of NHPs makes 
pathological changes varied. The obvious solution is to 
strengthen the system research of pathological features, 
such as lobe fibro-cavitary lesions, etc. Likewise, studies 
suggest that there is a high rate of HIV and MTB co- 
infection, which is thought to be one of factors 
contributing to the global resurgence of TB (Harrington, 
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2010; Kaushal et al, 2012). In theory, infection with HIV 
accelerates the activity of latent MTB and then increases 
the incidence of TB and, at the same time, TB seems to 
speed the HIV replication rate (Enserink, 2001), making 
them a potentially lethal complication. Establishing NHP 
models of TB then provides a foundation for the constr- 
uction of NHP models of MTB/HIV co-infection, which 
may prove more ideal than NHP models of MTB/SIV co- 
infection (Diedrich et al, 2010) currently in place. 


Influenza virus 

Of the many infectious disease pandemics, infl- 
uenza has been one of the worst (Kilbourne, 2006). This 
is especially poignant for China, as outbreaks of highly 
pathogenic H5N1 avian influenza viruses occurred in 
Hong Kong in 1997, killing 6 people (Mounts et al, 
1999), which was the first worldwide report of avian 
influenza virus directly infecting humans. Following the 
report, avian influenza A infection in humans was also 
reported in Southeast Asia, Europe and North Africa 
(Obenauer et al, 2006), and outbreaks of H7N9 in 
humans in China last year received local and global 
attention (Gao et al, 2013a, b; Hu et al, 2013; Liu et al, 
2013). Since rodent animals are not the natural hosts of 
influenza viruses and show no symptoms of upper- 
respiratory-tract infections (Zak & Sande, 1999), rodent 
models are not good systems for studying human 
influenza infections. NHP models are, however, excellent 
at mimicking human infections of influenza viruses. 
Foreign researchers have established NHP models 
infected with H1N1, H3N2, H5N1, etc (Baas et al, 2006; 
Carroll et al, 2008; Cilloniz et al, 2009; Fan et al, 2009; 
Itoh et al, 2009; Kobasa et al, 2007; Laddy et al, 2009; 
Rimmelzwaan et al, 2001). Recently, Chen et al (2009) 
infected Chinese RMs with the highly pathogenic avian 
influenza virus A/Tiger/Harbin/01/2002 (H5N1) intrana- 
sally and noted that RMs exhibited symptoms similar to 
humans, such as fever, loss of appetite, interstitial 
pneumonia, etc. These results indicated that Chinese 
RMs could be potentially successful infection models of 
the H5N1 viruses as they exhibit similar symptoms. To 
date, this model has some shortcomings, since H5N1 
only reproduced in the respiratory system and that RMs 
showed low sensitivity to H5N1 infection. In view of 
these limitations, Lü et al (2011) infected Chinese RM 
with A/SZ/406H/2006(H5N1) by endotracheal incuba- 
tion, and found that this method yielded a more sensitive 
NHP model of H5N1 virus infection, and that the virus 
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could replicate in many other organs outside the lung and 
better resembling severe cases seen in humans. 

Since the influenza virus is primarily transmitted via 
air, it is necessary to assess the contributions of small 
droplet aerosols and larger respiratory droplets to better 
understand epidemiological vectors for the disease and 
public responses to influenza virus outbreaks. Technol- 
ogical hurdles currently make it difficult, such as lack of 
inoculations for influenza viruses in NHPs and the diffi- 
culty in isolating and quantifying viable influenza viruses 
from dilute aerosols make this currently impractical. 


Hepatitis virus 

Compared to other infectious diseases, China has a 
long history of research into hepatitis. Over 30 decades 
ago, Ge et al (1989) performed long-term liver biopsies 
of M. assamensis, M. Speciosa, Nycticebus Coucang 
infected with HBV, successfully setting up NHP models 
of HBV infection. Zhang et al (1998) infected RMs with 
HGV and established the first RM model of HGV 
infection, while more recently Yang et al (2006) infected 
Chinese RMs with HEV gene type-IV via intravenous 
injection. Following Yang’s injection, Chinese RMs 
exhibited typical acute hepatitis, indicating Chinese RMs 
made viable animal models of HEV gene type-IV. 
Chinese researchers also successfully established 
Chinese RM models of HEV1 infection, and through 
these models found E2 IgG can serve as a marker of 
HEV infection (Zhang et al, 2003; 2004). 

Chinese researchers have used HAV, HBV, HEV 
and HGV to infect NHPs (Ge et al, 1989; Zhang et al, 
1998; Zhang et al, 2003; 2004). But there have been no 
reports of the establishment of domestic NHP models for 
HCV infection. Chimpanzee models of HCV infection 
have already been established by foreign researchers, 
though this is expected since chimpanzees are the only 
known animal to be susceptible to HCV infection (Bukh 
et al, 2001). However, ethical and financial restrictions 
largely limit the use of chimpanzees for such research, 
especially in the US and Europe (Akari et al, 2009). 
GBV-B, a close relative to HCV, provides an effective 
model and has been used to infect New World monkeys 
and the GBV-B infection models which closely resemble 
HCV infection in humans (Akari et al, 2009; Bukh et al, 
1999). In spite of the value of the GBV-B model, the 
establishment and development of an HCV/GBV-B 
chimeric virus model has higher application value for a 
preclinical study of antiviral vaccines and drugs, and 
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accordingly foreign researchers have created monkey 
models infected with an HCV/GBV-B chimera (Li et al, 
2014; Rijnbrand et al, 2005). In future, we should focus 
on developing NHP surrogate models for HCV infection. 


Enteroviruses 

Enteroviruses (EV) affect millions worldwide each 
year. The enterovirus is a genus of RNA viruses assoc- 
iated with several human diseases, and among them 
Enterovirus 71 (EV71) is highly contagious and causes 
foot and mouth disease (HFMD) , which has increasingly 
been a serious threat for Asian children (Lee & Chang, 
2010). CMs infected with EV71 via lumbar injection or 
intravenous route show some symptoms similar to those 
found in humans, but these CM models do not exhibit 
certain important infection symptoms, such as 
myocarditis and neurogenic pulmonary edema. Likewise, 
most of the EV71 applied to these CM models are host- 
adapted EV71 strains (Hashimoto & Hagiwara, 1982; 
1983; Hashimoto et al, 1978; Nagata et al, 2004; Nagata 
et al, 2002), and there can be marked differences 
between host-adapted EV71 strains and the epidemic 
strains. Unfortunately, NHPs are not susceptible to 
infection with epidemic strains of EV71 ( Liu & Zhang, 
2010). 

The study of EV71 infection NHP models began 
relatively late in China, with the first models appearing 
in 2010. Zhang et al (2011b) infected adult Chinese RM 
with EV71 FY-23 strain via intracerebral, intravenous, 
respiratory and digestive routes and showed similar 
symptoms to human cases but without herpetic lesions, 
the reason for which may be old-age of these 
experimental Chinese RMs. Perhaps neonatal monkeys 
would prove better animal models for the study of EV71 
pathogenesis, and accordingly Liu et al (2011) infected 
neonatal RMs with EV71 FY-23 strain via the respiratory 
tract by tracheoscopy and the results indicated these 
neonatal RM models more closely resembled infection of 
human neonates during lactation. Wang et al (2011a) 
likewise infected neonatal RMs with EV71 via nasal 
spray, which further confirmed that a respiratory 
infection route makes models more susceptible to EV71 
infection than the digestive tract infection route, 
providing some valuable information for future studies of 
childhood EV71 transmission. 


NHP models of cardio-cerebrovascular disease 
Myocardial infarction, a type of cardiovascular 
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disease, has high morbidity and mortality, and is com- 
mon in Europe and America (Yeh et al, 2010). Previous 
studies indicated that the treatment strategy—proved to 
be effective in a murine model—applied to NHPs failed 


to show any efficacy (Norol et al, 2003; Orlic et al, 2001). 


Yang et al (2011) recently established RM models of 
myocardial infarction by ligating the left anterior 
artery (LAD) and the 
histopathology and electrocardiogram demonstrated that 


descending changes of 
this model should be excellent for studying human 
myocardial infarction, but the main cause of myocardial 
infarction in human is cardiovascular blockage or spasm 
caused by the long-term effects of several factors 
(hypertension, hyperlipidemia, etc.). The pathogenesis of 
these RM models created by ligation of the coronary 
artery existed many differences compared with human, 
which currently limits drug discovery and development. 
It is necessary that new NHP models should be 
developed for studing human myocardial infarction in 
future. 

Cerebrovascular disease is the second leading cause 
of mortality in China. Ischemic cerebrovascular disease 
accounts for 80% of cerebrovascular diseases (Feigin et 
al, 2003). Pan et al (2006) injected polyvinyl alcohol into 
a branch of middle cerebral artery to make thrombed, 
and CT scanning and related behavioral changes showed 
significant damage to brain tissue in the area of vascular 
occlusion. These results indicate that NHP models of 
ischemic cerebral infarction were successfully establi- 
shed and they owned good stability and were suitable in 
the study of treating acute and chronic diseases. Zhu et al 
(2008) also established CM models of photochemical 
thrombosis by focusing the cold lamp-house on the 
precentral gyrus, and the construction of these models is 
simple and highly reproducible, making them excellent 
for the basic study and investigation of treatments for 
cerebral infarction, especially regarding the study of 
transplantation of stem cells. Xu et al (2012) recently 
generated Chinese RM models of middle cerebral artery 
occlusion by injecting auto-blood clots into middle 
cerebral artery through the femoral artery with a micro- 
catheter, a method that is easily repeatable, with high 
stability and low mortality. Likewise, the extent of 
embolization of the cerebral infarction could be timely 
detected. The aforementioned advantages have led to this 
technique becoming more common and ideal for the 
study of cerebral ischemia. Zhu et al (2009) previously 
established CM models of intracerebral hemorrhage by 
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injecting antiautologous blood slowly and gently into the 
caudate nucleus. This study indicated that the volume, 
site and speed of injection are key factors for establishing 
a successful model, which is instructive to development 
of the NHP models of intracerebral hemorrhage. 

While there has been significant progress in 
establishing NHP models of cardio-cerebrovascular 
disease, some shortcomings still exist. The NHPs chosen 
to be animal models for cardio-cerebrovascular disease 
are usually young and healthy, despite the fact among 
humans the elderly are at the highest risk for cardio- 
cerebrovascular disease (Baccini et al, 2008). Similarly, 
the patients with cardio-cerebrovascular disease usually 
suffer from other health conditions such as hypertension, 
hyperlipidemia or diabetes which may significantly 
affect cardio-cerebrovascular symptoms. Currently avail- 
able studies of NHP models of cardio-cerebrovascular 
disease that accompany these additional diseases is quite 
limited, to create viable NHP models of cardio-cerebro- 
vascular disease which is close to clinical nature will be 
our future task. 


NHP models of diabetes 

Diabetes is already quite prevalent in developed 
countries, mainly type 1 and type 2, and it is quickly 
rising in China now, with an estimated 11% of the 
population currently afflicted to some extent. Although 
rodents have been used as the primary model in studying 
diabetes, rodent models do not adequately mimic human 
diabetes and make them quite limited for clinical 
application (Anderson & Bluestone, 2005; Delovitch & 
Singh, 1997). NHPs, due to their similarities with human, 
are believed to be the most effective models for human 
diabetes. 

NHP diabetic models are relatively new in China. 
Liu et al (2009) divided the Chinese RM into a 125 
mg/kg STZ group, a 75 mg/kg STZ group and a 50 
mg/kg STZ group and then different doses of STZ were 
injected into corresponding subjects intravenously. The 
models induced by the dose of 125 mg/kg and 75 mg/kg 
STZ were not stable, and mortality was high due to the 
toxic effects of STZ. Most subjects in the 50 mg/kg STZ 
group retained a high level of blood glucose and their 
pancreas B cells were severely damaged, but their 
mortality was quite low, suggesting that 50 mg/kg of 
STZ may be optimal for inducing diabetes in Chinese 
RM models. Xu et al (2009) likewise injected STZ at the 
dose of 45mg/kg into the CMs intravenously and found 
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that drinking volume, food intake, urine output increased 
and body weight was reduced, while fasting glucose was 
increased and the pancreas showed atrophy and fibrosis, 
indicating that the CM diabetic models created with this 
method were viable. Given the serious toxic effects of 
STZ on experimental NHPs, Qiao et al (2009) and Jin et 
al (2010) established RM models of type 1 diabetes by 
the partial pancreatectomy combined with low-dose STZ. 
This strategy is generally safer and more reproducible as 
compared with single high-dose STZ or total 
pancreatectomy. 

Currently, the most commonly used method to make 
type 2 diabetic animal models is high fat feeding 
combined with drugs induction. Typical clinical features 
of type 2 diabetes was found by injection of STZ into 
CM at a dose of 35 mg/kg after the appearance of 
hyperlipidemia and obesity based on high-glucose and 
fat diet feeding, (Lu et al, 2013) and the models showed. 
Since pathogenesis of type 2 diabetes induced by STZ is 
unnatural as it significantly differs from the pathogenesis 
of human type 2 diabetes, theoretically, screening of 
monkeys for spontaneous diabetes would lead to a more 
natural pathogenesis in NHP models. It will make this 
model more valuable for translational medicine. Wang et 
al (2004) previously screened 3 RMs of spontaneous 
diabetes from 100 RMs by both the glucose tolerance 
and glycosuria tests. These 3 RMs were in middle and 
elderly age, indicating that onset age and the morbidity 
rate of spontaneous diabetes in RMs is consistent with 
that found in humans. After one year of observation, the 
similarity in clinical symptoms (urinating more and 
consuming more aggravated gradually and fasting blood 
glucose was significantly higher than normal) to humans 
suggested this method was feasible for screening RM 
models of spontaneous diabetes. Unfortunately, the 
current number of spontaneous diabetic NHP models is 
limited and cannot adequately meet the growing demand 
for effective scientific research. 

Given that the morbidity rate of human type 2 
diabetes increases annually as long as increases in the 
number of obese people (mainly caused by factors such 
as eating patterns, among others), monkey models of 
type 2 diabetes induced by diet can be accomplished 
more quickly, and compared with models induced by 
drugs, more closely resemble the pathogenesis of human 
type 2 diabetes. Zhang et al (2012) recently divided 36 
experimental CMs into 4 groups and fed them standard 
basic diet, high-sugar diet, high-fat diet and high-sugar 
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combined with high-fat diet respectively for 24 months. 
The progression of type 2 diabetes among CMs with 
high-energy diets showed that similar methods were 
successful in establishing NHP models of type 2 diabetes, 
and that it effectively modeled different stages of type 2 
diabetes. These findings also provided data on further 
optimizing dietary strategies that could establish even 
more effective animal models for human type 2 diabetes 
research. While type 2 diabetic mouse models have been 
made (Kadowaki, 2000), the successful establishment of 
transgenic NHP models of type 2 diabetes should vastly 
improve NHP study of type 2 diabete. 


NHP models of mental and neurological disorders 

To our knowledge, none of animal models can fully 
recapture the neurological symptoms for neurodegene- 
rative and psychological disorders due to their multiple 
origins. While NHPs are close to human in terms of 
genetics, there are numerous differences in the evolution 
and composition in brain. Currently, NHP models of 
mental and neurological disorders such as heroin 
addiction, morphine addiction, Parkinson’s disease (PD), 
depression syndrome, Alzheimer’s disease (AD), spinoc- 
erebellar ataxia, spinal cord injury, motor cortex lesion, 
epilepsy, etc (Chen et al, 2013; Chu et al, 2014; Ding et 
al, 2008; Liu et al, 1992; Lu et al, 1999; Qiao et al, 2007; 
Tian & Ma , 2014; Zha et al, 2006; Zhang et al, 1999 ) 
have been established in China. However, there is a gap 
in the types of NHP models for mental and neurological 
disorders as compared with similar models established in 
other countries. 

An increased understanding regarding natural and 
artificial repair of central nervous system injuries as 
resulted as more researchers have begun focusing on the 
study of nerve regeneration following spinal cord injuries. 
To model these injuries in NHPs, the main strategy is to 
inflict spinal cord injury via selective resection of the 
spinal cord with a sharp knife (Tian et al, 2010). Ni et al 
(2005) made NHP models of spinal cord hemisection by 
T-11 laminectomy and then resection of a 1-mm long 
hemispinal cord, while Zha et al (2006) made the C3-5 
exposed via a dorsal midline incision in the neck of RMs 
and transected the left hemispinal cord with a knife under 
the microscope after C4 laminectomy. The results 
showed this method to establish RM models of spinal 
cord injury were simple, effective, stable, accurate and 
repetitive, even if the method was quite different from 
naturally occurring human spinal cord injuries. For 
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epilepsy, models are usually created by the intentional 
application of penicillin, aluminum hydroxide, metrazol 
and coriaria lactone. Chen et al (2013) injected kainic 
acid into the right intra-hippocampus of RMs by the 
stereotacic technique, and the resulting seizure pattern, 
imaging observation and histopathologic tests showed 
this model closely resembled human temporal lobe 
epilepsy. Though promising, the long-term stability of 
this RM models induced by kainic acid requires further 
study. 

Modeling PD has likewise been met with several 
successes in China. Liu et al (2006) injected 1-methyl-4- 
phenyl-1, 2, 3, 6-tetrahydropyridine (MPTP) into right 
side internal carotid artery and then performed 
behavioral testing and imaging detection, which showed 
the method to be both safe and simple for establishing of 
NHP models of PD. However, this method made it 
difficult to control the extent of nigrostriatal damage 
(Herrero et al, 1991). In order to minimize the 
nigrostriatal damage, Yan et al (2014) established a 
chronic NHP model of PD by intramuscularly injecting 
MPTP. The time of appearance and the severity of each 
symptom of this model could be controlled by the daily 
dosage and frequency. This model gradually displayed 
PD symptoms and mimicked the development of human 
PD well. However, the PD symptoms were reversible in 
some MPTP-induced model, being quite different from 
the degenerative nature of PD in humans, and limiting 
the study of stem cell therapy for nervous system 
disorders. Interestingly though, substantia nigra neurons 
damage caused by 6-hydroxydopamine (6-OHDA) is 
irreversible, and the extent of nigrostriatal damage can be 
controlled by the injection amount and the injection sites. 
Accordingly Zhang et al (2006) generated RM models of 
PD by injecting 6-OHDA stereotaxically into the right 
substantia nigra. In testing potential treatments of PD, 
optimal treatment period is a key issue, and one that can 
be addressed using NHP models of early stage PD. 
Accordingly, Ding et al (2008) established NHP early 
stage models of PD by surgery to make the common 
carotid artery (CCA) of RMs exposed and then injecting 
low-dose MPTP into the exposed CCA. Compared with 
late-stage PD NHPs, NHP early stage models may better 
mimic the gradual development of human PD. One of the 


purposes of model development is to find new treatments. 


Yan et al (2014) reported that morphine could alleviate 
tremors in a MPTP-induced NHP model, suggesting the 
therapeutic effects of morphine complement those of L- 
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Dopa, a classic therapy, on PD symptoms. But due to 
MPTP toxicity to the digestive system, digestive system 
disorders of the NHP models induced by MPTP 
gradually appeared and would affect immune function 
and nutritional metabolism. Such models are difficult to 
maintain over long durations, making the future 
development of a long-term stable NHP PD models an 
urgent issue. 

Recently, studies have established a link between 
increased endogenous formaldehyde (FA), a methanol 
metabolite, and AD pathology (He et al, 2010; Li et al, 
2008; Tong et al, 2011, 2012). Studies conducted by Hu 
and He’s Labs in China on RM showed chronic feeding 
3% methanol ad libitum led to specific damages to the 
brain similar to those of Alzheimer’s disease. These were 
the first studies to show that methanol toxicity inflicted 
all the major hallmarks of AD on rhesus monkey: 
cognitive decline, tau phosphorylation, and amyloid 
plaques formation (Yang et al, 2014). 

Recent advances in optogenetic technology of have 
been proved as a revolutionary tool for studying neural 
systems. Optogenetics can be used for precisely targeting 
a specific area of neural systems as well as even in freely 
moving mammals. Researchers can then turn different 
neurons on or off, rapidly and safely, using this 
technology (Deisseroth, 2010; Liu & Tonegawa, 2010). 
Currently, research of optogenetic approaches has been 
done on NHPs by foreign researchers and some progress 
of this research has been made (Dai et al, 2013), but 
there is lack of such research being done domestically in 
China. 


NHP models of ophthalmic diseases 

Ophthalmic diseases present interesting complic- 
ations for animal models. Compared with other animals, 
the NHP visual system is closer to that of humans than 
any other animal, making experimental data from NHPs 
safer and more effective when translated into clinic 
settings. However, both cost restrictions and availability 
of experimental NHPs hamper the development of NHP 
models of ophthalmic diseases. Chinese researchers have 
established NHP models of 
ophthalmic diseases, such as corneal endothelial injury, 


successfully several 
amblyopia, glaucoma and so on (Dai et al, 2005; Zhu et 
al, 2013), but there are some gaps in the types and 
construction methods of NHP models of ophthalmic 
diseases between China and other developed countries. 
Zhu et al (2013) established NHP models of corneal 
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endothelial injury by phacoemulsification damaging 
corneal endothelial cells. This method proved easy to 
manipulate and is easily replicable, and it causes negli- 
gible damage to other tissues, making it a useful 
technique for constructing models that can accurately 
evaluate the effect of drugs and other factors on the 
(2005) 
established NHP models of chronic hypertensive glauco- 


endothelial healing. Similarly, Dai et al 


ma by semiconductor frequency-doubled 532 laser and 
argon laser, respectively, aiming at 360° functional trab- 
ecular meshwork. The light coagulation times of the 
successful induction of high intraocular pressure in RMs 
fluctuate greatly, so refining the existing methods of 
photocoagulation should improve modeling efficiency, 
though this requires greater study to verify. 


NHP models of reproductive diseases 

Reproductive diseases seriously affect human health 
and even psychological well-being. Infertility caused by 
reproductive diseases affect an average 8-12% of couples 
worldwide (Sellandi et al, 2012), and is quite pronounced 
in developed countries. Reproductive physiology of 
NHPs is quite similar to human, and the data gathered 
from studies on reproductive diseases in human using 
NHPs cannot be replaced using any other experimental 
animals. The progress in the establishment of NHP 
models of reproductive diseases, such as polycystic 
ovary syndrome (PCOS), oligozoospermia, metrorrhagia 
and postpartum hemorrhage (PPH), has been made by 
Chinese researchers. 

Tang et al (2012) gave RMs two cycles of subcut- 
aneous injections of propionic acid testosterone (PAT), 
and then muscle injections of human chorionic gon- 
adotropin (HCG). The results of this study indicated PAT 
combined with HCG could successfully induce RM 
PCOS, which was similar to human PCOS. Sun et al 
(2006) earlier set up an RM model of oligozoospermia 
induced by orally administered gossypol acetic acid and 
provided a suitable NHP model for the study of infertility 
in men. Until recently, most studies on PPH still were 
built on the retrospective analysis, with few reports 
evaluating or using animal models of PPH. Notably, 
Huang (2010) established CM models of PPH with 
uterine atony induced by medication such as oxytocin, 
which provided a useable platform for further research. 
Likewise, You et al (2003) had earlier set up RM models 
of metrorrhagia induced by endometritis and the 
symptoms of this model could closely resemble those of 
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human metrorrhagia caused by endometritis. 


Problems, Challenges and prospects 

While China is a key producer of NHP resources, it 
still lags behind other countries in its capacity for 
innovative NHP model research. As medical research and 
biotechnology advances—genetic modification, gene 
knockout and cloning, etc —1t is likely that a remarkable 
increase of NHP resources will be developed to meet the 
needs of more intensive and diverse application. It makes 
for some interesting possibilities in NHP research. Due 
to the economic crises and the subsequent “austerity” in 
budgets across developed countries in America and 
Europe, funding support for primate research is reduced 
or kept stagnant. Concurrently, growing criticism from 
animal-rights activists in these countries has brought a 
great deal of pressure and research institutions have to 
limit or even abolish primate research (Cyranoski, 2003; 
Wadman, 2011b). Given the advances and advantages 
provided by using NHP primates, researchers on NHP 
models from developed countries are increasingly 
seeking collaborations with primate research institutes in 
China and other developing countries in Southeastern 
2004). It is possible to bring 
unprecedented opportunities and challenges to life 
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sciences and biomedical research in China. The key issue 
is how to cope with these opportunities and challenges to 
further develop NHP resources effectively. 

Development of domestic Chinese NHP resources 
has made significant progress since its inception, 
however, in contrast to countries in America, Europe and 
East Asia (notably Japan and Korea), NHP resources in 
China are not widely or effectiely being used. In other 
words, while China has a rich diversity of species, many 
factors such as subspecies and geographical distribution 
of NHP resources, differences among subspecies in 
genetics and physiology and the occurrence of 
subspecies hybrid populations in breeding and over the 
course of experiments make the reproducability of 
experimental results poor in working towards etablishing 
of animal models, evaluating drugs, etc. Nontheless, 
some common experimental primate resources (espec- 
ially exotic species such as Cercopithecus aethiops) are 
rare or even non-existent in China, which restricting 
human diseases research and the development of 
corresponding drugs, but other more general factors such 
as technical skill, experience, and existing resources are 
also critical. We must realize that Chinese NHP research 
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can only be accomplished by bringing important 
experimental primates lacking in China—and researchers 
familiar with them—while simultanesouly strengthening 
the application existing experimental primates resources 
in China. 

Another factor making Chinese NHPs ineffective is 
a lack of the principle knowledge about NHPs. For 
example, compared with some experimental primate 
species such as Indian RMs, common experimental 
primate species such as Chinese RMs lack sufficient 
information regarding the genetic background of these 
primates, and it is also little known about the normal 
value of their important biological parameters (Xia et al, 
2009). Thankfully, recent genome sequencing and 
analysis of CMs and Chinese RMs offers a great deal 
more basic data and useful reference points for relevant 
NHP model research (Yan et al, 2011). It is necessary to 
complement basic biological data of experimental 
primates in China and establish subspecies breeding 
populations whose basic biological data (such as genetic 
backgrounds) are definite. Recently, the Kunming 
Primate Research Center is making efforts to raise a 
northern pig-tailed macaques (NPMs, M. leonina) colony 
and Prof. Zheng’s lab has collected many basic 
biological data of this colony of NPMs (Pang et al, 2013; 
Zhang et al, 2014), in order to promote the standar- 
dization and reliability of experimental primate data. 

Another key issue for NHP research is the demand 
for specific pathogen-free (SPF) experimental animals. 
SPF animals are in great demand right now. While it is 
available for many traditional models, experimental 
primates often carry zoonotic pathogens which not only 
pose health risks to experimenters but also partly 
interfere with experimental results. In the US, all 8 
NPRCs funded by grants through NIH can offer SPF 
experimental animals (Zheng, 2012), but to date China 
has not established SPF core groups of experimental 
primates. Presuming that research into drug development 
and microbe research on diseases such as HIV continues 
to grow, if China wants to be a leading player in NHP 
models, it must establish SPF core groups of experi- 
mental primates and all breeding efforts should be done 
only with careful caution and planning. 

In recent years, experimental primate husbandry in 
China has experienced rapid development and the 
product has been increased, but not always in tandem 
with the demand. As we mentioned earlier, the influence 
of the global economic crisis in developed countries in 
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2008 as well as interference from animal-rights activists 
have led to a plateau in primate imports from China. On 
the Chinese side, this unexpected lack of growth has led 
to rampant overproduction. Likely one reason for 
Chinese overproduction is that many NHP breeding 
facilities in China lack effective communication with 
their international counterparts, making it difficult to 
project and then match demand. To remedy the situation 
and ensure adequate supply of relevant and in-demand 
NHP resources, it is vital to establish an information and 
communication platform that can effectively offer NHP 
breeding information services and simultaneously carry 
out strategic plans for cooperation between NHP 
breeding facilities and relevant governing bodies and 
institutions. 

Two other factors could affect China’s position on 
NHP research and supply. First, animal welfare activists 
have had some effects on NHP research. While it is true 
that advocacy of animal welfare helps promote harmony 
between experimental animals and humans (and conse- 
quently decreases the danger of injuries to experimenters 
and often offers higher-quality experimental results), 
animal welfare has begun to attract extensive attention 
which will result in increased legislation on animal 
welfare. China has just recently begun to legislate animal 
welfare. On a cultural level, Chinese citizens’ awareness 
of animal welfare is weak. Awareness of animal welfare 
has placed economic pressure on NHP breeding facilities 
and research institutions. The dual-faced nature of this 
dilemma means that effective implementation of 
experimental primate welfare in China is difficult. To 
bring China more in line with international standards 
while continuing to develop China’s domestic research 
capacity, relevant policies should be made for 
strengthening propaganda, supervision, regulation and 
support of experimental primate welfare, in such a way 
to support the healthy development of the life sciences 
and medical industry in China. Second, it is a very 
crucial step to develop transgenic NHP resources. Since 
the birth of the first transgenic mouse, researchers have 
established a series of transgenic mice for modeling 
many types of human diseases. Transgenic mice have 
made greater contributions to the study and the treatment 
of human diseases, but their phylogenetic distance to 
humans hampers the usefulness of their results. 
Establishing primate transgenic models should be able to 
overcome these limitations and better predict the efficacy 
and toxicity of novel treatments in humans. The 
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development of these transgenic primates is still in its 
infancy; the first transgenic primate was born in 
2001(Chan et al, 2001), and over the past decade, few 
transgenic NHPs have been born (Chen et al, 2012). The 
major obstacles in producing transgenic NHPs are 
technical in nature. However, thanks to more efficient 
gene-editing techniques involving ZFN, TALEN and 
CRISPR/Cas9-based methods, researchers can avoid 
some shortcomings—notably low efficiency and high 
toxicity—of previous transgenic methods using retroviral 
or lentiviral vectors (Shen, 2013). Chinese researchers 
first applied TALENs and CRISPR/Cas9 to establish 
site-specific gene mutant RMs and CMs (Liu et al, 2014b; 
Niu et al, 2014). The myriad studies have demonstrated 
involving TALENs and 
CRISPR/Cas9, multiple and single genetic mutations can 


that using technologies 
be effectively performed in monkeys (Liu et al, 2014a; 
Niu et al, 2014). 

In addition to the above-mentioned problems, there 
are still many other key problems with NHP models to 
resolve, for example, high cost, difficulties regarding 
how to have a well-defined phenotype characterization, 
and a genotype-phenotype correlation for NHPs and how 
to use cutting-edge molecular imaging system to 
characterize the NHP brain structure and activity for 
better understanding mental disease. 
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Flow cytometric characterizations of leukocyte subpopulations in the 
peripheral blood of northern pig-tailed macaques (Macaca leonina) 
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Abstract: Pig-tailed macaques (Macaca nemistrina group) have been extensively used as non-human primate animal models for 
various human diseases in recent years, notably for AIDS research due to their sensitivity to HIV-1. Northern pig-tailed macaques (M. 
leonina) are distributed in China and other surrounding Southeast Asia countries. Although northern pig-tailed macaques have been 
bred on a large scale as experimental animals since 2012, the reference value of normal levels of leukocytes is not available. To 
obtain such information, 62 blood samples from male and female healthy northern pig-tailed macaques at different ages were 
collected. The normal range of major leukocyte subpopulations, such as T lymphocytes, B lymphocytes, natural killer (NK) cells, 
monocytes, and the expression levels of activation or differentiation related molecules (CD38, HLA-DR, CCR5, CD21, IgD, CD80 
and CD86) on lymphocytes were analyzed by flow cytometry. The counts of B cells decreased with age, but those of CD8" T cells 
and NK cells and the frequency of CD38°HLA-DR‘CD4" T cells were positively correlated with age. The counts of leukocyte 
subpopulations were higher in males than those in females except for CD4* T cells. Males also showed higher expression levels of 
IgD and CD21 within B cells. This study provides basic data about the leukocyte subpopulations of northern pig-tailed macaques and 
compares this species with commonly used Chinese rhesus macaques (M. mulatta), which is meaningful for the biomedical 
application of northern pig-tailed macaques. 


Keywords: Northern pig-tailed macaque; Flow cytometry; Leukocyte subpopulation; Age; Sex 


Pig-tailed macaques (Macaca nemistrina group, 
PTM), rhesus macaques (M. mulatta, RM) and cynom- 
olgus macaques (M. fascicularis, CM), all belonging to 
the Cercopithecidae family of Old World monkeys, have 
become more widely used models for AIDS pathogenesis 
in recent years (Lei et al, 2013; Zhu et al, 2010). As their 
name suggests, PTM are characterized by the short, 
semi-erect ‘pig-like’ tail and the flat vertex. They 
diverged from RM and CM approximately 5 million 
years ago, while the divergence between RM and CM 
was generated at one another 2.4 million years later 
(Baroncelli et al, 2008). 

PTM have been widely employed as animal models 
of sexually transmitted diseases due to the similarities 
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between the anatomic structure of the genital tract and 
the menstrual cycle with humans, the feature that they 
breed in all seasons and have a relatively large body size. 
By contrast, it is difficult to perform colposcopic 
examinations and multiple biopsies on CM as they have 
much smaller vaginal canals and cervix diameters than 
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PTM. The reproductive physiology of RM is also quite 
different from humans as their constant breeding season 
(Patton et al, 2009). More importantly, PTM have the 
overriding advantages on their abilities in acute human 
(HIV-1) 
exceptional susceptibility to simian immunodeficiency 


immunodeficiency virus-1 infection and 
virus (SIV) and simian-human immunodeficiency virus 
(SHIV). In fact, PTM are the only reported species of 
Old World monkeys that could be infected by HIV-1 
(Agy et al, 1992; Batten et al, 2006; Bosch et al, 2000). 

Although they have been extensively applied in 
biomedical areas, such as cognitive neuroscience, pharm- 
acology and infectious etiology, PTM did not obtain a 
clear taxonomic status until Gippoliti et al (2001) sepa- 
rated the initial taxonomic three subspecies (nemestrina, 
leonina, pagensis) of the M. nemistrina group at the 
species level. Sunda pig-tailed macaques (M. nemestrina) 
are found in the southern half of the Malay Peninsula 
(only just extending into southernmost Thailand), Borneo, 
Sumatra and Bangka Island. Mentawai macaques (M. 
pagensis) are located in the Mentawai islands. Northern 
pig-tailed macaques (M. leonina) range from about N8° 
in Peninsular Thailand through Burma and Indochina 
into Bangladesh, India extending north as far as to the 
Brahmaputra, and the southernmost of Yunnan, China 
(Malaivijitnond et al, 2012). 

Due to the lack of primate models that can be 
productively infected by HIV-1, especially given that 
HIV-1 infection is more efficiently restricted by tripartite 
motif protein 5a (TRIM5a), apolipoprotein B mRNA- 
editing enzyme and catalytic polypeptide-like 3 
(APOBEC3) proteins of macaques than those of humans, 
the pathogenic mechanisms of HIV-1 have not been fully 
elucidated. The early block to HIV-1 infection in monkey 
interference with TRIMSa 
expression (Stremlau et al, 2004). In our previous study, 
we for the first time reported that northern PTM do not 
express TRIMSa but rather TRIM5-CypA fusion gene, in 
which the B30.2/ SPRY domain of TRIM5a is substituted 


by retrotransposed CypA in the 3'-UTR, contributing to 


cells was relieved by 


the invalidation of restriction to HIV-1 replication in 
these animals (Kuang et al, 2009). Recent studies 
showed that PTM rather than RM or CM are susceptive 
to siman-tropic HIV-1 (stHIV-1), in which the HIV-1y74.3 
vif gene was replaced by the vif genes from SIV mac Or 
HIV-2zop to suppress the antiretroviral activity of several 
APOBEC3 proteins (Hatziioannou et al, 2009; 
Thippeshappa et al, 2011). Thus northern PTM may is a 
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desirable animal model for AIDS studies. 

Wild northern PTM mainly distribute in southwe- 
stern Yunnan, China and southeastern Tibet. Because of 
the scarcity of this primate species, they are under First 
Grade State Protection of China. In 2012, Kunming 
Primate Research Center, China, first introduced nort- 
hern PTM from Vietnam for large-scaled breeding for 
laboratory use. However, the limited knowledge of their 
immune systems has restricted their applications. In this 
study, we determined the phenotypic characteristics of 
leukocytes in northern PTM ranging from 2 to 11 years 
of age to provide detailed immunological data. 


MATERIALS AND METHODS 


Animals and sample collection 

Cage-bred northern PTM were maintained according 
to the regulations approved by the American Association 
for Assessment and Accreditation of Laboratory Animal 
Care (AAALAC) at the Kunming Primate Research 
Center, Kunming Institute of Zoology, CAS. Ethyle- 
nediamine tetraacetic acid (EDTA) stabilized blood sam- 
ples were collected from male (n=15) and female (n=14) 
juveniles (mean age=2.4 years (2-3 years)), as well as 
male (n=17) and female (n=16) adults (mean age=6.1 
years (5—11 years)) without apparent external symptoms 
(acute infections, trauma and severe diarrhea) for flow 
cytometric analysis. For comparison purposes, data from 
age-appropriate male (n=9) and female (n=8) juvenile 
(mean age=2.4 years (2—3 years)), as well as male (n=14) 
and female (n=12) adult (mean age=7.5 years (5—10 
years)) Chinese rhesus macaques (ChRM) were also 
collected. Animals’ age and sex information is shown in 
Table 1. All procedures were performed under the 
guidance of the Ethics Committee of Kunming Institute 
of Zoology. 


Antibodies 

Mouse anti-human CD molecule monoclonal and 
polyclonal antibodies that cross-reacted with PTM were 
used according to standard procedures and appropriate 
concentrations in this study. Anti-CD159a PE (clone 
Z199) mAb was obtained from Beckman Coulter. Anti- 
CD38 FITC (clone AT-1) mAb was obtained from STE- 
MCELL. Streptomycin-PE-Cy7 was obtained from Biol- 
egend. Anti-IgD biotin pAb used by conjunction with 
streptomycin-PE-Cy7 was obtained from Southern Biot- 
ech. Anti-CD14 APC (clone MSE2), anti-CD20 PerCP-CyS. 
5(clone 2H7), anti-CD21 APC (clone B-ly4), anti-CD3 
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Table 1 Study group of northern pig-tailed macaques and Chinese rhesus macaques 


Northern pig-tailed macaques 


Chinese rhesus macaques 











Age groups (years) Males Females Males Females 
n Age (years) n Age (years) n Age (years) n Age (years) 
Juvenile (2-3) 15 2.3+0.5 14 2.5+0.5 9 2.3+0.5 8 2.5+0.5 
Adult (4-11) 17 6.6+1.6 16 5.4£1.7 14 7.01.6 12 8.11.5 
Total 32 4.6+2.5 30 4.0+2.0 23 5.2+2.6 20 5.943.0 





APC-Cy7 (clone SP34-2), anti-CD4-PerCP-Cy5.5/ 
FITC (clone L200), anti-CD8a PE-Cy7 (clone RPA-T8), 
anti-CD80 PE (clone HB15e), anti-CD86 PE (clone 
FUN-1), anti-CCRS PE (clone 3A9) and anti- HLA-DR 
APC (clone L243) mAbs were all obtained from BD 
Pharmigen. 


Absolute quantification of major leukocyte subpopu- 
lations 

The methods of direct cell surface staining for whole 
blood and absolute number analysis were as described 
previously (Xia et al, 2009). Briefly, 50 uL of whole 
blood samples were added in TruCount tube (BD Bio- 
sciences) and incubated with fluorochrome-conjugated 
antibodies of CD3, CD4, CD8a, CD20, CD14 and 
CD159a mAbs for 30 min on ice. Erythrocytes were 
lysed with FACS lysing solution (BD Biosciences), and 
the samples were analyzed with BD FACSVerse 
cytometer directly. The absolute numbers of T cells 
(CD3*CD20), CD4* T cells (CD4°CD8a°CD3"), CD8* T 
cells (CD4°‘CD8a°CD3’), B cells (CD3°CD20’), NK cells 
(CD3°CD8a*CD159a’) (CD3 CD20 
CD14*) were calculated using the following formula: 


and monocytes 


Cell concentration (cell numbers/uL)=(events in gated 
region (n)xtotal number of TruCount beads (7))/(number 
of acquired beads (n)xsample volume (uL)) (Figure 1). 


Phenotype analysis by flow cytometry 

As shown in Figure 1, to determine CCR5* T cells, 
activated (CD38°HLA-DR’) T cells, activated (CD80/ 
CD86") B cells, resting (CD21°) B cells and IgD/IgM 
secreting (IgD") B cells, 100 uL of fresh whole blood 
used for each analysis were lysed with FACS lysing 
solution for 10 min at room temperature, followed by 
washing and resuspension with Dulbecco’s phosphate- 
buffered saline (DPBS) with 2% of new-born calf serum 
and 0.09% sodium azide (staining buffer). The suspe- 
nding leukocytes were then stained with the relevant 
directly conjugated mAbs for 30 min on ice and fixed 
using PBS containing 4% paraformaldehyde. 
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Flow cytometric acquisition was performed on the 
BD FACSVerse cytometer driven by the FACSuite 
software (version 1.0.3; BD). The acquired cell number 
was at least 100 000 CD3* T cells or 50 000 CD20° B 
cells. Analysis of the acquired data was performed using 
FlowJo software (version 7.6.1; TreeStar). 


Statistical analysis 

Statistical differences between two groups were 
determined by the nonparametric Mann-Whitney U rank 
sum test. The Spearman’s test was performed for corre- 
lation analysis. Two-tailed P<0.05 was considered statis- 





tically significant. Data were presented as mean+SD. All 
statistical analyses were performed via GraphPad Prism 
software (version 5.01; GraphPad Software). 


RESULTS 


The absolute number of leukocyte subpopulations in 
peripheral blood of northern PTM 

The absolute number of T cells, CD4* T cells, CD8* 
T cells, B cells, monocytes and NK cells in different 
groups are listed in Table 2. Numbers of B cells were 
negatively correlated with age (7=—0.426, P<0.001). 
CD8" T cell counts (7=0.336, P=0.008) and NK cell cou- 
nts (r=0.329, P=0.01) were slightly increased with the 
increasing of age. Adult group (1 117+548/uL) had 
dramatically less B cells than juvenile group (1 782+ 








741/uL, P<0.001), however, no significant cell count 
differences of CD8* T cells and NK cells were found 
between adult and juvenile groups. 

In adult groups, compared with females, males 
showed higher counts of T cells (male 3544+1341/uL, 
female 2 875+1192/uL; P=0.032), CD8* T cells (male 
1310+456/uL, female 1 0464533/uL; P=0.026), B cells 
(male 1594+621/uL, female 1251+789/uL; P=0.016;) 
and NK cells (male 701+617/uL, female 378+303/uL; 
P=0.004). However, no significant differences of age or 


























sex were found in their counts of CD4* T cells and 
monocytes. 
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Figure 1 Gating strategies for flow cytometric analysis 
A: Gating strategies for T cells, B cells, NK cells and monocytes. CD14/side-scatter scattergram used for gating monocytes, TruCount beads and remaining 
lymphocytes. CD3/CD20 dot polt gated on lymphocytes used to count T and B cells. T cells were further divided into CD4* and CD8" T cells. NK cells were 
defined as CD3°‘CD8a'CD159a" phenotype gating on lymphocytes. B: The expression levels of CCR5 were analyzed in each T cell subset. C: CD38'HLA-DR* 
cells in each T cell subset were defined as activated cells. D: The expression levels of CD21, CD80/CD86 and IgD were analyzed in B cells. 


Table 2 Cell counts (cell numbers/pL) of leukocyte subpopulations in different age and sex groups 




















Age (years) Sex T cells CD4" T cells CD8" T cells B cells NK cells Monocytes 
Juvenile M 3 43341 469 1 886874 2094533 1 7514698 5484411 4424319 
(2-3) F 3 524+] 322 1 933+701 2954633 1 815+809 4214368 5774330 
Adult M 3 642+1 254 19414662 3994597 1 4554526 8364741 6364412 
(4-11) F 2306+701 1 303+491 828+304 7574281 3404240 4144235 
Total M 3 54441 341 1 915+756 3104567 15944621 701+617 5454378 
F 2 87541 192 1 5974668 0464533 1251+789 3784303 4904290 
All 3 220+1 305 1 7614727 1824562 1 4284722 5634510 5184337 
Range 1 141-6 635 564-3 794 321-3 001 378-3 282 66-3178 154-1 844 
Expressions of activation or differentiation related increased progressively with the increaseing of age in CD4* 
molecules on lymphocytes of northern PTM T cells (=0.256, P=0.045), and more CD38°HLA-DR’ sub- 
The expression values of activation or differentiat- set within CD4* T cells (juvenile 1.3+0.4%, adult 2.0+0.9%; 











ion related markers on each lymphocyte subset in diffe- P=0.006) as well as CD8" T cells (juvenile 4.1+1.8%, adult 
rent groups were shown in Table 3. CD38 'HLA-DR* subset 5.0+1.9%; P=0.02) were observed in adult groups. 
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Figure 2 Comparisons between northern pig-tailed macaques and Chinese rhesus macaques 
Table 3 Expression levels (%) of activation or differentiation associated markers in different groups 
Age gee CD38"DR* CD38"DR* CCR5* CCR5* CD21* IgD* CD80/CD86* 
(years) CD4* T cells CD8* T cells CD4* T cells CD8* T cells B cells B cells B cells 
Juvenile M .250.31 3.85£1.46 3:53 £215 11.14+6.36 67.99 + 12.76 72.75 £12.93 44.794 12.51 
(2-3) F 3240.53 4.3242.21 3.48 + 1.98 11.8544.77 64.67 + 11.53 65.71 £13.97 48.19 +9.48 
Adult M 2.27 +0.92 5.49+2.18 4.07+ 1.45 12.02 +£3.73 66.81 + 12.05 69.97 + 11.94 4L.78411.11 
(4-11) F 63 £0.81 4.584 1.50 4.01 £2.20 14.3547.19 60.03 + 14.13 56.124 13.89 60.144 17.42 
Total M 79 £0.86 4.72 £2.03 3.824 1.81 11.61 +5.07 67.374 12.20 71.27 £12.29 43.19+ 11.69 
F .49 +0.70 4.46 1.84 3.764 2.08 13.19+6.21 62.194 12.98 60.60 + 14.53 54.57+ 15.30 
All .64 +0.79 4.59 + 1.93 3.79£1.93 12.37 £5.66 64.86 + 12.75 66.11414.35 48.70+ 14.61 
Range 0.49-3.84 1.48—10.70 1.00-8.79 3.42—28.80 33.80-90.50 35.80-93.10 17.00-80.10 
Males showed higher frequency of IgD* B cells em PTM 3 220+1 305/uL, ChRM 3 998+1 561/uL; 
(male 71.3+12.3%, female 60.6+14.5%; P=0.006) but P=0.006), less CD8* T cells (northern PTM 1 182+562/uL, 





lower frequency of CD80/86* B cells (male 43.2+11.7%, 
female 54.6415.3%; P=0.004) than females. However, 
no statistically significant age and gender differences in 
the frequencies of CCR5'CD4* T cells, CCR5°CD8* T 
cells and CD21'B cells were found. 





Comparisons of lymphocyte subpopulations between 
northern PTM and ChRM 


As shown in Figure 2, a flow cytometric analysis 
was performed to study the differences in lymphocyte 
subpopulations between PTM and ChRM. Compared 
with ChRM, northern PTM exhibited less T cells (north- 
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ChRM 1 634+921/uL; P=0.009), whereas higher freque- 
ncies of both CD4* T cells (northern PTM 55.1+7.3%, 
ChRM 47.3+10.1%; P=0.0001) and CD21" B cells 
(northern PTM 69.4+12.8% compared to 54.2+18.8%; 
P=0.001) as well as increased CD4/CD8 ratio (northern 
PTM 1.62+40.61, ChRM 1.38+0.72; P=0.016). Further 
analysis demonstrated that male northern PTM had less T 
cells (northern PTM 3 544+1 341/uL, ChRM 
4498+1458/1L; P=0.03), less CD8* T cells (northern 
PTM 1 310+567/uL, ChRM 1 9854991/uL; P=0.01; 
northern PTM 36.846.9%, ChRM 44.0413.8%; P=0.013), 
but higher frequency of CD4° T cells (northern PTM 
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54.347.2%, ChRM 45.9+11.0%; P=0.003) and CD4/CD8 
ratio (northern PTM 1.58+0.65, ChRM 1.25+0.75; 
P=0.01) than male ChRM. Whereas, females showed 
lower frequency of CD21 B cells (northern PTM 
37.8+13.0%, ChRM 54.6+15.4%; P<0.0001) compared 
to female ChRM. 




















DISCUSSION 


PTM group (M. nemestrina, M. leonine and M. 
pagensis) have served an important role in the field of 
biomedicine, notably in HIV/AIDS research, in recent 
years (Lei et al, 2013). However, no detailed taxonomic 
description of these animals was found in most studies 
that use the collectively called “pig-tailed macaques” 
(Agy et al, 1992; Batten et al, 2006; Bosch et al, 2000). 
We previously demonstrated that northern PTM express 
TRIM5-CypA fusion gene, resulting in its sensitivity to 
HIV-1 and improving the application value of these new 
laboratory animals in HIV/AIDS research (Kuang et al, 
2009). We also provided the reference values of immu- 
noglobulins, complements, C reactive protein (CRP), 
hematological and biochemical indexes of north PTM, 
and elucidated that gender, age and weight can influence 
these indexes (Pang et al, 2013; Zhang et al, 2014). 
However, the lack of the immune system characteristics 
of northern PTM may have restricted its application as 
AIDS animal models. Here we evaluated the absolute 
number and the expression of activation or differen- 
tiation related markers of major lymphocyte subpop- 
ulations in 62 male and female northern PTMs ranging 
from 2 to 11 years of age, and compared the lymphocyte 
subpopulations of PTM with those of the most widely 
used ChRM. 

Serving as a major risk factor of morbidity and 
mortality for many infections caused by various patho- 
gens, age influences the immune system in both human 
and non-human primates (NHP) (High, 2004). Wiener et 
al (1990) found that the absolute number of both CD19” 
B cells and CD3* T cells decreased while the absolute 
number of CD4* and CD8" T cells and the CD4/CD8 
ratio were well maintained with increasing age. The 
following Swedish OCTO immune longitudinal study 
showed decreased CD4 subset, increased CD8 subset, 
and a lower CD19" B cell percentages in aged population 
(Ferguson et al, 1995). Fagnoni et al (1996) reported a 
decrease in the absolute number of both CD4* and CD8* 
T cells in people older than 60 years. However, no age 
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correlations of the absolute or relative cell numbers of 
lymphocyte subpopulations were found in recent study 
(Klose et al, 2007). Although the results were contr- 
adictory in these studies, a widely accepted theory 
suggests that aging is characterized by the decline in B 
cell numbers and an accumulation of CD8" T cells rather 
than a loss of CD4 T cells in peripheral blood (Frasca et 
al, 2008; Koch et al, 2007). On the other hand, Studies 
have suggested that aging is associated with numerous 
alterations in innate immunity which is the first line of 
defense against pathogens and plays a key role in 
regulating the responses of adaptive immunity (Shaw et 
al, 2010). Positive correlations were observed between 
age and CD16°CD56* NK cell numbers, and it is 
generally accepted that the CD56dim peripheral NK cell 
population expands with age (Jiao et al, 2009; Mahbub et 
al, 2011). Tollerud et al (1989) evaluated the influence of 
age, race, and gender on the immune system, and did not 
observe age related effects for CD14" cells. However, 
other recent research reported that the number of 
monocytes was significantly higher in the elderly than in 
the young group (Della Bella et al, 2007). 

Because they share greater genetic and physiological 
similarities with humans than rodent models, NHP has 
been used for biomedical research for several decades. 
Most of our understanding of the immune system in Old 
World monkeys comes from studies utilizing rhesus 
macaques (Messaoudi et al, 2011). However, only few 
reports focused on the age-related changes in the 
immune system of PTM, especially northern PTM. 
Asquith et al (2012) measured age-related changes in T 
cell homeostasis in Indian rhesus macaques (InRM) 
ranging from 1 to 30 years old, and indicated that the 
frequency of CD8* T cells expands with age while the 
frequency of CD4” T cells remains stable, whereas age- 
related changes of the absolute number of both CD4* and 
CD8* T cells were not found in InRM (Didier et al, 
2012). The CD4/CD8 ratio in juvenile InRM is 
approximately 2, and starts to decline after the first 2 
months of life, whereas no further decline is detected 
between 5 and 7 years of age (Dykhuizen et al, 2000). 
Similarly, B cells in RM are identified based on the 
expression of CD20 and HLA-DR, and show a sharp 
decrease during aging (Didier et al, 2012). Early studies 
indicated that the number of circulating CD56°CD16° 
NK cells in RM reduces with age (Coe & Ershler, 2001). 
However, highly different from humans, NK cells in 
peripheral blood of macaques express high levels of 
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CD8a and almost all of them express NKG2A, but do 
not express CD56, which make them a CD3°CD8a° 
NKG2A\ phenotype (Hong et al, 2013). A recent study of 
rhesus macaques ranging from 2 to 24 years old did not 
find any correlation between the circulating NK cells and 
age (Didier et al, 2012). Though the frequency of 
circulating monocytes increases in InRM older than 5 
years, the numbers of these cells reduces with aging 
(Asquith et al, 2012). 

It has been reported that CD38°HLA-DR™ subset 
within CD8* T cells has a strong positive correlation with 
HIV-1 disease progression and CD38°HLA-DR* CD4* T 
cells are highly susceptible to infection with R5-tropic 
virus because of elevated CCRS5 expression (Kestens et 
al, 1992; Meditz et al, 2011). CCRS5 plays an essential 
role in HIV pathogenesis as a main co-receptor 
responsible for viral transmission. CCR5*CD4* T cells 
are the target of R5-tropic virus and peripheral blood 
CCR5*CD8* T cells accumulate during HIV infection 
(Brenchley et al, 2004). In addition, a loss in the 
expression of CD21 and IgD, and increased levels of 
activation markers, including CD80 and CD86, on circul- 
ating B cells are suggested be linked to HIV viremia 
(Moir & Fauci, 2008). Klatt et al (2012) demonstrated 
that higher pre-infection levels of immune activation 
than RM may contribute to rapid disease progression in 
PTM. Thus it is also necessary to examine these in both 
healthy humans and laboratory animals. The frequency 
of CCR5* cells within CD4° T cells in northern PTM 
(3.79+1.93%) was the lowest compared to Sunda PTM 
(14.5446.3%) and RM (7.76+4.8%) (Klatt et al, 2012). 
The frequency of CD38°HLA-DR’ cells was 4.4+2.5% 
for CD4* T cells and 14.5+7.0% for CD8* T cells in 
healthy human as compared with 1.64+0.79% and 
4.59+1.93% for CD4* and CD8* T cells, respectively, 
from northern PTM and 3.03+2.84% for CD8" T cells 
from rhesus macaques (Brignolo et al, 2004; Onlamoon 
et al, 2005). Elevated levels of CD21 B cells are 
believed to account for the poor proliferative responses 


























of B cells. This population is significantly expanded in 
most patients with systemic lupus erythematosus (SLE) 
and HIV-1 and its normal level is approximately 15% in 
healthy people, 45% in CM, 46% in RM and 36% in 
northern PTM (Das et al, 2011; Kling et al, 2011; 
Zandieh et al, 2013). The B cells in the peripheral blood 
of healthy people have an average 75% IgD* subset, 
while northern PTM have lower frequency (66%) (Klein 
et al, 1998). 
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ChRM infected with SIV commonly have lower 
plasma viral loads and develop AIDS more slowly than 
InRM, whereas PTM are susceptible to many lentiviral 
infections, including infection with SIV, HIV-1, HIV-2, 
SHIV and stHIV-1, and their progress to AIDS after SIV 
infection are also more rapid than RM (Ling et al, 2002). 
A comparison of lymphocytes between ChRM and 
northern PTM may provide clues to explaining this 
phenomenon. Previous studies demonstrated that ChRM 
experience significant age-related changes in leukocyte 
subpopulations, including decreased lymphocyte 
numbers, CD4* T cell frequency and CD4/CD8 ratio, 
increased CD8* T cell frequency, but stable monocytes 
(Xia et al, 2009; Zheng et al, 2014). Similarly, northern 
PTM showed increased numbers of NK cells and CD8* T 
cells, reduced B cells, but stable CD4* T cells and 
monocytes in our study. In addition, northern PTM also 
showed higher CD4* T cell frequency but fewer CD8* T 
cells and CD21 B cells than ChRM. 

The sex effects on circulating lymphocytes have 
been reported in both humans and NHP. Uppal et al 
(2003) demonstrated that women have a higher number 
and frequency of CD4° T cells, whereas men have a 
higher frequency of CD8° T cells. However, no 
differences in circulating B cells, NK cells and CD14° 
cells were found between men and women (Choong et al, 
1995; Tollerud et al, 1989). However, the influence of 
sex on the leukocytes in macaques is less clear. 
Tryphonas et al (1996) revealed significant gender 
differences in the CD4* T cell percentage (females> 
males) and CD4/CD8 ratio (females>males) of infant 
CM. Recent studies on circulating lymphocytes in InRM 
did not find differences between males and females 
(Asquith et al, 2012; Didier et al, 2012). However, our 
previous study revealed that female ChRM have higher 
counts of CD4* T cells, CD8* T cells, B cells and NK 
cells than males (Xia et al, 2009). In this study, northern 
PTM also showed numerous gender differences in 
leukocytes and activation markers, and males may have 
stronger immune response due to their higher levels of 
CD8* T cells, B cells and NK cells than females. Howe- 
ver, further work is required to determine the character- 
istics of immune function between males and females. 

In summary, the present study not only provided 
data on the immunological characteristics of the 
leukocyte subpopulation, but also demonstrated age and 
sex effects on those in northern PTM. Moreover, the 
comparison of northern PTM with ChRM provided clues 
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to understanding the big differences between these two 
species when infected with HIV-1 or SIV. 
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Abstract: Seasonal variation in environmental factors is vital to the regulation of seasonal reproduction in primates. Consequently, 
long-term systematic data is necessary to clarify the birth seasonality and pattern of primates in highly seasonal environments. This 
study indicated that black-and-white snub-nosed monkeys (Rhinopithecus bieti) at Mt. Lasha exhibited strict birth seasonality with a 
pulse model. Infants were born with a certain degree of synchronization. Birth distribution showed three birth peaks, and the birth 
pattern showed a “V” style in even-numbered years and a gradual increase in odd-numbered years. The beginning date, end date and 
median birth date were earlier in even-numbered years than those in odd-numbered years. The higher latitude of their habitats, earlier 
birth date, shorter birth period, fewer birth peaks and stronger birth synchrony might be adaptations for strongly seasonal variation in 
climate and food resources. After the summer solstice when daylight length began to gradually shorten, R. bieti at Mt. Lasha started 
to breed during the period with the highest environmental temperature and food availability, which implied that photoperiod may be 
the proximate factor triggering the onset of estrus and mating. It appears that R. bieti coincided conception and mid-lactation with the 
peak in staple foods, and weaning with the peak in high quality of foods. Thus, food availability was the ultimate factor regulating 
reproductive seasonality, and photoperiod was the proximate factor fine-turning the coordination between seasonal breeding and food 
availability. 


Keywords: Rhinopithecus bieti, Birth seasonality; Birth pattern; Birth synchrony 


In wild environments, birth seasonality exists in 
many animals. Primates show reproduction continuum 
from seasonality to un-seasonality like other mammals 
(Brockman & van Schaik, 2005). The typical reproduc- 
tion patterns of primates include seasonality (births 
within a specific period), peaks (birth peaks in certain 
months) and irregularity (births within the entire bree- 
ding cycle without regularity) (Struhsaker & Leland, 
1987; Brockman & van Schaik, 2005). In high and 
moderate birth seasonality, >67% and 33%—67% of infa- 
nts were found to be born within three months, respecti- 
vely (van Schaik et al, 1999). External environmental 
factors (e.g. temperature, precipitation, food resources, 
photoperiod, seasonal climate, and food availability) and 
biological rhythms affect reproductive seasonality (Bron- 
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son & Heideman, 1994; Negus & Berger, 1972; Andelm- 
an, 1986; Bronson & Heideman, 1994; Brockman & van 
Schaik, 2005; Tecot, 2010; van Schaik & Pfannes, 2005). 
In natural environments, seasonal food resources are the 
most important factor that determine birth pattern, while 
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precipitation, temperature and daylight cycle function as 
proximate factors (Wikleski et al, 2000). Animals 
synchronize their critical reproductive period with the 
optimal phase of environments to maximize survival 
(Negus & Berger, 1972). Photoperiod and temperature in 
temperate regions, as well as precipitation in tropical 
areas, affect seasonal food resources, and eventually 
induce reproductive seasonality (Vivien-Roels & Pevet, 
1983). Under some circumstances, however, photoperiod 
has no influence on conception, and food availability is 
the only determining factor of reproductive seasonality 
(Huang et al, 2012; Carnegie et al, 2011; Kowalewski & 
Zunino, 2004). With increases in both latitude and 
altitude, climate and food seasonality change signific- 
antly. In certain primate species (e.g. Macaca thibetana), 
the median birth date occurs earlier with the increase in 
altitude (M. fuscata: Fooden & Aimi, 2003; P. entellus: 
Newton, 1987). 
individual’s weight gain indicates the accumulation of 


During the breeding season, an 
nutrients, which not only guarantees the nutritional req- 
uirements of mating individuals, infants and nursing 
activities, but also ensures that the critical period of 
infant survival falls into the next autumn with enriched 
food resources (Zhao & Deng, 1988). Reproductive 
seasonality is more obvious among species living in 
areas with high altitude and latitude (Janson & Verdolin, 
2005; R. brelichi: Yang et al, 2009; R. roxellana: Qi et al, 
2008; Ren et al, 2003; Zhang et al, 2000; P. entellus: 
Newton, 1987; C. polykomos: Dasilva, 1989; P. pileatus: 
Stanford, 1990). In addition, reproductive peaks can be 
significantly correlated with climate and food seasonality 
(e.g. P. entellus: Newton, 1987), while relatively stable 
food acquisition (e.g. provision or crop stealing) can 
induce reproductive irregularity (Bishop, 1979). 
Black-and-white snub-nosed monkeys (Rhinopi- 
thecus bieti) are endemic to the Trans-Himalayas, which 
are bound by the upper rivers of the Yangtze to the east 
and the Mekong to the west (N26°14'—N29°20') (Long et 
al, 1994). From north to south, the altitudes of groups of 
this species decrease, habitat transits from dark coni- 
ferous into dark coniferous-broadleaf mixed forest, and 
gradient trends manifest in photoperiod, temperature, 
precipitation and food availability (Xiang & Sayers, 
2009; Ding & Zhao, 2004; Huang et al, 2012). The birth 
period of the Xiaochangdu group (N29°15') was found to 
occur from February 4 to March 14, with the median 
birthdate on February 24 and standard deviation of 6.5 
days (Xiang & Sayers, 2009). The birth period of the Mt. 
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Lasha group (N26°20') was found to occur from 
February 15 to April 7 in 2009-2010, with the median 
birth date on March 27 (Huang et al, 2012), but was from 
February 19 to April 12 in 2011, with the median birth 
date on March 17 (Wang et al, 2012). It appears, 
therefore, that reproductive patterns are influenced by 
differences in habitat altitude and latitude, with monkeys 
in high altitude and latitude areas showing earlier birth 
dates but shorter breeding periods, as well as by 
differences year to year. Mating behaviors are triggered 
by various factors, such as photoperiod (Sadlier, 1972), 
so for some species living in highly seasonal 
environments, maintaining high energy levels may not 
always be necessary for conception (Brockman & van 
Schaik, 2005; Drent & Daan, 1980). 

Reproductive seasonality in black-and-white snub- 
nosed monkey groups (especially the Xiaochangdu group) 
is not induced by photoperiod, but is determined by the 
enriched food resources during the conception period 
(Huang et al, 2012; Xiang & Sayers, 2009). Reproduction 
of R. beiti is also characterized by birth synchrony, with 
spatial-temporal differences (Huang et al, 2012; Wang et al, 
2012). Reproductive synchrony refers to the way species 
maximize food quality and availability with high energy 
consumption periods (e.g. nursing and weaning) (van 
Schaik & van Noordwijk, 1985). Black-and-white snub- 
nosed monkeys, also called snow monkeys (Elliot, 1912), 
inhabit areas with high altitude (Long et al, 1994), as well 
as significant seasonal climate changes and food resources 
(Kirkpatrick et al, 1998; Xiang & Sayer, 2009; Huang et al, 
2012). However, limited short-term data have been unable 
to provide a strong interpretation of the reproductive 
seasonality and corresponding environmental factors in the 
breeding of black-and-white snub-nosed monkeys. In this 
study, reproductive data of the Mt. Lasha black-and-white 
snub-nosed monkey group were collected to determine the: 
(1) reproductive seasonality; (2) birth pattern (peaks and 
variations); and (3) mechanisms of regulating reproductive 
seasonality. 


MATERIALS AND METHODS 


Study areas and animal subjects 

Our study site is located at Mt. Lasha (N26°20', 
E99°15’) in the Yunling National Reserve, Lanping 
County, Nujiang Prefecture, Yunnan Province, China. 
The highest peak is 3 854 m in elevation. Only one 
black-and-white snub-nosed monkey population with 
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130 individuals has been observed at Mt. Lasha, and is 
comprised of eleven one-male multi-female units (OMUs) 
and two all-male-units (AMUs) (Huang et al, 2012; 
Wang et al, 2012). Home range of this population covers 
an area of 14 km? and is surrounded by villages and 
farmland at low altitude (<2800 m), as well as alpine 
pastures at high altitude (>3 600 m). Its habitat is patchy 
with sporadic pastures, fire grasslands, and charcoal 
burning lands (from oak). The vegetation from low to 
high altitude forms zonal distributions of coniferous 
forest, mixed broadleaf-conifer forest, broad-leaved 
deciduous forest and dark coniferous forest (Huang et al, 
2012; Wang et al, 2012). 


Annual average temperature at Mt. Lasha is 





12.42+6.81 °C, with lowest and highest monthly average 
temperatures of —1.88+7.15 °C (—4.15—-1.00 °C, January) 
and 25.84+5.97 °C (15.86-34.06 °C, October), respect- 
ively. An extreme temperature fluctuation occurred from 
August 2012 to July 2013 (4.15-34.16 °C). Annual 
precipitation is 767.60 mm, with the highest proportion 








falling from March to September (695.60 mm), and can 
be subdivided into precipitation peak (July to September) 
and sub-peak (March to May). The snow season occurs 
from November to the following March. 


Data collection and statistical analyses 

Data collection was scheduled ahead of the monkey 
breeding season (from February to April) (Huang et al, 
2012; Wang et al, 2012). At panoramic view points (e.g. 
ridges and rocks) opposite monkey activation areas, we 
recorded numbers of infants and female adults when they 
are crossing an open area and foraging in the wide-field 
forests using a monocular telescope (Leica, T77). 
Newborns were characterized by overall white hair 
(expect a hint of dark hair on the head), black tip of the 
tails (one third), as well as red and hairless auricles and 
extremities. Newborns less than 1.5 months old are 
carried by adult females. To ensure the accuracy of the 
data, only observations with more than half of adult 
females of the group were used for further analysis since 
its interbirth interval is two years (Huang et al, 2012; 
Wang et al, 2012). When multiple observations were 
carried out in one day, only the most reliable (e.g. best 
view, most individuals included) was used in statistical 
analysis. Reproductive seasonality was analyzed by 
considering the numbers of infants and the ratio of 
infants to female adults (I/F). 

We used each 7-day period to calculate all births, to 
infer the median birth date (Md), mean birth date (Mn) 
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and standard deviation (SD) (Caughley, 1977, but see 
Eisenberg et al, 1981). We then transformed each day 
into degrees of a circle (360°) for testing birth seasonality 
by circular statistics, the mean vector length r was 
calculated to measure birth season dispersion, ranging from 
0 (uniform) to 1 (clustered) (Batschelet, 1981). The 
Rayleigh test (Z=nr*) was used to determine whether the 
birth data were unevenly distributed throughout the year 
(Batschelet, 1981; Zar, 1999; Huang et al, 2012). 

The I/F data of the Xiaochangdu group obtained 
from either close tracking recording or field scope 
scanning were used to determine birth pattern. When an 
I/F was obtained by both methods at the same time, the 
former had priority (Xiang & Sayers, 2009). 


RESULTS 
Birth timing 

The first R. bieti infant at Mt. Lasha was born on 
February 20, 2013. The I/F value indicated that the 
breeding season ended on April 14. The numbers of new- 
born infants peaked (13) on April 8, and thereafter, 
stabilized on 13. Because the I/F value was rather unstable 
than the J-value, and the F-value also included adult 
females that had given birth in the last year, the I-value 
was the more reliable index reflecting the birth timing. 
The end of the breeding season inferred from the I-value 
was April 8 (Figure 1A). Therefore, in 2013, there were 13 
new-bom infants at Mt. Lasha and the breeding season 
was from February 20 to April 8. By referring to the I- 
values, the birth data of R. bieti at Mt. Lasha from 2009 to 
2012 was analyzed (Figure 1). Although the beginning of 
the breeding season was consistent in the two estimations 
based on the I-value and I/F, the latter was overall 
postponed with deviations (—1 to 8 days) (Table 1). 


Table 1 Breeding seasons of Mt. Lasha monkey group 
inferred from different indexes 








Jei Beginning of End of breeding season 
breeding season I/F I Deviations (day) 

2009 2-21 4.02 4-03 zi 

2010 2-15 4.04 4-04 0 

2011 2-19 4.20 4-12 8 

2012 2-17 4.07 4-06 1 

2013 2-20 4.14 4-08 6 


I: Infants; F: Female adults; I/F: Ratio of infants to female adults. 


Birth seasonality 
Strict birth seasonality with a pulse type (SD<30 
days) was observed in R. bieti at Mt. Lasha. The 
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breeding season lasted 57 days. Reproductive variables 
had certain similarities every two years (e.g. beginning 
phase, ending phase, and Md). The breeding season 
began and ended earlier in even-numbered years than 
those in odd-numbered years, and Md was postponed by 
ten days (Table 2). 
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exhibited certain synchrony, and showed three peaks on 
February 15-26, March 11—19 and March 20—April 12, 
respectively (Figure 2). The data on infant births using 3- 
day intervals indicated that the third peak occupied 50% 
of the entire breeding season, whereas the first and the 
second peaks accounted for 25% and 23.5%, respectively. 
The highest birth peak in both odd- and even-numbered 















Birth pattern years fell in the third peak period, whereas in odd- 
The birth of new-born infants at Mt. Lasha numbered years, the sub-peak fell in the second peak 
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Figure 1 Values of I/F and I of Rhinopithecus bieti at Mt. Lasha from 2009 to 2013 


A: 2013; B: 2012; C: 2011; D: 2010; E: 2009. 


Table 2 Indexes of birth seasonality of Rhinopithecus bieti at Mt. Lasha 


T/F 


I/F 





Year Breeding period (day) F I I/F Mn Md SD (day) R Z 
2009 2-21—4-03 (42) 24 11 0.46 3-26 3-21 11.58 0.974 10.44” 
2010 2-15—4-04 (49) 27 16 0.59 3-23 3-31 9.48 0.938 14.08” 
2011 2-19—4-12 (53) 29 13 0.45 3-24 3-17 7.23 0.960 11.98" 
2012 2-17—4-06 (50) 33 15 0.45 3-28 3-28 7.27 0.966 14.00" 
2013 2-20—4-08 (48) 34 13 0.38 3-24 3-19 6.90 0.961 12.00" 
2009-10 2-15—-4-04 (49) - 27 7 3-24 3-27 6.95 0.953 24.51” 
2009-11 2-15—4-12 (57) = 40 = 3-25 3-22 6.72 0.955 36.47” 
2009-12 2-15—4-12 (57) - 55 = 3-25 3-24 6.73 0.958 50.46" 
2009-13 2-15—4-12 (57) - 68 - 3-25 3-22 6.86 0.958 62.43" 
Odd-numbered year 2-19-4-12 (53) = 37 = 3-25 3-19 5.81 0.964 34.37" 
Even-numbered year 2-15-4-06 (51) = 31 = 3-25 3-29 8.19 0.952 28.07" 


I: Infants; F: Female adults; Mn: Mean birthdate; Md: Median birthdate; SD: Standard deviation; ™s P<0.001. 
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Figure 2 Birth patterns of Rhinopithecus bieti at Mt. Lasha from 2009 to 2013 


Disconnected points mean no data was available or fewer infants were recorded than three days earlier. 


period and in even-numbered years, the sub-peak fell in 
the first peak period (Figure 2). 


Birth pattern comparison of the northern and south- 
ern monkey groups 


The Xiaochangdu and Mt. Lasha monkey groups 
were 3° apart in altitude. With increasing altitude, the 
beginning phase, ending phase and Md of the monkeys 
shifted to an earlier date; moreover, the breeding season 
shortened and birth synchrony was strengthened (SD 
decreased by ten days) (Table 3). The birth pattern of R. 
bieti at Mt. Lasha was triple-pulse-type, whereas, 


that of the Xiaochangdu group was double-pulse-type 
(peak: February 27—March 7, sub-peak: February 6—11) 
(Figure 3). However, no differences were found in the 
percentage of new-born infants between birth peaks of 
the two groups (Xiaochangdu vs. Mt. Lasha: 62.2% vs. 
50%, proportion test: Z=1.20, P=0.11). 


DISCUSSION 


Regulation of birth seasonality 
Reproduction is strongly influenced by biological 
rhythms and external stimulation (e.g., photoperiod, 


Table 3 Birth indexes of Rhinopithecus bieti at Xiaochangdu and Mt. Lasha 





Populations Breeding season (day) Mn Md SD (day) Altitude (m) Latitude Reference 
Xiaochangdu 2-04-3-13 (38) 2-28 2-24 6.50 3 500-4 250 N29°15’ Xiang & Sayers, 2009 
Mt. Lasha 2-15—4-12 (57) 3-25 3-22 16.86 2 850-3 800 N26°20’ This study 
Mn: Mean birthdate; Md: Median birthdate; SD: Standard deviation. 
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Figure 3 Accumulative numbers of Rhinopithecus bieti infants at Xiaochangdu and Mt. Lasha. Data on Xiaochangdu group was 
collected in 2005 and data on Mt. Lasha group was from 2009 to 2012 
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temperature, precipitation and food abundance) (Amroso 
& Marshall, 1960; Lancaster & Lee, 1965). Energy 
requirements and allocation have been considered the 
ultimate regulation factors in reproductive patterns. 
Other factors may play roles in regulating birth seaso- 
nality, including temperature and food abundance, to 
meet the requirements of infant survival, nursing and 
development (Brockman & van Schaik, 2005; Janson & 
Verdolin, 2005). Reproduction in Xiaochangdu and Mt. 
Lasha groups showed strong seasonality. However, 
because Xiaochangdu group’s habitat was 3° higher in 
altitude than that of the Mt. Lasha group, its reproduction 
was more concentrated over a shorter breeding season 
with an earlier beginning date (Table 3). 

More obvious seasonality was found in the climate 


and food resources in Xiaochangdu area (Table 4). The 
mating period of both Xiaochangdu and Mt. Lasha 
groups was from July to October (Xiang & Sayers, 
2009; Huang et al, 2012; Wang et al, 2012; Li, 
unpublished). The breeding season of black-and-white 
snub-nosed monkeys began after midsummer with 
shortening daytime, high temperature, abundant food 
and they started mating in July. This indicated that 
photoperiod initiated estrus and mating. Huang et al 
(2012) found that low temperature, especially the 
lowest temperature during midsummer, triggered repro- 
ductive behavior in the Mt. Lasha group. Cozzilino et al 
(1992) also found that temperature and its variations 
initiated the beginning of estrus and mating in Macaca 
fuscata. 


Table 4 Seasonality of food resources in the Xiaochangdu and Mt. Lasha monkey groups 


Xiaochangdu (Xiang et al, 2007; Xiang & Sayer, 2009) 


Mt. Lasha (Huang et al, 2012; Li, unpublished) 





Climate Food 


Climate Food 





Maximum precipitation from June Usnea thrives, other vegetation buds 


to August; in mid-May; 


Buds, leaves, flowers and fruits/nuts 
Monthly 


(10.2-12.5 °C); 


average temperature 


are in abundance from June to 


August; 


Monthly average lowest 


temperature (3.7-5.4 °C) 


Leaves fall in mid-October 


Photoperiod, temperature and precipitation affect 
food abundance, and, in turn, the reproduction period 
(Bronson, 1988; Sadleier, 1969). Based on birth date 
(Xiang & Sayer, 2009) and lactation period (195-204 
days) of black-and-white snub-nosed monkeys (Ji et al, 
1998), the conception period of Xiaochangdu group was 
calculated to be from July 24 to September 2 (median 
August 13), whereas, the conception period of the Mt. 
Lasha group was from August 6 to September 30 
(median September 9). R. bieti at Xiaochangdu began to 
conceive after one month of food abundance peak, and 
monkeys at Mt. Lasha did after two months of food peak. 
Although the conception of Xiaochangdu and Mt. Lasha 
groups began at different time points, they all terminated 
before the end of the food summit. Moreover, the median 
date of conception of the two monkey groups occurred in 
the middle of the last month of the peak in food 
availability. The long cold winter with limited food 
consumed the energy stores of the monkeys. Breeding 
individuals recovered, restored energy and completed 
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conception during the peak in food availability. 
Therefore, complete conception during the peak in food 
availability played an important role in regulating 
reproductive seasonality, and food abundance instead of 
photoperiod initiated ovulation and conception (Huang et 
al, 2012; this study). 

The decrease in both food abundance and tempe- 
rature during the conception period means that energy 
requirements increased. Infants were born during a 
period of low temperature and relatively limited food. R. 
bieti infants at Mt. Lasha began to wean at 13—14 months 
old (March-April of the following year) and ended at 18- 
old (June of the 
unpublished). This “bottle phase” (weaning) fell in the 


months following year) (Li, 
period with increased growth in “soft” and high nutrition 
foods (buds, leaves and flowers). Compared with the Mt. 
Lasha area, the peak in food availability at Xiaochangdu 
was one month shorter, buds formed 2.5 months later, 
and leaves fell one month earlier. Although we lacked 
weaning information on R. bieti infants at Xiaochangdu, 
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according to the relatively short food supply and the fact 
that infants at Mt. Lasha rarely foraged independently at 
six-months-old (Li, unpublished), we assumed that the 
chances of Xiaochangdu infants being weaned during the 
same year were quite low, and more likely occurred in 
late-May to early-June of the following year when abun- 
dant young vegetation was available. If no appropriate 
food is available during weaning, reproduction will 
ultimately fail and lead to the demise of an entire popula- 
tion, or perhaps even species. Thus, the weaning of R. 
bieti during a period of abundant high nutrition food is 
not only a prerequisite for successful reproduction, but 
also an adaptive strategy to environmental seasonality. 
Early- and late-lactation period of this species 
occurred during relatively low temperature and limited 
food periods, whereas, highly energy-consumptive mid- 
lactation overlapped with the peak in food availability. 
Monkeys stored energy during pre-conception and mid- 
lactation to alleviate energy requirements during early- 
and mid-to-late-lactation. Compared with the Mt. Lasha 
monkeys, Xiaochangdu infants were born one month 
earlier and weaned 2.5 months later. Their lactation 
period was prolonged with relatively lower temperatures, 
and higher energy consumption phase (mid-lactation) fell 
during the second peak of food resources. We concluded 
that the high energy requirement of lactation may have 
induced the high infant mortality, which was 54% at 
Wuyapuya (N28°35’) (Kirkpatrick et al, 1998), whereas 
that of the Mt. Lasha group was much more lower (6.7% 
in 2012) (Li, unpublished). To avoid energy overdraft of 
adult females during lactation (especially during winter) 
and increase infant survival, high frequency of male- 
infant-caretaking was found in Xiaochangdu group 
(Xiang et al, 2010), whereas such behaviors were only 
observed at low frequency in the Mt. Lasha group (Wang 
et al, 2012). The short breeding season with early start in 
the Xiaochangdu group appears to be an adaptive stra- 
tegy to high environmental seasonality. Research on M. 
fuscata groups at higher altitude showed that infants born 
early but with a prolonged developmental period were 
able to survive the low temperature and lack of food 
experienced during through winter (Food & Aimi, 2003). 
The reproduction of many high altitude species is 
affected by photoperiod, with the length of daytime and 
its variation triggering related neuroendocrine mechan- 
isms (Van Horn, 1980; Bronson & Heideman, 1994; 
Goldman, 2001). The habitat of Xiaochangdu monkeys 
was located at a higher altitude than that of Mt. Lasha 
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monkeys; therefore, more obvious seasonality in both 
climate and food was shown. Although there were diffe- 
rences in birth patterns of two monkey groups, their birth 
peaks stabilized during the following 6 to 7 months after 
the peak in food availability. Thus, the relative stability 
of the breeding season seems to be an adaptation to 
seasonal environments, which was critical for its succe- 
ssful reproduction (Table 3). Previous research has sho- 
wn that 100% of wild black-and-white snub-nosed mon- 
keys and 96.3% of captive individuals start reproduction 
after mid-summer (Cui et al, 2006), indicating that phot- 
operiod is the triggering factor of their estrus and mating 
behaviors. Although the stable food supply decreased the 
seasonality and synchrony of infants born, a birth peak 
was still observed in captive Rhinopithecus monkeys (R. 
bieti: Cui et al, 2006; R. roxellana: Zhang et al, 2000). 
These findings indicate the existence of boilogical 
reproductive rhythms that are not influenced by either 
photoperiod or temperature (Fooden & Aimi, 2003). This 
study indicated that conception and mid-lactation of R. 
bieti overlapped with the peak in staple food availability, 
and weaning overlapped with the abundance of high 
quality foods (new suitable vegetation). Therefore, food 
availability was the ultimate factor determining the 
reproductive seasonality of R. bieti, while other enviro- 
nmental factors, such as precipitation, temperature and 
photoperiod, functioned as proximate factors for fine- 
turning coordination between seasonal breeding and food 
availability (Wikleski et al, 2000). However, whether 
conception or weaning has priority in this regulation of 
reproductive seasonality still needs long-term, systematic 
observations of monkeys in different habitats. 


Birth synchrony 

Restrictions in environmental conditions induce 
seasonal synchrony (Power, 2013). When reproduction 
of Northern Plains Gray Langur (Semnopithecus entellus) 
is not tied with seasonality, synchrony is lost. Obvious 
reproductive signals (such as menstruation) in females 
enable males to control mating, and thereafter decrease 
food competition and increase their fitness. However, the 
more restricted reproduction is by seasonality, the more 
significant the reproductive synchrony. Females may 
conceal reproductive signals to confuse paternity and or 
decrease the chance of adult males initiating infanticide 
2013). In 
seasonality results in greater synchronicity in female 
ovulation and birth (Power, 2013). R. bieti births at 


(Power, addition, strong reproductive 


Volume 35 Issue 6 


482 LI, et al. 


Xiaochangdu were more synchronized than those at Mt. 
Lasha (SD was 6.5 and 16.9 days, respectively). This 
indicated that higher habitat altitudes resulted in stronger 
birth synchronicity, and may also be correlated with 
highly synchronized ovulations in female monkeys. 
Reproductive synchrony of females, including 
estrus, birth and length of breeding season, prevents 
males from monopolizing mating (Carnegie et al, 2011). 
The more females share reproductive synchrony, the 
more difficult it is for males to control the reproduction 
process (Power, 2013). Breeding patterns, such as 
polygamy and multi-male-female groups, are determined 
by reproductive seasonality because even the strongest 
male is unable to totally monopolize multiple females in 
synchronized pregnable states (Srivastava & Dunbar, 
1996). As a result, reproductive synchrony is prone to 
induce monkey groups with multiple males and females, 
whereas polygenous groups are more common with 
reproductive non-synchrony, which is beneficial in 
ensuring paternity (Power, 2013). The Xiaochangdu 
monkeys showed a double peak in births, while the 
monkey at Mt. Lasha displayed three peaks. These 
phenomena resulted from the strong reproductive 
synchrony (brief breeding season) of the high altitude 
monkey groups, the rather large group size (Xiaochangdu 
group was 210 in 2003 and the Mt. Lasha group was 130 


in 2012), and the birth synchrony among different OMUs. 


Annual differences in birth patterns 

The birth peak of the Mt. Lasha group was “V” 
shaped in even-numbered years, but exhibited an 
increasing trend in odd-numbered years. In odd- 
numbered years, the Md shifted ten days earlier than that 
in even-numbered years. Moreover, the commencement 
and termination of the breeding season occurred earlier. 
These phenomena may be correlated with annual 
differences in food resources or strategies adopted by 
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OMUs to decrease food competition; however, more 
evidence from long-term food detection and systematic 
research is needed. 


Decrease in reproductive rate 

Adult female R. bieti at Mt. Lasha gradually 
increased; however, no such trend was found in infants. 
From 2009-2011, all females successfully reproduced. In 
2012, there were 33 adult females (including four 
females that reached maturation in 2012) in the group, 
but there were only 15 infants. In 2013, the numbers of 
adult females and infants were 34 and 13, respectively. 
These results indicate that either some adult females did 
not reproduce (e.g., in 2012, five females did not 
reproduce and in 2013, six did not), or some infants died. 
Further research is required to determine if this low 
reproductive rate is caused by individual aging, female 
sexual competition, habitat saturation, or by infant deaths. 


Breeding season determination 

Black-and-white snub-nosed monkeys inhabit 
original coniferous forests at high altitude in temperate 
regions. Due to the steep terrain and harsh climates of 
their habitats and the nature of monkeys to avoid humans, 
it is difficult to perform close observation or achieve 
individual recognition. Under such circumstances, only 
partial I/F data was acquired with one observation, and 
only I/F can be used to estimate breeding season (Xiang 
& Sayers, 2009). However, if all adult females and 
infants could be recognized in multiple observations, 
then the I-value alone is sufficient to infer accurate 
breeding seasonality and birth pattern (this study, Huang 
et al, 2012). 
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Abstract: Coxsackie virus A16 (CA16) is commonly recognized as one of the main human pathogens of hand-foot-mouth disease 


(HFMD). The clinical manifestations of HFMD include vesicles of hand, foot and mouth in young children and severe inflammatory 
CNS lesions. In this study, experimentally CA16 infected tree shrews (Tupaia belangeri) were used to investigate CA16 pathogenesis. 
The results showed that both the body temperature and the percentages of blood neutrophilic granulocytes / monocytes of CA16 
infected tree shrews increased at 4—7 days post infection. Dynamic distributions of CA16 in different tissues and stools were found at 
different infection stages. Moreover, the pathological changes in CNS and other organs were also observed. These findings indicate 
that tree shrews can be used as a viable animal model to study CA16 infection. 


Keywords: Coxsackie virus A16; Infection; Tree shrew 


Hand-foot-mouth disease (HFMD) is a pathogenesis 
that mainly affects infantile populations. Although most 
infection cases show only mild symptoms (Mao et al, 
2013), there are a small number of individuals suffering 
from neurological symptoms of encephalitis and severe 
heart and lung failure (Tan et al, 2014). From the view of 
the pathogen, HFMD is mainly caused by enterovirus 
type 71 (EV71) and Coxsackie virus group A typel6 
(CA16) infection (Yang et al, 2014). Currently, EV71 
vaccine has completed phase III clinical trial (Chen et al, 
2014; Li et al, 2014; Zhu et al, 2014), whereas a CA16 
vaccine is still under development. Due to the lack of 
systematic analysis on the infection and immunization of 
CA16 virus, an effective animal model for CA16 
infection is urgently necessary. The conventional animal 
models of CA16 infection, such as mice or rats (Liu et al, 
2014), can neither fully reflect the etiological or 
pathological characteristics of infection and immuni- 
zation, nor meet the comprehensive assessment require- 
ments of vaccine safety and efficacy. Non-human prim- 
ate animal models are of special importance in vaccine 
evaluation, e.g., Zhang et al (2014) used infant rhesus 
monkeys (Macaca mulatta) to explore the efficacy and 
safety of EV71 vaccine. However, high cost and limited 
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animal resources has significantly limited its application. 
Due to a variety of unique characteristics, e.g., small 
adult body size, short reproductive cycle and life span 
(Wang et al, 2013), and especially a close affinity to 
primates, tree shrews (Tupaia belangeri) have been widely 
applied in biomedical research, especially in virology 
research (Fan et al, 2013). In this study, 2-month-old tree 
shrews were used to investigate GX18 infection and their 
potential as an animal model for CA16 infection was 
assessed. 


MATERIALS AND METHODS 


Virus and cells 

The CA16 virus used in this study was an isolated 
GX18 strain from an infected child with mild clinical 
symptoms in Guilin, China, 2012. The GX18 strain was 
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identified as a CA16 B2 sub-genotype virus by RT-PCR. 
The virus grown in Vero cells (cultured by Institute of 
Medical Biology) was harvested for freezing at —20 °C 
with a concentration of 10°° cell culture infectious doses 
(CCIDs59)/mL. The Vero cells were maintained in 5% 
MEM containing fetal bovine serum (provided by 
Institute of Medical Biology) and grown to a confluent 
monolayer in a T-25 vessel or in 96-well plates for viral 
isolation and antibody neutralization assays. 


Virus titration and neutralization assay 

CA16 harvested from cell culture or isolated from 
different organs and tissues of infected tree shrews were 
analyzed by a microtitration assay using a standard 
protocol (Arita et al, 2006). The neutralization test was 
performed according to the standard protocol (WHO, 
1988). A mixture of diluted serum containing anti-CA16 
antibodies (provided by Institute of Medical Biology) 
and the virus at a titer of 500-1 000 CCIDso in 100 uL of 
PBS (provided by Institute of Medical Biology) was 
incubated at 37 °C for 1 h. The cellular pathogenic effect 
(CPE) of the virus was examined by inoculating the 
mixture onto Vero cells grown in 96-well plates. The 
plaque assay was performed in 6-well plates containing 
90% confluent Vero monolayer cells. The cells were 
inoculated with samples from blood, throat swabs or feces. 
After 1 h incubation at 37 °C, the infection media were 
aspirated off and each well was then covered with 3 mL of 
agar overlay medium and incubated for 2 days at 37 °C in 
a CO, incubator. At the end of the incubation, the cells 
were fixed with formaldehyde and stained with 0.1% 
Crystal Violet as a standard protocol (Hung et al, 2010). 


Tree shrews 
Fifteen healthy female tree shrews (80+10 g, 2- 
month-old) were divided into a infection group (n=12) 





and a mock-infection control group (n=3). All animal 
procedures were approved by the Office of Laboratory 
Animal Management of Yunnan Province, China. 
Animals were individually caged and fed according to 
the guidelines of the Committee of Experimental Anim- 
als at the Institute of Medical Biology, Chinese Academy 
of Medical Sciences (CAMS). Individuals were conf- 
irmed free of antibodies against CA16 prior to the 
experiment via neutralization test. 


CA16 infection 
Our previous study on EV71 infection in rhesus 
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monkeys showed that virus distributions in blood and 
other organs can be observed via intravenous and 
intratracheal inoculation routes (Zhang et al, 2011). Ong 
et al (2008) found that infections via intravenous route in 
cynomolgus monkeys (Macaca fascicularis) can not yield a 
similar distribution of the virus. 

In this study, 12 tree shrews were infected with 
CA16 (10*° CCID50) via the respiratory tract by nasal 
spraying. The other 3 tree shrews were mock sprayed by 
saline as uninfected negative controls. The animals were 
monitored daily, and their body temperatures were 
measured rectally by a digital stick thermometer (MC- 
BOMR, Omron Co.) from day 2 post infection (p.i.). In 
parallel, venous blood samples were taken into EDTA 
coated capillary tubes daily for routine detections of 
biological indicatorsand viral loads (Veterinary Multi- 
species Hematology System, Hemavet 950FS, Drew 
Scientific Co.). Pathogenic and histopathological examinat- 
ions were performed on all the tissues and organs of tree 
shrews sacrificed by anesthesia at day 4, 7, 10 and 14 p.i.. 


Viral RNA extraction and quantitative RT-PCR 

Viral RNA was extracted from whole blood (100 uL), 
fresh tissue homogenate (10%, 100 mg) or fecal homo- 
genate (5%, 100 mg) by TRNzol-A‘, according to 
manufacturer’s protocols (Tiangen, China). The viral 
RNA was eluted in a final volume of 20 uL. For 
quantification, a single-tube, real-time Taqman RT-PCR 
assay was performed using the Taqman l-step RT-PCR 
Master Mix in the 7500 Fast Real-time RT-PCR system 
(Applied Biosystems, Foster City, CA, U.S.). The 
reaction mixture (20 uL) contained the Taqman Univ- 
ersal PCR master mix, primers (each at 20 pmol/L), the 
FAM/TAMRA labeled probe (10 pmol/L) (Takara Co., 
Ltd. Dalian, China), and RNA (2 uL). The sequences of 
the CA16-specific primers and probe were as following: 
forward primer (5'-ACACTCCATTACCCTGAGGGT 
GTA-3'); probe (5'-ATGAGAATCAAACACGTCAG 
GGCATGGAT-3'); and reverse primer (5'-ACACTCC 
ATTACCCTGAGGGTGTA-3'). The following protocol 
was used for all PCR assays: 5 min at 42 °C and 10 s at 
95 °C, followed by 40 cycles at 95 °C for 5 s and 60 °C 
for 30 s. The standard reference curve was obtained by 
the measurement of the serially diluted virus RNA 
generated by in vitro transcription from a DNA construct 
that contains the Vpl region. Viral copies were 
quantified according to vitro-synthesized RNA by 
spectrophotometry quantification, and the quantity was 
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expressed as a relative copy number, determined by the 
equation: [(ug of RNA/uL)/(molecular weight)] xAvog- 
adro's number= viral copy number/uL. 


Histopathological examination and 
chemical analysis 


The tissue samples from different organs were fixed 
in 10% formalin in PBS, dehydrated in ethanol gradients 
and embedded in paraffin before obtaining 4 um sections 


immunohisto- 


for further H-E staining. Histopathological analysis of 
the tissue sections from each organ was performed under 
a light microscope. For the immunohistochemical analy- 
sis, tissue samples were embedded in an optimal cutting 
temperature (OCT) compound (Miles Inc., Elkhart, Ind.) 
and frozen in liquid nitrogen. The frozen tissues were 
then cut into 4 um sections, placed on poly-l-lysine- 
coated glass slides and fixed in 3.7% paraformaldehyde. 
The endogenous peroxidase activity of the tissues was 
inhibited by treating with hydrogen peroxide (2.5%). The 
CA16 antigen was detected by mouse anti-CA16 sera 
and horseradish peroxidase (HRP)-conjugated anti- 
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mouse IgG antibodies (Sigma, Deisenhofen, Germany) 
followed by color development with diaminobenzidine 
for the detection of the antigen-antibody reaction. 


Statistical Analysis 
All the data were expressed as the mean of three 
samples from all experiments (mean + SE). 


RESULTS 


Clinical observations of CA16-infected tree shrews 

No typical papules or vesicles on limbs or mouths 
were found in infected tree shrews, but increased body 
temperatures p.i. (Figure 1A), as well as increased 
percentages of neutrophilic granulocytes and monocytes 
in some animals (Figure 1B, C) from day 5 to 7 p.i. with 
the decrease of lymphocyte (Figure 1D) were observed. 
These findings suggest that similar clinical CA16 
infection symptoms as humans (Lo et al, 2011), such as 
fever and inflammatory responses (reflected in blood cell 
analysis), can be found in tree shrews. 
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Figure 1 Clinical symptoms of CA16 infected tree shrews via respiratory route 
A: Body temperatures of CA16-infected tree shrews; B: The percentage of neutrophilic granulocytes in leukocytes of CA16-infected tree shrews; C: The 


percentage of monocytes in leukocytes of CA16-infected tree shrews; D: The percentage of lymphocytes in leukocytes of CA16-infected tree shrews. 


Dynamic profiles of virus in blood and stool of CA16- 
infected tree shrews 


We analyzed viral load dynamic profiles in blood 
and stool at different stages of infection. The peak levels 
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of viral load were observed in the blood (500 copies/100 
uL) at day 6 p.i. and in stool (2 000 copies/100 mg) at 
day 7 p.i., respectively. These findings indicate that the 
viremia of the infected tree shrews develops from day 3 
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to day 7 p.i., whereas, the virus shedding in feces is 
detectable from day 3 to day 14 p.i. (Figure 2B). These 
virus dynamic profiles imply that the clinical symptoms 
of CA16-infected tree shrews are correlated with the 
viremia manifestations. 


Viral distributions in tissues from CA16-infected tree 
shrews 

The target organs of CA16 infection in humans 
involve several different tissues (Wang et al, 2004). To 
understand the CA16 distributions in tree shrews, the 
viral loads in various tissues were determined at different 
stages of infection. CA16 virus was detectable in the 
lymph node, lung, spleen, kidney, thigh muscle and brain 
from day 7 p.i.. Relatively high viral loads were found in 
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the spleen (356 copies/100 mg), cerebellum (290 
copies/100 mg), lung (294 copies/ 100 mg), lymph 
node (200 copies/100 mg), parotid gland (155 copies/100 
mg) and thigh muscles (290 copies/ 100 mg), whereas no 
detectable viral loads were observed in the same tissues 
of non-infected control tree shrews. In general, viral 
replication peaked on day 6 or 7 p.i. in the lymph nodes 
and glands (Figure 3B, C), whereas, peaked on day 10 p.i. 
in the CNS and spleen, lung and muscles (Figure 3A, D), 
which might be the major target organs of infection in 
tree shrews. 

It was reported that in certain cases, CA16 is able to 
target CNS and lead to viral encephalitis (Chang et al, 
1999; Goto et al, 2009). In this present study, the total 
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Figure 2 Dynamic distributions of CA16 in the blood and stool of infected tree shrews via respiratory route 


A: Viral RNA extracted from blood specimens; B: Viral load in feces. 
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Figure 3 Viral distributions in tissues of CA16-infected tree shrews 
A: Viral load in the major organs of infected tree shrews; B: Viral load in glands of infected tree shrews; C: Viral load in the lymph nodes of infected tree 


shrews; D: Viral load in CNS of infected rhesus monkeys. 
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distributions of viral loads increased from day 4 to day 
10 p.i. (Figure 3B), and higher loads were found 
specifically in the cerebellum (291 copies/100 mg) and 
lobus frontalis (294 copies/100 mg) at day 10 p.i., which 
were consistent with the observation that viral loads 
peaked in muscles and lungs on day 10 p.i.. 


Viral antigen detection in tissues from CA16-infected 
tree shrews 

To investigate CA16 replications, antigen expres- 
sions of CA16 in the related tissues were examined 
immunohistochemically. High expressions of viral 
antigen in the heart (Figure 4A), spleen (Figure 4B) and 
lung (Figure 4C) were observed. Viral antigens were also 
observed in the CNS including cerebellum (Figure 4d). 
These findings are consistent with the patterns of viral 


load distributions of widespread infections. 


Figure 4 Viral antigen expressions in tissues from CA16 
infected tree shrews 
A: Sample of the heart collected on day 10 p.i.; B: Sample of the spleen 


collected on day 10 p.i.; C: Samples of the lung collected on day 10 p.i.; 
D: Samples of the cerebellum collected on day 10 p.i.; Arrows indicate the 


stained CA16 antigen; Images are shown at 200x magnification. 


Pathological changes in tissues from CA16-infected 
tree shrews 

In this study, the pathogenic changes of CA16 
infected tree shrews were evaluated by systematic 
pathologic analysis on various tissues, including the liver, 
kidney, spleen, heart, lung, muscles, spinal cord and 
brain. Pathological responses, such as cell damage and 
inflammatory cell infiltration (Figure 5A—D) were found 
in the lung, kidney and muscles, correlating with the 
presence of high viral loads. The observed neuropath- 
ological lesions in the cerebellum indicate that the 
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Figure 5 Pathological analyses of CA16-infected tree shrews 


The typical features of pathological changes as infiltration of inflammatory 
cells (black arrow) in the lung (A), liver (B), kidney (C) and muscles (D) of 
CA16 infected tree shrews. Images show the proliferations of neuroglial and 
neuronal lesions in the thalamus (e) andspinal cord (f) of CA16 infected tree 


shrews. Images are shown at 200x magnification. 
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pathogenic changes in CNS may contribute to the degen- 
eration of neural cells, infiltration and inflammatory cell 
aggregation (Figure 5E, F). Moreover, as expected, all 
the observed pathological changes were correlated with 
the viral antigen expressions. These results demonstrate 
that CA16 has a tropism to the CNS and lung and 
thereafter may induce severe lesions of these tissues. 


DISCUSSION 


The establishment of an viable animal model plays 
critical roles in understanding infectious diseases, such 
as HFMD, and vaccine evaluations. Animal models of 
EV71 infection include suckling mice and neonatal 
rhesus monkeys (Wang et al, 2004; Zhang et al, 2011). 
The conventional animal model in the studies of CA16 
infection is still neonatal mice. The mouse model 
provides evidence of the protective immunity of CA16 
vaccine (Dong et al, 2011) but fails to provide convi- 
ncing data on the infection pathology and the compreh- 
ensive evaluation of vaccine. According to previous reports, 
tree shrews can be widely used as a potential animal model 
for many infectious diseases, such as hepatitis C, measles, 
atypical pneumonia (SARS), tuberculosis, etc (Han et al, 
2011; Xu et al, 2013; Yang et al, 2013). 

In the present study, experimentally CA16-infected 
tree shrews have been used to investigate the etiology 
and pathology of HFMD. In the CA16-infected tree 
shrews, although no typical vesicles or papules on limbs 
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Abstract: Endosymbionts influence many aspects of their hosts’ health conditions, including physiology, development, immunity, 
metabolism, etc. Tree shrews (Tupaia belangeri chinensis) have attracted increasing attention in modeling human diseases and 
therapeutic responses due to their close relationship with primates. To clarify the situation of symbiotic bacteria from their body 
surface, oral cavity, and anus, 12 wild and 12 the third generation of captive tree shrews were examined. Based on morphological and 
cultural characteristics, physiological and biochemical tests, as well as the 16S rDNA full sequence analysis, 12 bacteria strains were 
isolated and identified from the wild tree shrews: body surface: Bacillus subtilis (detection rate 42%), Pseudomonas aeruginosa 
(25%), Staphlococcus aureus (33%), S. Epidermidis (75%), Micrococcus luteus (25%), Kurthia gibsonii (17%); oral cavity: Neisseria 
mucosa (58%), Streptococcus pneumonia (17%); anus: Enterococcus faecalis (17%), Lactococus lactis (33%), Escherichia coli 
(92%), Salmonella typhosa (17%); whereas, four were indentified from the third generation captive tree shrews: body surface: S. 
epidermidis (75%); oral cavity: N.mucosa (67%); anus: L. lactis (33%), E. coli (100%). These results indicate that S. epidermidis, N. 
mucosa, L. lactis and E. coli were major bacteria in tree shrews, whereas, S. aureus, M. luteus, K. gibsonii, E. faecalis and S. typhosa 
were species-specific flora. This study facilitates the future use of tree shrews as a standard experimental animal and improves our 
understanding of the relationship between endosymbionts and their hosts. 


Keywords: Tree shrew; Microbial; Separation; Identification 


Symbiotic bacteria are bacteria living in symbiosis with 
their hosts. They influence almost every aspect of the 
physiological prosess of the host, including the growth 
and development, physiology and biochemistry, gene 
expression, metabolism, immunology, etc. Recent years, 
tree shrews (Tupaia belangeri chinensis) have been 
rapidly and widely applied in biomedical research as a 
novel experimental animal model, especially in studies 
on human diseases (Wang et al, 2010; Huang et al, 2013; 
Xu et al, 2013). Domestic tree shrews of China include 
northern tree shrews (Tupaia belangeri) and six 
suspecies (Simpson, 1945). Currently, laboratory tree 
shrews are mainly from the field, and therefore their 
unclarified genetic background, physical condition and 
situations of symbiotic bacteria have brought many 
difficulties into their application in research into the 
development of novel antimicrobial medications, the 
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bacterial diseaseas, etc. Although scientists have carried 
out studies on the symbiotic bacteria of healthy wild tress 
shrews (Gao et al, 2009; Wang et al, 1987; Wang et al, 
2011; Xing et al, 2012; Zhang et al, 2009), there is still 
no report on symbiotic bacteria from the body surface, 
oral cavity and anus of captive tree shrews, especially the 
third generation captive tree shrews. The most commonly 
used method in studying symbiotic bacteria is the meta- 
genomic methods basing on the second generation sequ- 
encing tools, which means, the extracted DNA of flora 
will connect with T-vector before Sanger sequencing. 
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In this study, the isolated symbiotic bacteria from the 
body surface, oral cavity and anus of wild and third 
generation captive tree shrews were identified based on 
their morphological, cultural, physiological and bioche- 
mical characteristics, as well as the results of 16S rDNA 
whole sequence analysis. These findings not only pro- 
vide basic data in setting microorganism standard of 
laboratory tree shrews but also are metagenomic refe- 
rence of studies on the symbiotic bacteria in tree shews. 


MATERIALS AND METHODS 


Experimental animals 

Adult (12-month-old, body weight=130—150 g), 
healthy wild (n=12, six males and six females) and the 
third generation captive tree shrews (n=12, six males 
and six females) were provided by the Laboratory 
Animal Center of Kunming Institute of Zoology, CAS. 
Wild tree shrews were captured from the western suburb 
of Kunming, Yunnan, China, and their symbiotic bacteria 
were sampled on the same day of capture. The rearing 
environment of captive tree shrews was well ventilated, at 
16-25 °C in temperature, 40%—80% in humidity, 12 h/ 
12 h (0800h—2000h lights on) in light/dark cycle. Tree 
shrews were fed a grain premixture (main ingredients 
include corn meal, wheat flour, fish meal, milk powder, 
bone meal, sugar, salt, yeast, dregs of beans and 
vitamins), fruits (apples and bananas), terzebrio molitors 
and cooked eggs. The laborotary animal production and 
utilization numbers are SCXK (Yunnan) K2012-0001 
and SYXK (Yunnan) K2012-0003, respectively. 


Experimental reagents and facilities 
Columbia blood agar bases, Salmonella Shigella agar 
bases and Mac Conkey agar bases were all from OXIOD 
(UK); Biochemical identification tubes were from Biome- 
rieux (France); Luria Bertani (LB) agar bases were homemade. 
The instruments and facitities used in this study incl- 
ude PCR instrument (Biorad Mycycler), CO. incubator (The- 
rmo scientific Series 8000DH), thermostatical rocking plate 
(COS-211C), sterile operation platform (BAKER A2), auto- 
matic autoclave (THA-3560C), transmission electron micr- 


oscope (HITACHI, H-7650), spectrophotometer (mode 721). 


Experimental methods 

Symbiotic bacteria from the body surface, oral 
cavity and rectum of wild and third generation captive 
tree shrews were sampled, cultured and isolated. The 
isolated bacteria stains were identified and systematically 
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classified based on their morphological, cultural, 
physiological and biochemical characteristics, combining 
the results of 16S rDNA whole sequence analysis. The 
flowchart of experimental procedures (approve number: 


SYBW20110416-1) is shown in Figure 1. 


Sampling under sterile 
conditions 


Enrichment culture Purification 


| 


Pure culture 


















16S rDNA sequence 
determination 





PCR amplification 








Gram-staining microscopic 
examination of purely 
cultured bacteria 






Sequence analysis 


Species/subspecies 
identification 


Figure 1 Flowchart of experimental procedures 

















Physiological and biochemical 
characteristics determing 








Bacteria sampling and pure culture 

Next to the alcohol burner inside of the sterile 
operation platform, symbiotic bacteria from the body 
surface (sampling a bit of fur by a sterile biceps), (using 
sterile cotton swab) and rectum (using sterile cotton 
swab) of 24 tree shrews were sampled. Body surface and 
oral specimens were inoculated onto the Columbia blood 
agar bases and LB agar bases, whereas, rectal specimens 
were inoculated onto the SS agar bases and Mac Conkey 
agar bases. Bacterial species can be initially and roughly 
identified according to their growth situations on differ- 
ent agar bases. Innoculated agar bases were cultured in 
37 °C incubator. The morphology, color, hemolysis and 
pigments of the colonies were observed 18-72 h later. 
Single colony (numbered from hs1 to hs12) was picked 
by a sterile toothpick and placed on a fresh surface for 
futher indentification. 


Morphological identification 

Gram-staining used in this study was ammonium 
oxalate crystal violet staining (Katznelson et al, 1964). 
Specimens were fixed in 2.5% glutaraldehyde, stained 
with 0.1mol/L tungstophosphoric acid and then observed 
under a transmission electron microscope. 


Biochemical characteristics identification 

Main biochemical indexes of isolated bacteria strain 
were determined via the biochemical identification tubes. 
Strains cultured overnight at 30 'C were inoculated into 
the biochemical identification tube via the sterile 
inoculating loop.Every strain inoculation were triplicated. 
Negative control tubes were void of bacteria. Results 
were read within 12-72 h. Positive results went through 
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a three-generation continuous inoculation. Indentification 
conclusions were obtained by comparing the data against 
the characterics of other close bacteria strains. 


Physiological characteristics identification 

The physiological characters of bacteria strains 
mainly include growth temperature, pH tolerance range 
and optimal pH, as well as salt endurance. 

Fresh bacteria fluid (5%) was inoculated onto the 
improved LB liquid agar bases, shaking cultured by 
thermostatical rocking plates (180 r/min) at 5, 10, 25, 30, 35, 
37, 40, 42, 45 and 50 °C, respectively. ODgo9 light 
obsorption values were obtained at 24, 48, 72 h, respectively, 
to determine the growth range at different temperatures. 

To obtain the optimal pH value, fresh bacteria fluid 
(5%) were inoculated onto the improved LB liquid agar 
bases with pH at 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0 and 
10.0, respectively, and shaking cultured at 30 ‘C, 180 
p/min. OD¢oo light absorption values were obtained at 24, 
48, 72 h, respectively, to determine the pH tolerance range. 

Fresh bacteria fluid were inoculated onto the improved 
LB liquid agar bases with different NaCl concentrations 
(0%-10%) and cultured at 30 °C for 48 h, then their OD¢00 
light absorption values were obtained to determine their salt 
endurance range and the optimal salt concentration. 


16S rDNA sequence analysis 
The unindentified bacteria strains were under pure 
cultures and small scale amplifications. Universal primers 
were designed based on the high conservativeness of 165 
rDNA sequence in prokaryotes and produced by Sangon 
Biotech (Shanghai, China) (forward: 5'-AGAGTTTGA 
TCCTGGCTCAG-3'; reverse: 5'-AAG GAGGTGATCC 
AAGCCGCA-3'). Bateria DNA were templates and the 
PCR products were larger than 1 400 bp. The reaction 
solution (50 uL) included 2.5 ng templated DNA, 10 uL 
forward primer, 10 uL reverse primer, 0.25 uL 
TaKaRa (5 U/uL), 5.0 uL 10x PCR Buffer (Mg” Free), 
4.0 uL MgCl, dNTP mixture and 35.75 uL sterile 
distilled water. PCR reactions were pre-denaturation at 
94 °C for 4 min, denaturation at 94 ‘C for 1 min, 
renaturation at 52 °C for 1 min, extension at 72 C for 2 
min, incubate at 72 °C for 10 min after 30 cycles and 
then kept at 8 °C. PCR products were detected by 1% 
agarose gel electrophoresis and then sequenced by 
Genscript (Nanjing, China). 
16S rRNA sequences were obtained via Chromas 
and DNASTAR.Lasergene.v7.1, and then went through 
the phylogenetic analysis via Blast (NCBI) (http: 
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www.ncbi. nlm.nih.gov/Blast/). 


Bacteria identification 

The isolated bacteria stains were identified and syst- 
ematically classified based on their morphological, cultu- 
ral, physiological and biochemical characteristics, combi- 
ning the results of 16S rDNA whole sequence analysis. 


RESULTS 


Identification results of the isolated symbiotic strains 
The results of staining microscopy, physiological 
characteristics and 16S rDNA whole sequence analysis 
of the isolated symbiotic bacteria were shown in Table 1. 
High homologies were found by comparing the 16S 
rDNA sequences of these isolated strains with the gene 
sequences of the online registed stains (http: //www.ncbi. 
nlm.nih.gov). When the results of phylogenetic tree show 
that an isolated unknown strain has more than 99% 
homology with a known strain, then according to its 
morphology, physiological and biochemical characte- 
ristics, this isolated unknown strain can be identified. 


Physiological and biochemical characteristics of the 
isolated bacterial strains 

The physiological and biochemical characteristics 
of the isolated bacterial strains, from hs1 to hs12, were 
shown in Table 2 to Table 13, respectively. 


Detected bacteria from the wild (Table 14) and third 
generation captive tree shrews (Table 15) 

In this study, 12 bacteria species were identified from 
the wild tree shrews, body surface: B. subtilis (42%), P. 
aeruginosa (25%), S. aureus (33%), S. epidermidis 
(75%), M. luteus (25%) and K. gibsonii (17%); oral 
cavity: N. mucosa (58%), S. pneumoniae, (17%); anus: E. 
faecalis (17%), L. lactis (33%), E. coli (92%) and S. 
typhosa (17%). Four bacteria specises were identified 
from the third generation captive tree shrews, body 
surface: S. epidermidis (75%); oral cavity: N. mucosa 
(67%); rectum: L. lactis (33%) and E. coli (100%). S. 
epidermidis, N. mucosa, L. lactis and E. coli were the 
major symbiotic bacteria strains of tree shews, whereas, S. 
aureus, M. luteus, K. gibsonii, E. faecalis and S. typhosa 
are species-specific bacteria strains of tree shrews. 


DISCUSSION 


Symbiotic bacteria affect many aspects of the host’s 
physiological conditions, including the growth and 
development, gene expression, metabolism, biochemical 
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Table 1 Results of staining microscopy, physiological characteristics and electrophoregram of the isolated symbiotic bacteria 
pH 





R : 6 
Microscopy Temperature (°C) tolerence range NaCl agar base (%) 
Bacteria z Bid a r a 7 Electrophoregram 
ram’s roo rowt x rowt i 
dian Morphology body Flagellum Capsule range Optimal range Optimal 





Rhabditiform; paired 
B. subtilis + orin chains; round Y 5-50 25-30 6.5-9.5 0.6-1.2 1.0 
or squre ends 


Brevibacterium; 

E. coli = median sized, blunt. N 5-40 20-30 6.5-9.5 0.1-1.0 0.6 nl 
ends; scattered or z 
paired 

E. faecalis + Globular or oval N Y N 10-45 37 4.2-4.6 0.1-6.5 4.0 


Thallus slender and 
in various sizes; 










i 7 — shape in rod or N Y 25-42 25-30 4.2-4.6 0.1-6.5 4.0 
aeruginosa ery Boa ; 
linear; paired or in 
short chains 
K. gibsonii + Blunt ends; scattered N Y 25-42 25-30 4.2-4.6 0-6.0 4.0 
Small thallua, 
N.mucosa —  reniformed; N Y 25-38 30 7.0-7.5 0-6.0 4.0 
in dual rank 
2000bp 
S.epidermidis + N N N 20-45 37 6.0-9.0 10-100 5.0 ai 
500bp 
100bp 
D am en $ Slobular or sliemly Z N N 25-40 37 45-98 10-85 3.0 some 
epidermidis oval; aciniformed see 
L. lactis + Globular or oval N N 5-40 30 6.5-9.8 0.1-7.5 2.0 
Cellular spherical; 
M. luteus + paired, in quadruplet N 5—40 25-37 5.5-9.0 0.1-7.5 25 
or clustered 
Š g ‘Globular oroval; N N N 5-40 37 65-85 0.1-8.5 4.0 
pneumoniae paired or in chains 
S. typhosa =, © Brevibacterium. N N 5-40 35 65-95 O1-10 06 od 


blunt ends; scattered 
sovep — 


100 —| 


+: Positive in Gram’s stain; —: Negative in Gram’s stain; Y: Presence; N: Absence. 
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Table 2 Physiological and biochemical characteristics of Table 5 Physiological and biochemical characteristics of 
hs1 hs4 
Characteristics Results Characteristics Results Characteristics Results Characteristics Results 
Catalase + Nitrate-reducing te Maltose z Glutin + 
Carbamide > Amylohydrolysis F Mannite = Mushroom sugar = 
Gelatin liquefaction + Citrate + Sucrose = Xylose = 
Sucrose + D-Xylose + Galactose = Acetamide ak 
D-mannitol + D-Mannose tf Carbamide + H2S = 
Sorbierite = Maltose = Arabinose = Nitrate-reducing F 
D-raffinose = L-pectinose + Fructose = Oxidase F 
+: Positive or capable of utilizing; —: negative or incapable of utilizing; +/—: Indole = Arginine hydrolysis + 
Subtle differences between strains of one genus. Glucose; acid-producing + Citrate ES 


Glucose; aerogenesis 5 
Table 3 Physiological and biochemical characteristics of 








+: Positive or capable of utilizing; — : Negative or incapable of utilizing. 








hs2 
Claases Rawls Ohtani Resulis Table 6 Physiological u a characteristics of 
S 
Mannite $ V-P test = 
Characteristics Results Characteristics Results 

Lactose + Phaseomannite = 
: Arabinose = Mannite = 

Esculiw F D+ (—) glycogen = 
V-P test = Seminose z 

Maltose + Salicin = 
Esculin hydrolysis 5 Indole = 

Glucose + Methyl red + 
Rhamnose = Fucose = 

Raffinose + Mushroom sugar + 
Aaina re Peroxidase y Lactose 3 A-methyl-D-glucoside z 
Citrate ra Tyrosine ts Methyl red = Nitrate-reducing = 
Xylose + Indole + Fructose = Acetate S 
Carbamide = Melibiose + Sucrose - Sorbose at 
Phenylalanine = Oxidase = Raffinose = Phaseomannite = 

Casein protein 
H2S = eae = Maltose = B-galactosidase = 
Fructose + Amylohydrolysis + Galactose 2 Melibiose a 
Sorbierite + Gelatin liquefaction = Sorbierite = Glycerol = 
+: Positive or capable of utilizing; —: Negative or incapable of utilizing. Glucose; acid-producing Citrat 

= itrate = 


Glucose; aerogenesis 
Table 4 Physiological and biochemical characteristics of 




















Urease = Gelatin liquefaction = 
hs3 . 
Amylase 7 Oxidase = 
Characteristics Results Characteristics Results f 
Melampyrite = Catalase + 
Lact +  Hemolyti + 
ene aan +: Positive or capable of utilizing; —: Negative or incapable of utilizing. 
Maltose + Catalase S 
Mannite + Sorbierite n Table 7 Physiological and biochemical characteristics of 
; hs6 
Glucose te 10% gail ik 
Raffinose —~ 40% gail Ye Characteristics Results Characteristics Results 
Xylose 5  {9alicit a Sucrose + Nitrate + 
Arabinose —  B-Galactose glucoside enzyme + Glucose Y Nitrite a 
Rhamnose — Oxidase z Peroxidase + Oxidase + 
Amylase — Arginine hydrolysis T Maltose: acid-producing + Indole = 
Esculin hydrolysis + — Sucrose + Asrogenesis 
Nitrate-reducing + Lactose = 
+: Positive or capable of utilizing; —: Negative or incapable of utilizing. +: Positive or capable of utilizing; —: Negative or incapable of utilizing. 
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Table 8 Physiological and biochemical characteristics of 


Table 12 Physiological and biochemical characteristics of 














hs7 hs11 

Characteristics Results Characteristics Results Characteristics Results Characteristics Results 
Robiocina + Sucrose + Starch + Synanthrin = 
Plasma-coagulase + lactose + Mannite mi Lactose a 
Oxidase + maltose + 50% gall + Maltose + 
Carbamide aH glucose + Sucrose + Arabinose + 
Methyl red + arginine F Glucose z Sorbierite $ 
Mannite + V-P test Weak + Raffinose = 
Gelatin : 7 he SE A BG : , Bie, 

: : + nitrate + +: Positive or capable of utilizing; —: Negative or incapable of utilizing. 
Liquefaction 


+: Positive or capable of utilizing; —: Negative or incapable of utilizing. 


Table 9 Physiological and biochemical characteristics of 











hs8 
Characteristics Results Characteristics Results 
Gelatin liquefaction + Glucose F 
V-P test + Mannite a 
Urease = Sucrose + 
Plasma-coagulase = Lactose + 
H2S = Maltose + 
M. R + Fructose + 


+: Positive or capable of utilizing; — : Negative or incapable of utilizing. 


Table 10 Physiological and biochemical characteristics of 








hs9 
Characteristics Results Characteristics Results 
Arginine hydrolase + Maltose + 
Glucose a Melibiose = 
Amylase z Galactose + 
Indole = Lactose F 
Methyl red + Ribose + 
Oxidase a Fructose + 
Glutin = Melizitose = 
Catalase m Raffinose x 


+: Positive or capable of utilizing; — : Negative or incapable of utilizing. 


Table 11 Physiological and biochemical characteristics of 








hs10 
Characteristics Results Characteristics Results 
Mannite = Lactose = 
Gelatin hydrolysate + Inorganic nitrogen agar + 
Glycerol = Glucose = 
Oxidase + Arginine - 
Catalase + Nitrate = 
Esculin hydrolysis + Citrate = 


+: Positive or capable of utilizing; —: Negative or incapable of utilizing. 
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Table 13 Physiological and biochemical characteristics of 











hs12 
Characteristics Results Characteristics Results 

Arabinose — Sucrose = 
Ornithine , Mannite: acid-producing + 

Aerogenesis + 
Maltose +  Carbamide 7 
D-tartrate = 
Lactose —  A-methyl-D-glucoside = 
V-P test —  Melampyrite = 
Sorbierite — Methyl red + 
Glucose: acid- + 
producing Phenylalanine = 
Glucose: aerogenesis + 
Fructose —  Glutin = 
Indole — Citrate + 
+: Positive or capable of utilizing; —: Negative or incapable of utilizing. 


changes and immunity. Therefore, it is important to 
understand the classifications of symbiotic bacteria. 


The symbiotic bacteria hosted in tree shrews are 
affected by the living environment of microorganisms 

The symbiotic bacteria hosted in tree shrews are 
diversified and complicated, and are greatly affected by 
the living environment of microorganisms. In this study, 
12 bacteria strains were isolated and identified from the 
wild tree shrews. B. subtilis widely exsits in earth and 
soiled organisms, and multiplies quickly in hay infusion. 
In this study, B. subtilis was identified only from the 
body surface of wild tree shrews, with detectable rate at 
50%, indicating it is a commonly hosted bacterium in 
wild tree shrews. P. aeruginosa, which was indentified 
from wild tree shrews with detectable rate at 25% in this 
study, is one of the most commonly found bacteria in 
earth. S. aureus is highly pathogenic and exsits in air, 
water, dirt and excretions, whereas, S. epidermidis breeds 
on body surfaces and is a normal flora. The detectable 
rates of these two Staphlococcus bacteris, which were 
33% and 75%, respectively, are probably due to their 
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Table 14 Detected bacteria of 12 wild tree shrews 
R Sampling area Numbers of the detected bacteria (n) Detectable rate 
Genus Species 

Body surface Oral Rectal Q 3 (%) 
Bacillus B. subtilis e o o 3 2 42 
Escherichia E. coli o o e 6 5 92 
Enterococcus E. faecalis o o ° 1 1 17 
Pseudomonas P. aeruginosa e o o 1 2 25 
Kurthia K. gibsonii e o o 1 1 17 
Kurthia N. mucosa o e o 3 4 58 
Staphlococcus S. aureus e o o 2 2 33 
Staphlococcus S. epidermidis e o G 5 4 75 
Lactococus L. lactis o o e 2 2 33 
Micrococcus M. luteus e o 1 2 25 
Streptococcus S. pneumoniae o e o 1 1 17 
Salmonella S. typhosa o o e 1 1 17 


e: Positive; o: Negative. 


Table 15 Detected bacteria of 12 third generation captive tree shrews 


Sampling area 


Numbers of the detected bacteria (n) 








Genus Species Detectable rate (%) 
Body surface Oral Rectal Q 3 
Escherichia E. coli o o e 6 6 100 
Neisseria N. mucosa o e o 4 4 67 
Staphlococcus S. epidermidis e o o 5 4 75 
Lactococus L. lactis o o e 2 2 33 


e: Positive; o: Negative. 


wide distrbutions in natural environments. M. luteus can 
be found in air, earth, water, as well as on the body 
surface of animals and plants. This conditioned pathogen 
may cause local or severe infections in tissues. K. 
gibsonii aften grows in animal excretions or meat 
products and so far, there are no reports of K. gibsonii as 
a pathogen. In this study, 17% K. gibsonii were detected 
from the body surface of wild tree shrews. The detectable 
rate of oral N. mucosa was 67%. It is mainly hosted in 
the oral mucosa of mammals and is a non-pathogenic 
bacterium. S. pneumoniae exists in natural environments, 
animal excretions and the pharynx nasalis of healthy 
human beings and was detected in two wild tree shrews. 
E. faecalis was only found in wild tree shrew rectums. L. 
lactis can be found in diary and plantation products. Our 
results show that its detectable rates in both wild and 
third generation tree shrew were both 33%, indicating its 
probiotic advantages in tree shrews. E. coli are widely 
distributed in natural environments and mostly are non- 
panthogenic. Although the detectable rates of Escherichia 
coli were extremely high in this study (92%—100%), we 
consider them as intestinal tract normal flora, because all 
the experimental tree shrews were healthy and did not 
show any clinical disorders. S. typhosa is one of the 
bacteria that can easily cause infections and exists in 
almost every natural environment, including the air, 
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water, food, vegetable, dirt and animal excretions. We 
only found Salmonella typhosa in two wild tree shrews, 
indicating it is a rare strain in tree shrews. 


The influences of artificial environments to the sym- 
biotic bacteria hosted in tree shrews 

The succeful artifical reproductions of tree shrews 
require several aspects the artifical environments, 
including the facilities, nutritions, as well as the proper 
managing and reproducing methods (Li et al, 2009; Jiang 
et al, 2011; Zhao et al, 2013). The long-term and frequent 
application of sterilization measures determined that less 
bacteria exist in artifical environments than in wild. 

In this study, four bacteria strains were isolated and 
identified from the third generation of captive tree shre- 
ws, suggesting that the living environments are signif- 
icantly correlated with the species of microorganisms. 

As a novel animal to study human diseases and from 
the angle of evolution, to understand the similarities and 
differences of symbiotic bacteria from the same boby 
part of human and tree shrewsis necessary. Now, Human 
Microbiome Project (HMP) has carried out detailed 
investigations on the bacteria flora of human skin, oral 
cavities and intestinal tracts to determine the factors 
influencing the growth of symbiotic bacteria, such as 
gender, race, geography distribution, diet and body weight. 
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CONCLUSIONS 


The tree shrews used in this study were all in good 
healthy condition at time of sampling and 30 days after 
sampling. Therefore, we assume that the indentified 
bacteria are major parasitic bacteria of tree shrews. The 
high detectable rate indicates that E.coli is the major 
normal parasitic bacteria, which is consistant with the 
study of Gao et al (2009). In other conventional 
laboratory animals, Staphlococcus bacteria is rare; 
however, in both wild and third generation captive tree 
shrews, its detectable rates were high (75%). Our results 
also show a high occurrence of pathogens in wild tree 
shrews (17%-33%), whereas no pathogens were found 
in the third generation captive tree shrews. 

As a novel laboratory animal, the applications of 
tree shrews are still lack of national standards (Shen et al, 
2011). For example, for the conventional animal model 
rhesus monkeys (Macaca mulatta), it is clearly pointed 
out that laboratory rhesus monkeys have to be clean of 
Salmonella spp. Pathogenic dermal fungi and Cam- 
pylobaceter jejuni. Wang et al (1987) reported that C. 
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Acoustic signal characteristic detection by neurons in ventral nucleus 
of the lateral lemniscus in mice 
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Abstract: Under free field conditions, we used single unit extracellular recording to study the detection of acoustic signals by 
neurons in the ventral nucleus of the lateral lemniscus (VNLL) in Kunming mouse (Mus musculus). The results indicate two types of 
firing patterns in VNLL neurons: onset and sustained. The first spike latency (FSL) of onset neurons was shorter than that of 
sustained neurons. With increasing sound intensity, the FSL of onset neurons remained stable and that of sustained neurons was 
shortened, indicating that onset neurons are characterized by precise timing. By comparing the values of Qio and Q30 of the frequency 
tuning curve, no differences between onset and sustained neurons were found, suggesting that firing pattern and frequency tuning are 
not correlated. Among the three types of rate-intensity function (RIF) found in VNLL neurons, the proportion of monotonic RIF is 
the largest, followed by saturated RIF, and non-monotonic RIF. The dynamic range (DR) in onset neurons was shorter than in 
sustained neurons, indicating different capabilities in intensity tuning of different firing patterns and that these differences are 
correlated with the type of RIF. Our results also show that the best frequency of VNLL neurons was negatively correlated with depth, 


supporting the view point that the VNLL has frequency topologic organization. 


Keywords: Ventral nucleus of the lateral lemniscus; Firing pattern; Frequency tuning; Intensity tuning; Mouse 


The ventral nucleus of the lateral lemniscus (VNLL) 
is the largest nucleus of the lateral lemniscus (LL) (Batra 
& Fitzpatrick, 2002). As a station in the auditory 
ascending pathway, the VNLL plays a vital role in 
acoustic signal processing by receiving afferent 
projections from the contralateral ventral cochlear 
nucleus (VCN), ipsilateral trapezoid body (TB), superior 
olivary complex (SOC) and contralateral dorsal cochlear 
nucleus (DCN), and casts efferent projections into the 
inferior colliculus (IC) (Benson & Cant, 2008; Schofield 
& Cant, 1997; Kelly et al, 2009; Huffman & Covey, 
1995; Prather, 2013). 

Electrophysiological studies on the VNLL indicate 
that VNLL neurons can detect, process and encode the 
characteristics of acoustic signals and have various 
responding reactions. The firing patterns of VNLL 
neurons include onset firing (instantaneous firing to the 
onset of acoustic signals) and sustained firing (sustained 
firing during the whole interval of acoustic signals) 
(Zhang & Kelly, 2006). Recently, it was reported that the 


firing pattern and ability of VNLL neurons to detect 
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acoustic signal characteristics are closely correlated. For 
example, the first spike latency (FSL) of onset neurons is 
short and with increasing sound intensity, the FSL shor- 
tens insignificantly, whereas those of the sustained 
neurons are opposite (Zhang & Kelly, 2006). The 
minimal threshold (MT) of onset neurons is also higher 
than that of sustained neurons (Aitkin et al, 1970). 
Regarding firing tuning in VNLL neurons, the value of 
Qio (Qio=best frequency/the band width at 10 dB above 
threshold) is used to assess the sharpening of the 
frequency tuning cures (FTC). Studies on different 
animals show that the value of Qio in onset neurons is 
smaller than in sustained neurons (Aitkin et al, 1970; 
Covey & Casseday, 1991). Regarding intensity tuning in 
VNLL neurons, the dynamic range (DR, according to the 
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rate-intensity function, RIF, of neurons, DR is usually set 
at 10% higher than the smallest firing spikes and 10% 
lower than the highest firing spikes) is used to assess 
intensity tuning abilities. The study on Dutch-belted 
rabbits (Oryctolagus cuniculus) indicates that the DR of 
onset neurons is smaller than that of sustained neurons 
(Batra & Fitzpatrick, 1999). 

We assume that the firing pattern, intensity and 
frequency of VNLL neurons have some sort of conne- 
ction. Despite being an animal model for acoustic studies, 
detection characteristics and frequency topologic organ- 
ization of VNLL neurons to acoustic signals in Kunming 
mice (Mus musculus) remain unclear. Here, using single 
unit extracellular recording, we observed the firing 
patterns of VNLL neurons and compared FSL, frequency 
and intensity tuning characteristics under different firing 
helpful 
understanding correlations between VNLL neuron firing 


patterns. Our results are in more fully 
patterns and acoustic signal detection characteristics, and 
exploring the physiological functions of VNLL neurons 


and the effects of ascending projections on the IC. 


MATERIALS AND ANIMALS 


Animal surgery and electrode preparation 

Kunming mice (7=54, 18-28 g) with normal 
hearing were obtained from the Laboratory Animal 
Center, Institute of Hubei Preventive Medicine, China. 
Animal surgical procedures are referenced in previous 
studies (Wang et al 2007; Tang et al, 2007). All experi- 
ment procedures were performed strictly in accordance 
with the guidelines published in the NIH Guide for the 
care and use of Laboratory Animals. Anesthesia was 
induced by injection of mebumal sodium (45-50 mg/kg, 
bw). After the mouse was immobilized in a stereotaxic 
apparatus, the VNLL was exposed (4.16-4.48 mm 
posterior to the bregma, 1.10-1.95 mm lateral to the 
midline, 2.85-4.10 mm in dorsal ventral depth) by 
removing skull and dura over the VNLL (Paxinos & 
Franklin, 2001). An electrode was inserted into the small 
hole (200-500 um in diameter) over the VNLL. After 
surgery, the mouse was transferred onto a shake proof 
bench in the anechoic shielded chamber and its head was 
immobilized on a homemade bracket with the aid of a tack 
(1.8 cm in length) stuck on the skull. The animal’s eyes, 
ears and the center of the speaker were at the same level. 
Background noise inside the anechoic shielded chamber 
was 29 dB SPL and the temperature was set at 25-28 °C. 
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During the experiment, 0.6% mebumal sodium solution 
was used to keep the animal anesthetized and paralyzed. 

By referencing previous studies (Cheng et al, 2013; 
Mei et al, 2012), a homemade bipolar tungsten electrode 
made of two insulated tungsten wire electrodes (FHC Inc, 
Bowdoin, ME, USA) was used to record neuronal responses 
of the VNLL and to confirm the recording site after 
experimentation. The electrode was <10 um in diameter, 
1-5 MQ in impedance and the two poles were 30-80 um 
apart. 


Acoustic signal system and neuronal response recor- 
ding 

Under free field conditions, the speaker was placed 
at 0° in the pitching position, and 60° in the contralateral 
meridional direction (opposite to the animal’s external 
auditory canal opening). Acoustic signals were produced 
by the Tucker-Davis Technology System 3 (Alachua, 
USA). Parameters, such as frequency, intensity and 
duration were controlled by computer and calibrated by a 
sound meter (2610, B&K, Denmark) and 1/4-inch 
microphone (4936, B&K, Denmark). Acoustic signals 
were tone bursts with a duration of 40 ms, rise and fall of 
5.0 ms, and a frequency of 2 stimuli/s. The bipolar 
tungsten electrode was vertically inserted into the VNLL 
via a hydraulic drive (Kopf 640, USA). Tone burst was 
adopted to locate acoustic signal sensitive neurons. The 
responding characteristics of the neurons of the dorsal 
nucleus of the lateral lemniscus (DNLL) (unpublished 
data) were used as parallel comparisons to improve 
accuracy when locating the nucleus. Neuronal responses 
were amplified by bioelectric amplifier (ISO-DAM, WPI, 
USA) via the bipolar tungsten electrode, displayed on the 
oscilloscope (PM3084, FLUKE, USA), monitored by the 
monitoring apparatus (AM9, GRASS, USA) and then 
stored in a computer. The depth, best frequency (BF) and 
MT of the acoustic signal sensitive neuron were recorded. 
Each tone burst stimulus was repeated 32 times and the 
induced neuronal firings were profiled into a peri- 
stimulus time histogram (PSTH). Under the BF, starting 
from the MT, firing times with the intensity of acoustic 
signal increasing every 10 dB segment were recorded to 
infer the RIF. Meanwhile, the highest and lowest 
frequencies of the neuron were found and their FTCs 
were inferred. 


Electrical burning and histological examination 
At the end of the experiment, the recording site in 
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the anesthetic mouse was burned for 10 s by a 0.1 mA 
direct current for precise location later. Current stimula- 
tion was generated by a stimulator and introduced to the 
mouse via a direct current stimulation isolator (Model 
SEN-7203, Nihon Kohden Co, Shinjuku, Tokyo, Japan) 
and the bipolar tungsten electrode. After a 2 h perfusion 
with 200 mL saline water and 4% paraformaldehyde, the 
brain was isolated and fixed in 4% paraformaldehyde for 
24 h. Regular paraffin sections (40 um) were prepared 
with a vibration slicing machine (VT 1000S, Leica CO, 
USA) and under Nissl’s staining (toluidine blue) for 
immunohistochemistry examination. The burned location 
on the section was observed under a light microscope 
(Leica, DM4000B, Germany) and photographed. 


Statistical analysis 

All statistical analyses were completed using SPSS 
13.0 and all figures were generated with Sigmaplot 10.0. 
Data are expressed as mean+SD. One-way ANOVA and 





t-tests were used to compare differences. 


RESULTS 


Characteristics of VNLL neurons 

For the 98 VNLL neurons we recorded in this study, 
the dorsal ventral depth was 3 3684156.2 um (2 844-4 
090 um); BF was 19.89411.12 kHz (2.6-39 kHz); MT 
was 52.66+13.11 dB SPL (13.83-71.61 dB SPL); and 
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FSL was 5.78+2.79 ms (2—-18.5 ms). Linear regression 
analysis showed linear relationships between BF and 
dorsal ventral depth (R=—0.36, P<0.001) (Figure 1A). 
We found that neurons sensitive to high frequencies are 
located in the dorsal VNLL. From dorsal to ventral, the 
BF of VNLL neurons decreased, indicating frequency 
topological organization along the dorsal-ventral-axis of 
the VNLL. No linear relationships were found among 
FSL, MT and dorsal ventral depth (R=0.1, P>0.05; 
R=0.06, P>0.05) (Figure 1B, C). 


Temporal firing pattern of VNLL neurons 

As one of the basic features of central auditory 
neurons, the temporal firing pattern reflects firing 
changes in neurons before and after stimulation. With 
acoustic signals at BF and MT+20 dB, the firing patterns 
of 98 VNLL neurons were recorded. Following the 
research of Zhang & Kelly (2006), our observed firing 
patterns were categorized into two types: onset and 
sustained. The onset type can be sub-categorized into: (1) 
phasic: only one or two spikes to the onset of the aco- 
ustic signals; (2) phasic burst: several spikes to the 
onset of the acoustic signals; and (3) double burst: fixed 
intervals can be found between the spikes to the onset of 
the acoustic signal. Likewise, the sustained type can also 
be sub-categorized into: (1) onset plus sustained: significant 
response can be found to the onset of the acoustic signals, 
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Figure 1 Correlations between recording depth and BF (A), FSL (B) and MT (C) in VNLL neurons and Nissl-stained coronal section 
of the VNLL and surrounding structures (D) 


n: Numbers of neurons; R: correlation coefficient. 
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then decreases and lasts through the entire duration; 
(2) tonic: a fixed firing rate lasts through the entire 
duration; and (3) primary-like: the firing rate to the onset 
of the acoustic signals quickly reaches maximum, then 
decreases to a certain level. The representative post- 
stimulus time histograms (PSTH) of these six firing 
patterns and neuron numbers are shown in Figure 2. 
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Among them, phasic burst had the highest percentage, 
followed by tonic, phasic and onset plus sustained, 
respectively; primary-like and double burst had the 
smallest percentages. Although differences were found 
among the percentages of various firing patterns, the 
quantities of onset (52/98, 53.06%) and sustained (46/98, 
46.94%) patterns were comparable. 
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Figure 2 Types of onset (A-1, A-2, A-3) and sustained (B-1, B-2, B-3) patterns of recorded VNLL neurons 
n: Numbers of neurons; Recording depth (um), MT (dB SPL) and BF (kHz) of these neuron are 3118, 53.83, 16.2 (A-1); 3739, 64.28, 35 (A-2); 3970, 66.79, 
10.6 (A-3); 2908, 38.14, 13.7 (B-1); 4090, 54.22, 3.4 (B-2); 3042, 3.14, 17.8 (B-3), respectively. 


Under intensity at MT+20 dB, the FSL of onset 
VNLL neurons was 5.47+1.86 ms (2-11 ms), with 
67.31% of the FSL within the range of 2.5-5 ms; the 
FSL of sustained VNLL neurons was 7.1+3.13 ms 
(2-19.5 ms) with a relatively wide distribution (Figure 
3A, C). The FSL of onset VNLL neurons was shorter 
than that of sustained neurons (P<0.001, unpaired t-test). 








When the intensity of the acoustic signal increased every 
10 dB segment to MT+50 dB, no significant changes in 
the FSL of onset VNLL neurons were observed (one-way 
ANOVA, P>0.05), whereas, the FSL of sustained 
neurons decreased (one-way ANOVA, P<0.001) (Figure 
3B, D). 
Frequency tuning of VNLL neurons 

Frequency is a basic parameter of acoustic signals. 
Auditory neurons analyze and encode external acoustic 
stimulation, which play an important role in the behavior 
of humans and other animals. FTC are commonly used in 
assessing the frequency tuning of auditory neurons. 
According to Suga (1997) and our previous data, the 


FTC of the 71 VNLL neurons recorded here were 
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divided into six types: (1) V-shape: with increasing 
intensity, the verges of high and low frequency are 
broadening; (2) U-shape: the changes in FTC are 
independent of intensities [the V- and U-shape can be 
differentiated by bandwidth level intolerance (BLI, 
BLI=BWmax/BW10), if BLI<3, it is U-shape, if BLI23, 
it is V-shape]; (3) Lower-tail-upper-sharp (LTUS): the 
low frequency verge is characterized with a long tail and 
the high frequency verge is very sharp; (4) Upper-tail- 
lower-sharp (UTLS): the high frequency verge is charac- 
terized with a long tail and the low frequency verge is 
very sharp [in LTUS and UTLS, the Slopemt=+30-mT+10230%, 
otherwise, if Slopewr+3o-mr+10<30%, the FTC is V- or U- 
shape]; (5) W-shape: two peaks can be found in FTC and 
the two frequency tuning segments are partitioned by an 
auditory completely insensitive area only under low 
intensity, but not under high intensity; and (6) Double V- 
shape: two peaks can be found in FTC and the two 
frequency tuning segments are partitioned by an auditory 
completely insensitive area under both low and high 
intensities (Figure 4). 
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Figure 3 Distribution of FSL in onset neurons (A) and sustained neurons (C) and FSL-intensity function in onset neurons (B) and 


sustained neurons (D) 


Abscissa of B and D were normalized by each neuron’s MT as the control of tone intensity; numbers above the vertical bars in B and D are the numbers of 


neurons calculated statistically. 
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Figure 4 Types of frequency tuning curves of recorded VNLL neurons 
n: Numbers of neurons; Recording depth (um), MT (dB SPL) and BF (kHz) of these neuron are 3114, 32.22, 16.8 (A); 3925, 49.63, 11.2 (B); 3001, 61.66, 25.5 


(C); 3513, 55.81, 11.5 (D); 3009, 58.03, 40.3 (E); 3179, 67.09, 51.8 (F), respectively. 


Among the 71 neurons, 38 were onset and 33 


were sustained. The MT of the onset neurons was 
55.33410.83 dB SPL (27.09-69.34 dB SPL) and that 





of the sustained neurons was 48.834 
(13.83-73.54 dB SPL) (P<0.05, 





12.85 dB SPL 


unpaired t-test) 


(Figure 5 A). The value of Qio of the FTC of onset 
neurons was 2.49+1.67 (0.65-7.01) and that of the 
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sustained neurons was 3.43+3.83 (0.8-19.26) (P>0.05, 
unpaired t-test), whereas, the value of Q3o of the FTC 





of onset neurons was 0.944 


0.52 (0.26—2.45) and that 





of the sustained neurons was 0.96+0.26 (0.13-1.99) 
(P>0.05, unpaired t-test) (Figure 5 B, C), indicating 


that the firing patterns and 
neurons are irrelevant. 
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Figure 5 Differences in the MT (A), Qio value (B) and Q3 value (C) between onset neurons and sustained neurons 


*: P<0.05. 


Intensity tuning of VNLL neurons 

Intensity reflects the amplitude of acoustic vibration. 
Within certain range, most of the auditory neurons 
intensity. Under 
relatively high intensity, different neurons show different 


increase spiking with increasing 
changes in firing rate function. With an acoustic signal at 
BF, 40 ms, and intensity increases every 10 dB segment, 
the firing numbers of VNLL neurons changed and 
showed various characteristics of intensity tuning. Accor- 
ding to previous work and our own data (Tang et al, 2007; 
Zhou & Jen, 2000), under increasing intensity, if the 
firing numbers keep increasing, neurons are recognized 
as monotonic RIF; if the firing numbers do not change or 
decrease to <25%, neurons are recognized as saturated 
RIF; and if the firing numbers decrease to=25%, neurons 
are recognized as non-monotonic RIF. The representative 
curves of these three types of RIF and numbers of 
neurons are shown in Figure 6. The monotonic RIF had 
the highest percentages, followed by saturated and then 
non-monotonic RIF. 

Of the 91 VNLL neurons we recorded in this study, 
49 were onset and 42 were sustained. The percentages of 
monotonic, saturated and non-monotonic RIF of onset 
neurons were 48.98% (24/49), 42.86% (21/49) and 





8.16% (4/49), respectively; of sustained neurons the 
percentages were 61.9% (26/42), 21.43% (9/42) and 
16.67% (7/42), respectively (Figure 7A). Therefore, in 
onset and sustained neurons, the distribution patterns of 
RIF from high to low were monotonic, saturated and 
non-monotonic, respectively. The DR of the onset 
neurons was smaller than that of sustained neurons 
(P<0.01, unpaired t-test). In neurons with monotonic and 
saturated RIF, the DR of the onset neurons was smaller 
than sustained neurons (P<0.01; P<0.05, unpaired t-test); 
in neurons with non-monotonic RIF, no differences were 
found between onset and sustained neurons (P>0.05, 
unpaired t-test) (Figure 7B). 


DISCUSSION 


Firing patterns and temporal features of VNLL 
neurons 


The percentages of onset and sustained patterns 
were comparable, 53.06% (52/98) and 46.94% (46/98), 
respectively, consistent with previous findings on big 
brown bats (Eptesicus serotinus) (Covey & Casseday, 
1991), rats (Zhang & Kelly, 2006) and rabbits (Batral & 
1999). The 2.04% (2/98) double-burst 
pattern we found has only been reported in rats (Zhang & 


Fitzpatrick, 
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Figure 6 Types of rate-intensity function of recorded VNLL neurons 
n: Numbers of neurons; Recording depth (um), MT (dB SPL) and BF (kHz) of these neuron are 3403, 50.52, 33 (A); 3259, 27.09, 24.8 (B); 3211, 38.89, 19.8 


(C), respectively. 
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Figure 7 Differences in percentage of RIF (A) and DR (B) between onset neurons and sustained neurons 


M: monotonic; Sat: saturated; NM: non-monotonic. *: P<0.05; **: P<0.01. 


Kelly, 2006) and its biological meaning in acoustic signal 
processing remains unclear. 

The FSL of onset neurons was shorter than 
sustained neurons. With the increase in intensity, the FSL 
of onset neurons remained stable (Figure 3B), indicating 
a feature of precise timing of onset neurons. This specific 
feature of onset VNLL neurons may be attributable to 
projections of octopus cells in the cochlear nucleus (CN) 
(Adams et al, 1997; Pollak et al, 2011). The octopus cell 
have large, irregularly shaped cell bodies with three to 
five thick primary dendrites emerging from only one side 
of the cell body. The sizable low voltage activated 
potassium channels on its membranes endow the cell 
with a low impedance and fast time constant, which 
bestow its projections to the VNLL with temporal 
advantages (Johnston et al, 2010; Bal & Ocrtel, 2007). 
The precise timing of onset VNLL neurons is not only 
beneficial in the encoding of acoustic signal onset, but 
also affects the process of IC neurons to temporal 
information via its ascending projections. For example, 
studies on rats (Nayagam et al, 2005), Mexican free- 
tailed bats (Tadarida brasilensis mexicana) (Xie et al, 
2007) and big brown bats (Voytenko & Galazyuk, 2008) 
show that an inhibitory postsynaptic potential (IPSP) can 
be observed ahead of the excitations of some IC neurons, 
which means an IPSP occurs before depolarization. It is 
assumed that because of the temporal advantage of 
VNLL neurons possess, the inhibitory projection from 
VNLL neurons will reach the IC ahead of the excitatory 
projection, thereafter, an early IPSP is initiated. 

The FSL of sustained VNLL neurons decreased 
with increasing intensity (Figure 3D). This phenomenon 
is not specific to sustained VNLL neurons (Zhang & 
Kelly, 2006) and has also been observed in the CN 
(Goldberg & Brownell, 1973), SOC (Sanes & Rubel, 
1988), IC (Tan et al, 2008) and auditory cortical (AC) 
(Heil, 1997). Therefore, it is assumed that FSL encodes 
the intensity of the onset of the acoustic signal. The larger 
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the intensity, the shorter the FSL, and the quicker the 
signal can be recognized. However, correlations between 
FSL and intensity are unclear. Pharmacological studies 
show that neither y -aminobutyric acid (GABA) nor glyc- 
ine (Gly) affects correlations between FSL and intensity in 
AC and IC neurons (Zhang et al, 2005), indicating that 
these correlations may be encoded by auditory neurons 
below the IC, but this requires further exploration. 


Frequency tuning and frequency topologic organiza- 
tion of VNLL neurons 

Although the percentages of six different types of 
FTC are consistent with previous studies on rats (Zhang 
& Kelly, 2006) and rabbits (Batra & Fitzpatrick, 1999), 
the degrees of sharpening are quite different among 
different species. For example, in this study, the average 
value of Qio of the FTC of VNLL neurons in Kunming 
mice was 2.9, but was 7.1 (Zhang & Kelly, 2006) in rats 
and 9.1 (Covey & Casseday, 1991) in bats, respectively, 
indicating that different species are characterized with 
different capabilities of frequency tuning. The sharp- 
ening of FTC and GABAergic inhibitory projection are 
correlated with the formation of the inhibitory area of 
frequency tuning (Chen & Jen, 2000). The FTC of VNLL 
neurons in mice have smaller Qio and low degrees of 
sharpening. Because the VNLL in mice receives less 
inhibitory projections from the ipsilateral trapezoid body, 
but has a wide frequency receptive range to excitatory 
projections from the contralateral cochlear nucleus, 
perhaps the FTC of VNLL neurons is not characterized 
by significant sharpening. 

No differences were found in the Qjy of the FTC in 
onset and sustained neurons (Figure 5). However, studies 
on big brown bats and cats show that the Qio in onset 
neurons is smaller than sustained neurons (Covey & 
Casseday, 1991; Aitkin et al, 1970). Based on these 
discrepancies, it is assumed that different species have 
different neural circuits corresponding to specific behaviors. 
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For example, because the frequency tuning of the onset 
VNLL neurons in Mexican free-tailed bats affects the reor- 
ganization of IC neurons to some special communication 
signals, they show certain corresponding firing patterns and 
features of frequency tuning (Pollak et al, 2011). 

Frequency topological organization of VNLL neu- 
rons remains controversial. When the IC of rats was 
infected with retrograde tracing agents, low frequency 
neurons in the VNLL distributed along the band, whereas, 
high frequency neurons were evenly distributed (Kelly et 
al, 1998). Electrophysiological experiments on bats 
(Covey & Casseday, 1991) and cats (Aitkin et al, 1970) 
show that low frequency neurons are distributed dorsally 
of the VNLL, whereas, high frequency neurons are 
distributed ventrally of the VNLL. Although differences 
exist in the distribution patterns of different species, the 
logical distribution of low and high frequency neurons 
supports the view of frequency topologic organization in 
the VNLL. However, an electrophysiological experiment 
on rats (Nayagam et al, 2006) and the study on infecting 
retrograde tracing agents in the IC of cats (Glendenning 
& Hutson, 1998) and the primary auditory field (AI) of 
Mongolian gerbils (Budinger et al, 2013) indicate a lack 
of topologic organization in the VNLL. Possible reasons 
of these discrepancies are differences in animal models 
and experimental methods. In this study, the BF and 
dorsal ventral depth of VNLL neurons were negatively 
correlated (Figure 1A). Neurons with high frequencies 
were located dorsally of the VNLL and from dorsal to 
ventral, the 
frequency topologic organization along the dorsal- 
ventral-axis of the VNLL. 


BF of neurons decreased, indicating 


Intensity tuning of VNLL neurons 

In the 91 VNLL neurons recorded here, the 
percentage of non-monotonic RIF (12.09%) was lower 
than that in the advanced auditory nuclei of mice, such as 
the IC (Tang et al, 2007) and AC (Qi et al, 2013). The 
occurrence of non-monotonic type is due to inhibitory 
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Abstract: Viviparidae are widely distributed around the globe, but there are considerable gaps in the taxonomic record. To date, 18 
species of the viviparid genus Cipangopaludina have been recorded in China, but there is substantial disagreement on the validity of 
this taxonomy. In this study, we described the shell and internal traits of these species to better discuss the validity of related species. 
We found that C. ampulliformis is synonym of C. lecythis, and C. wingatei is synonym of C. chinensis, while C. ampullacea and C. 
fluminalis are subspecies of C. lecythis and C. chinensis, respectively. C. dianchiensis should be paled in the genus Margarya, while 
C. menglaensis and C. yunnanensis belong to genus Mekongia. Totally, this leaves 11 species and 2 subspecies recorded in China. 
Based on whether these specimens’ spiral whorl depth was longer than aperture depth, these species or subspecies can be further 
divided into two groups, viz. chinensis group and cathayensis group, which can be determined from one another via the ratio of spiral 
depth and aperture depth, vas deferens and number of secondary branches of vas deferens. Additionally, Principal Component 
Analysis indicated that body whorl depth, shell width, aperture width and aperture length were main variables during species of 


Cipangopaludina. A key to all valid Chinese Cipangopaludina species were given. 


Keywords: Cipangopaludina; Gastropoda; Bellamyinae; Anatomy; Taxonomy; Chinese 


The freshwater gastropod Viviparidae is currently 
divided into three subfamilies (Bouchet & Rocroi, 2005): 
Viviparinae Gray, 1847, Bellamyinae Rohrbach, 1837 
and Lioplacinae Gill, 1863, which together comprise 
approximately 150 species and 31 genera recognized to 
date (Franke et al, 2007). Typically, viviparid species 
inhabit lakes, ponds, and lentic rivers of temperate to 
tropical region (Strong et al, 2008; Ying et al, 2013), 
explaining their wide global distribution across every 
continent save for South America and Antartica (Brown, 
1994). In China, the existing taxonomy for Viviparidae 
includes approximately 61 recognized species in 9 
genera (Du et al, 2011; Liu, 1991; Wang & Xie, 2005; 
Yen, 1943): Viviparus Montfort, 1810, Bellamya Jousse- 
aume, 1886, Filopaudina Habe, 1964, Mekongia Crosse 
et Fischer, 1876, Angulyagra Benson, 1836, Margarya 
Nevill, 1877, Rivularia Heude, 1890 and Trochotaia 
Brandt, 1974, Cipangopaludina Hannibal, 1912. The last 
genus, Cipangopaludina is particularly interesting as a 
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diverse genus mainly distributed in China, Japan, Korea, 
Thailand, Vietnam, Laos, India, Burma and Malaya 
(Brandt, 1974; Liu et al, 1993) that is comprising of 35 
species and subspecies (Global Names Index; http:// 
gni.Globalnames.org), there has been substantial debate 
on the taxonomy of Cipangopaludina. Hannibal (1912) 
had previously proposed Cipangopaludina as a subgenus 
of Idiopoma, based on invasive specimens collected in 
California which were identified as Paludina malleata 
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Reeve. Meanwhile, Annandale (1920) proposed Lecyth- 
oconcha as a distinct genus for Paludina lecythis Benson 
while Prashad (1928) argued that Lecythoconcha was 
invalid and adopted the name Cipangopaludina. 

Even within China, there is no clear verdict on the 
validity of each species or how many are even distributed 
in the area. Heude (1890) listed 9 species from China, 
viz. C. chinensis (Benson (non Gray)), C. lecythoides 
(Heude), C. diminuta (Heude), C. longispira (Heude), C. 
fluminalis (Heude), C. leucostoma (Heude), C. cath- 
ayensis (Heude), C. ventricosa (Heude) and C. aubryana 
(Heude). Later, comprehensive work by Kobelt (1909) 
revised Viviparidae, including practically all the forms of 
the family known at the time, which was one species and 
12 subspecies in China, with most of species being 
designated varieties of the species C. chinensis (Gray): 
cathayensis (Heude), compacta (Nevill), ventricosa 
(Heude), (Heude), fluminalis (Heude), 
hainanensis (Kobelt), diminuta (Heude), leucostoma 


longispira 


(Heude), aubryana (Heude), patris (Kobelt), wingatei 
(Smith), lecythoides (Benson) and its form latissima 
(Dautzenberg et H. Fischer). Several decades later, Yen 
(1943) revised species of Chinese Viviparidae, and 
pointed out that leucostoma and diminuta were synon- 
yms of chinensis and that lecythoides (Heude, non 
Benson) distributed in Zhoushan Island, compacta (Ko- 
belt, non Nevill) distributed in Hainan, and lecythis 
crassior (Annandale) distributed in western Yunnan were 
all synonyms of patris (Kobelt), and Jongispira (Heude) 
was argued to be synonymous with haasi. Liu (1991) 
later clarified the 
Cipangopaludina and raised 15 of the subspecies to the 


distribution and taxonomy of 


species level, viz. C. chinensis, C. cathayensis, C. 
fluminalis, C. aubryana, C. latissima, C. lecythoides, C. 
ussuriensis (Gerstfeldt), C. yunnanensis, C. dianchiensis, 
C. ampullacea, C. haasi, C. lecythis, C. longispira, C. 
ventricosa and C. menglaensis. Finally, Wang & Xie 
(2005) evaluated the extinction risk of Cipangopaludina 
species, listing 14 distinct species: Cipangopaludina 
hainanensis, C. leucostoma, C. menglaensis, C. latissima, 
C. dianchiensis, C. lecythis C. lecythoides, C. haasi, C. 
ampullacea, C. ventricosa, C. ampulliformis, C. 
aubryana, C. longispira and C. ussuriensis. Curiously, 
these latest works by Liu (1991) and Wang & Xie (2005) 
made no mention of two species, C. patris and C. 
wingatei, and ignored some of the earlier work done by 
Yen (1943). 


Despite discrepancies and arguments as well as the 
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elevation and demotion of species to subspecies and vice 
versa, the existing literature names 18 potential species 
of Cipangopaludina distributed in China: C. ampullacea, 
C. ampulliformis, C. aubryana, C. cathayensis, C. 
chinensis (synopsis C. leucostoma and C. diminuta), C. 
C fluminalis, C. haasi (synopsis C. 
longispira), C. hainanensis (synopsis C. compacta), C. 


dianchiensis, 


latissima, C. lecythoides, C. lecythis, C. menglaensis, C. 
patris, C. ussuriensis, C. ventricosa, C. wingatei and C. 
yunnanensis. Similarly, the disparate results of the 
previous studies on the topic also highlight how difficult 
and complex the identification of Cipangopaludina is, 
largely due to intraspecific variations in shell shape. This 
has left the taxonomic position and phylogenetic 
relationships of these species largely unknown. In this 
study, we opted to perform a fresh check on the 
morphology and anatomy of Chinese Cipangopaludina 
species currently being held at repositories in Beijing and 
Kunming in order to arrive at a more complete species 
list. Morphometric analyses of shell variations were 
conducted on Cipangopaludina species through a 
Principal Component Analysis. 


MATERIALS AND METHODS 


Sample collection and identification 

Four hundred and ninety-six, including 11 Cipang- 
opaludina species were collected from Jilin, Jiangxi, 
Zhejiang, Guangxi and Yunnan Provinces. These specim- 
ens were preserved in Kunming Institute of Zoology, 
Chinese Academy of Sciences (KIZ). Additionally, 145 
Cipangopaludina specimens, including 18 Cipangop- 
aludina species or its synonyms listed in the existing 
literature preserved in the Invertebrate Museum of the 
Institute of Zoology, Chinese Academy of Sciences, 
Beijing (IOZ) were measured. A morphological analysis 
was conducted on these specimens and the validity of the 
currently identified species was determined. Specimens 
were identified mainly based on original description, 
type specimens, Liu et al (1979, 1993), Kobelt (1909) 
and Prashad (1928). 


Characters selected for morphometric analysis 

To select characters for use in morphometric anal- 
ysis, we reviewed the available taxonomic descriptions 
from literature (Liu et al, 1979). Five shell characters 
were measured with callipers accurate to 0.1 mm: Shell 
depth (D), the maximum dimension parallel to the axis of 
coiling; Shell width (W) the maximum dimension 
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perpendicular to D; Length of the aperture (LA), the 
maximum dimension from the junction of the outer lip 
with the penultimate whorl to the anterior edge of the 
aperture (oblique to coiling axis); width of the aperture 
(WA), the maximum dimension perpendicular to LA; 
depth of the body whorl (BW), the dimension from the 
lower margin of the aperture to the upper suture 
delimiting the first whorl; and N, the number of whorls. 
Embryonic shells were measured to 0.1 mm using 
an ocular micrometer. Anatomy was studied using a 
microscope with drawing apparatus. Radulae and embry- 
onic shells were studied via scanning electron micro- 
scopy (SEM). The radulae were cleaned enzymatically 
with proteinase K, as described by Holznagel (1998), 
then sonicated and mounted on aluminium specimen 
stubs with adhesive pads. Embryonic shells were clea- 
ned mechanically, sonicated, and mounted on adhesive 
carbon-coated pads. Both radulae and embryonic shells 
were coated with gold-palladium and studied on AMR- 
AY 1000B scanning electron microscope at 30 kV. The 
description of shell features and visceral hump follows 
the general terminology suggested by Simone (2004). 


Allometric shell growth 

Due to ontogenetic effect of allometric growth have 
usually caused confusion in morphometric analyses, 
these individuals within the range of allometric growth 
were eventually excluded before multivariate analyses 
(Armbruster, 1995; Chiu et al, 2002; Valovirta & Vaisa- 
nen, 1986). However, in the present-day study, it is 
difficult for us to obtain snails of all size classes from 
each population per species and examine the allometric 
shell growth. Under this circumstance, the consequence 
of Chiu et al (2002) was followed that excluded 
individuals with aperture lengths of less than 1.0 cm. 
Morphometric data was explored through a Principal 
Component Analysis by SPSS 10.0 (SPSS for Windows, 
Chicago, IL, USA). 


RESULTS 


Systematics Cipangopaludina Hannibal, 1912 

Idiopoma (Cipangopaludina) Hannibal 1912, Proc- 
eedings of the Malacological Society of London 10: 194 
(type species, Paludina malleata Reeve, original desig- 
nation) 

Lecythoconcha Annandale 1920, Records of the 
Indian Museum. 19: 111 (type species, Paludina lecythis 
Benson, monotypic); Annandale 1921, Records of the 
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Indian Museum, 6: 401. 

Viviparus (Cipangopaludina), Clench & Fuller, 
1965, Occasional papers on Mollusks, Museum of Com- 
parative Zoology, 2: 385—412. 

Bellamya (Cipangopaludina), Smith, 2000, Nautilus, 
114(2): 31-37. 

Diagnosis: Shell medium to large (up to 70 mm), 
oval, brown or greenish in colour. Apex pointed or 
obtused. Suture narrow. Umbilicus, big or narrow or 
closed. Sculpture lacking except growth lines and axial 
undulations. Aperture subcircular, inner lip whitish blue 
and outer lip black in colour, easily broken. Young 
specimens (from brood pouch) with up to 4 whorls; 
protoconch smooth, three primary rows of chaetae on last 
whorl and other whorls have two rows of chaetae. Male’s 
testis with two or three main branches. 

Cipangopaludina chinensis was the only species 
given a detailed description, while diagnosis charac- 
terstics are listed for all other species. All morphometric 
and meristic characters of Cipangopaludina species are 
given in Table 1. 


Cipangopaludina ampullacea (Charpentier, 
(Figure 1A) 

Paludina ampullacea Charpentier: Küster, 1852, Syste- 
matisches Conchylien—Cabinet. 2 ed.p. 23, pl. 4, fig 2-3. 

Material: 25 ex. (2 female and 3 male were anat- 
omy). KIZ 000110-000117, collected from Qiubei Cou- 
nty, Yunnan in August 2005 by Cui GH and Du LN; KIZ 
000453-000470, collected from Beihai, Tengchong 
County, Yunnan in April 2006 by Chen XY. 

Diagnosis: shell larger, up to 53 mm, brown- 


1863) 


greenish, with 5 whorls. Apex obtused. Spiral whorl 
depth is less than aperture length. Body whorl inflated, 
shell width about 80.1%-—90.6% shell depth. Body whorl 
depth is less than shell width. Spiral whorls inflated and 
round, form a flat shoulder. Osphradium ridge-like, 
shorter or equal with length from anterior tip of 
osphradium to mantle border. The base of gill filaments 
wide and apex narrow and round, curve to left. Kidney 
arch-shape, the dorsal length of kidney right plane about 
1.5 times than ventral left plane and posterior plane, and 
left plane and posterior plane very curve (Figure 2C-II). 
Vas deferens with 2 main branches and the first anterior 
branch with 4 parallel secondary branches. Vas deferens 
opens in last 1/4 of testis. Pouch filled with less than 44— 
48 embryos in various stages of development. Embryos 
with 4 whorls, body whorl with 3 chaetae and other 
whorls have 2 rows of chaetae. 
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Figure 1 Shell of Cipangopaludina 
A: C. lecythis ampullacea, KIZ 000453; B: C. lecythis, KIZ000284; C: C. leucostoma, KIZ 000270; D: C. chinensis, IOZ—FG166277; E: C. chinensis fluminalis, 


KIZ 000443; F: C. haasi, KIZ 000101; G: C. dianchiensis, KIZ 000536; H: C. hainanensis, IOZ-FG168161; I: C. latissima, IOZ—-FG84656; J: C. lecythoides, 
IOZ-FG166538; K, L: C. menglaensis, IOZ—-FG0142 (Holotype); M: C. yunnanensis, IOZ—-FG0152 (Holotype); N: C. aubryana, IOZ-FG 170556; O: C. 


ussuriensis, IOZ—-FG84819; P: C. ventricosa, KIZ000227. Scale bars=1 cm. 
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Figure 2 Anatomy and genital system of Cipangopaludina 


A: Shell removed showing external features of animal; B: Male with mantle cavity opened mid—dorsally; C: Two type of kidney, I triple shape and II arch shape. 


Abbreviations: aa — anterior aorta; ag — albumen gland; an — anal papilla; au —auricle; bp — brood pouch; cm—columellar muscle; ct — ctendium; df — 


dorsal fold of buccal mass; eg — egg-shell gland; es — esophagus; ey — eye; ft — foot; int — intestine; jw 





jaw; kd — kidney; 1 — liver; me — mantle edge; mo 





— mouth; ne — nephrostome; op — operculum; pa — posterior aorta; pe—penis; pg—prostate gland; re—rectum; sn—snout; te—tentacle; ur—ureter; ve— 


ventricle; vd—vas deferens. 


Distribution: Southwest of China, Yunnan, Sichuan 
provinces and Macao (Liu, 1991). 

Remarks: Cipangopaludina ampullacea can be 
distinguished from C. chinensis, C. fluminalis, C. haasi, 
C. menglaensis, C. dianchiensis and C. ussuriensis by 
spiral depth shorter than aperture depth. Cipangop- 
aludina ampullacea can be distinguished from other 
species by the following characters: osphradium length 
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shorter or equal with length from anterior tip of 
osphradium to mantle border, vs. longer in C. ventricosa, 
C. patris and C. aubryana; shell larger, apex eroded vs. 
shell medium and apex pointed in C. yunnanensis; 
umbilicus small or closed vs. big and round in C. 
latissima; apex obtused vs. pointed in C. cathayensis and 
C. lecythoides. Cipangopaludina ampullacea can be 
distinguished from C. lecythis by apex eroded vs. pointed; 
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the upper angle of aperture acute vs. straight. However, 
the anatomy of C. ampullacea and C. lecythis are not 
obviously different, so it is suitable to treat C. 
ampullacea as a subspecies of C. lecythis. 


Cipangopaludina aubryana (Heude, 1890) (Figure 1N) 

Paludina aubryana Heude, 1890, Memoires conce- 
mant L’Histoire naturelle e L’empire chinois par des 
peres de la compagnie de Jesus, 175, pl. 39. 

Vivipara (chinensis) aubryana Heude, Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p 118, pl. 
20, fig 6-7. 

Viviparus chinensis aubryana (Heude), Yen, 1939, 
Diese Arbeit erscheint gleichzeitig in den Abhandlungen 
der Senckenbergischen Naturforschenden Gesellschaft 
Abhandlung, 444, p. 35, taf. 3. 

Viviparus chinensis aubryana (Heude), Yen, 1942, 
Proceedings of the Malacological Society of London, 190, 
200, pl. 13, fig 32. 

Cipangopaludina chinensis aubrayana (Heude), 
Pradshad, 1928, Memoirs of the Indian Museum, Calc- 
utta, 8: 168. 

Cipangopaludina lecythoides aubryana (Heude), 
Yen, 1943, Nautilus, 56(4): 124-130. 

Cipangopaludina aubryana (Heude), Liu, 1991, 
Proceeding of Tenth International Malacological Con- 
gress, 590. 

Materials: 20 ex. (4 females and 4 males were anat- 
omy). IOZ-FG166735, IOZ-FG166737, from Wang Co- 
unty, Sichuan Province in June 1964; IOZ-FG170556, 
IOZ-FG170553, IOZ-FG170551, IOZ-FG 170549, from 
Hanshou County, Hunan province in March 1963; KIZ 
000087—100, from Pu’er City, Yunnan in April 2011 by 
Zheng LP. 

Diagnosis: shell medium, brown green, with 5 who- 
rls. Body whorl inflated, but other whorls not, without 
shoulders. Apex pointed. Spiral whorl depth less than 
aperture depth. Osphradium ridge-like, longer than 
length from anterior tip of osphradium to mantle border. 
The base of gill filaments wide and apex narrow and 
pointed, does not curve to right. Kidney triple-shape, the 
dorsal length of kidney right plane about 1.2 times than 
ventral left plane, 2 times than posterior plane, the left 
plane and posterior plane do not curve much (Figure 2C- 
I). Vas deferens with 3 main branches and the first 
anterior branch with 4 parallel secondary branches. Vas 
deferens opens in last 1/3 of testis. Pouch filled with less 
than 10 embryos in various stages of development. 
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Embryos with 4 whorls, body whorl with 3 chaeta and 
other whorls have 2 rows of chaetae. 

Distribution: Sichuan, Hunan, Guangdong, Guizhou 
and Yunnan Provinces, China (Liu, 1991). 

Remarks: Cipangopaludina aubryana can be disti- 
nguished from C. chinensis, C. haasi, C. fluminalis, C. 
dianchiensis, C. ussuriensis and C. menglaensis by spiral 
depth less than aperture length and shell medium. 
Cipangopaludin aubryana can be distinguished from 
other species by the following characteristics: kidney 
triple-shape vs. arch-shape in C. ampullacea, C. patris 
and C. ventricosa, acute spiral whorl and without 
shoulder (vs. spiral whorl inflated and form shoulder in C. 
lecythoides and C. yunnanensis); shell medium and body 
whorl not very inflated vs. shell larger and body whorl 
very inflated in C. latissima; apex pointed, not eroded vs. 
apex eroded, with only 4 whorls in C. hainanensis, vas 
deferens opens in last 1/3 of testis vs. 1/4 in C. lecythis 
and C. l. ampullacea. 


Cipangopaludina cathayensis (Heude, 1890) (Figure 1C) 

Paludina catayensis Heude, 1890, Memoires conce- 
mant L'Histoire naturelle e L’empire chinois par des 
peres de la compagnie de Jesus, 175, pl. 39. 

Vivipara chinensis cathayensis Heude, Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p 112, pl. 
18, fig 5—6. 

Cipangopaludina chinensis cathayensis (Heude): 
Prashad, 1928, Memoirs of the Indian Museum, Calcutta, 
8: 168; Yen, 1943, Nautilus, 56(4): 124-130. 

Cipangopaludina cathayensis (Heude), Liu, 1991, 
Proceeding of Tenth International Malacological Cong- 
ress, 590. 

Materials: 34 ex. (2 males and 2 females were 
anatomy) KIZ000016-17, collected from Lake Xingyun, 
Yunnan in July 2008 by Du LN; KIZ000422-434, 
collected from Lake Dianchi, Yunnan in September 2006 
by Du LN, Jiang YE and Cui GH; KIZ 000076 collected 
from Debao County, Guangxi in August 2009 by Yu GH 
and Li Y, KIZ000190-196, collected from Liaoyuan City, 
Jilin in October 2012 by Du LN. 

Diagnosis: shell shin, medium, up to 60 mm, with 6 
inflated whorls, spiral whorl with shoulder. Apex pointed. 
Spiral whorl depth less than aperture depth. Osphradium 
length similar with length from anterior tip of osphr- 
adium to mantle border, anterior 2/3 part of osphradium 
ridge-like and posterior 1/3 part inflated. Gill filaments 
same with C. chinensis. Kidney triple-shape. Vas defer- 
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ens with 3 main branches and the first anterior branch 
with 3—4 parallel secondary branches. Vas deferens opens 
in last 1/3 of testis. Female pouch filled with more than 
60 embryos in various stages of development. Embryos 
with 3 whorls, body whorl with 3 chaeta and other 
whorls have 2 rows of chaetae. 

Distribution: Jilin, Shanxi, Hebei, Shandong, Anhui, 
Jiangsu, Zhejiang, Hunan, Sichuan, Guizhou and Yunnan 
provinces (Liu, 1991). 

Remarks: In the original description, this species 
was spelt catayensis, but later reviewers all spelt it 
cathayensis, according to the International Code of 
Zoological Nomenclature (1999). Catayensis is then 
treated as an inadvertent error, and thus this spelling is 
incorrect. Cipangopaludina  cathayensis can be 
distinguished from C. chinensis, C. haasi, C. fluminalis, 
C. dianchiensis and C. menglaensis by spiral depth less 
than aperture depth. Cipangopaludina cathayensis can be 
distinguished from other species by the following 
characteristics: anterior 2/3 of osphradium ridge-like and 
posterior 1/3 inflated; female’s pouch with more than 60 
embryos, and embryo’s have 3 whorls. 


Cipangopaludina chinensis (Gray, 1834) (Figure 1D) 

Paludina chinensis Gray: in Griffith & Pidgeon, 
1834, Animal Kingdom, pl. 1, fig 5. 

Paludina chinensis Gray: Gray in Griffith & 
Pidgeon, 1834, Animal Kingdom, 599. 

Vivipara chinensis (Gray), Kobelt, 1909, Abbildun- 
gen Nach de Natur mit Beschreibngen, p111, pl 23, fig 1. 

Cipangopaludina chinensis (Gray): Prashad, 1928, 
Memoirs of the Indian Museum, Calcutta, 8: 168. 

Viviparus chinensis chinensis (Gray): Yen, 1939, 
Diese Arbeit erscheint gleichzeitig in den Abhandlungen 
der Senckenbergischen Naturforschenden Gesellschaft 
Abhandlung, 444, 35. Taf. 3. 

Viviparus chinensis (Gray): Yen, 1942, Proceedings 
of the Malacological Society of London, 190, 200, pl. 13, 
fig 32. 

Cipangopaludina chinensis leucostoma (Heude): 
Yen, 1943, Nautilus, 56(4): 124—130. 

Cipangopaludina chinensis diminuta (Heude): Yen, 
1943, Nautilus, 56(4): 124-130. 

Bellamya (Cipangopaludina) chinensis 
Smith, 2000, Nautilus, 114(2): 31-37. 

Cipangopaludina chinensis (Heude): Prashad, 1928, 
Memoirs of the Indian Museum, Calcutta, 8: 168; Liu, 
1991, Proceeding of Tenth International Malacological 


(Gray): 
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Congress, 590. 

Material: 40 ex. (1 males and 2 females were 
anatomy). IOZ-FG166277, Shangsi Co. Guangxi Prov- 
ince in 19" April 1974; IOZ84521, Zhaotong City, Yun- 
nan Province in May 1932; I0Z166773, IO0Z166777, 
Qionghai Lake, Xichang County, Sichuan Province in 4" 
July 1964; 10Z166733, Luofang village, Guangdong 
Province in March 1982; KIZ000164—165, Lake Erhai, 
Yunnan, collected by Du LN, Jiang YE and David 
Aldridge in 14 April 2006; KIZ 000294-310, Lake north 
mountain, Jilin City, Jilin Province, collected by Du LN 
on 13 October 2012; KIZ 000132-133, KIZ 000149-150, 
West Lake, Dali City, Yunnan, collected by Du LN in 13 
April 2006, Jiang YE and Aldridge D; KIZ 000157-158, 
Lake Xinyun, Jiangchuan County, Yunnan, collected by 
Du LN, Jiang YE and Cui GH in September 2009; KIZ 
000013015, KIZ 000038040, KIZ 0000363-366, KIZ 
000421, Lake Dianchi, Kunming City, Yunnan, collected 
by Du LN and Aldridge D. in March 2006. 


Description 

Shell: large (up to 70 mm), oval, brown or greenish 
in colour, with 6—7 whorls. Apex pointed. Suture narrow. 
Umbilicus narrow or closed. Sculpture lacking except 
growth lines and axial undulations. Aperture subcircular, 
inner lip whitish blue and outer lip black in colour, easily 
brown. Young specimens (from brood pouch) with up to 
4 whorls; protoconch smooth, three primary rows of 
chaetae on last whorl and other whorls have two rows of 
chaetae. 

Operculum: Occupying entire shell aperture. Centre 
is brown and outer parts greenish. Nucleus sub-central, 
located closer to inner margin. Sculptured with conce- 
ntric growth lines. Inner surface glossy; scar sub-circular, 
close to inner margin, occupying about 2/3 of opercula 
area. 

External morphology (Figure 2B): head and foot 
black. Snout cylindrical, anterior margin flat. Length of 
tentacles about 1.3 times longer than snout length, base 
at side of snout base. Ommatophore short, located 
between basal and middle third of outer surface of each 
tentacle. Foot large, sole simple. Opercular pad larger 
than its base in dorsal foot surface. Columellar muscle 
thick. 

Radula: Central tooth with wide rounded major 
denticles and 4 smaller triangular denticles on each side. 
Lateral tooth with tongue-shaped major denticles and 3 
minor denticles on each side. Inner marginal tooth with 
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tongue-shape major denticles and 3 minor denticles and 
marginal teeth with 10 almost equal-sized denticles. 

Mantle organs (Figure 2A-B): Mantle border simple 
and thick. Osphradium ridge-like, very close to ctenid- 
ium, about same length from anterior tip of osphradium 
to mantle border. Ctenidium long and narrow, about 
same length as pallial cavity, anterior end at mantle 
border. Gill filaments very tall and narrow, arched 
towards right and apex slightly pointed, close to food 
groove. 

Hypobranchial gland lies on left of gill margin. 
Right margin of mantle cavity filled by oviduct in 
females. Ureter runs along mantle cavity right margin in 
males or edging oviduct in females. Rectum on dorsal 
and left sides of ureter. At right end of mantle border, in 
females, presenting three openings, most posterior and 
smaller is ureter pore, anus most anterior, and female 
pore larger, located between anus and ureter pore. 

Alimentary canal: Oval mouth, bounded by fleshy 
lips, ventral at anterior end of snout. Jaw plates ridge-like 
(Figure 2A). Esophagus with pair of dorsal, longitudinal 
folds, about 1.5 times of buccal mass length. Esophagus 
runs along the hypobranchial gland to posterior of 
columellar muscle, and then turns right and passes 
upwards in floor of pericardial chamber to reach 
digestive gland in upper part of visceral hump where it 
curves round to open into stomach. Stomach located 
halfway posterior to pallial cavity, posterior half 
immersed in digestive gland. Esophagus is inserted in 
posterior gastric extremity, and as gradually enlarges 
curves towards left. After this curve, stomach abruptly 
expands becoming almost as wide as whorl, posterior 
narrowing, without clear separation with intestine. 
Intestine forms loop overlying pericardial cavity and 
when it reaches digestive gland turns sharply forwards to 
continue as rectum. Rectum passes forwards on right of 
mantle cavity to open at anus lying just behind mantle 
edge. 

Circulatory and excretory systems: Heart located at 
anterior ventral of kidney. Auricle connected with 
ctenidial vein just posterior to gill. Ventricle posteriorally 
connected to auricle by a thin tube. The size of ventricle 
is same as the auricle. Kidney triangular (Figure 2 C-I), 
with four almost plane surfaces: 1) posterior surface 
towards pericardium; 2) left with pallial cavity; 3) 
ventral right with ureter; 4) dorsal with mantle. The 
dorsal length of the kidney’s left plane is similar to the 
ventral right plane, and about two times longer than 
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posterior plane. Nephrostome is a round, small pore, 
located in centre basally inferior region of ventral surface 
of kidney. Ureter runs between rectum and oviduct (in 
females) or right pallial cavity edge (in males). Ureter 
pore just posterior to anus (or just posterior to female 
pore). 

Genital system. 

Male (Figure 2A): Testis compact, semi-lunar, on 
right of mantle cavity. Extends to upper end of mantle 
cavity, where apex close to pericardial cavity also 
connected by thin fold of membrane with lower surface 
of digestive gland. Testis flattened laterally, right surface 
abuts columellar mussel while left surface abuts right 
wall of ureter. Vas deferens very narrow, running on 
columellar and inferior ventral margin of testis, with 2 
main branches, and in the anterior branch with 5-6 
parallel secondary branches, the branches with many 
very narrow branches from different portions of testis. 
Vas deferens opens in last 1/3 of testis and runs to 
prostate gland in mantle cavity along mantle cavity floor 
for about 3/4 of its length; left 1/4 length of vas deferens 
abruptly narrows and is surrounded by very thick, 
muscular walls. Vas deferens runs all along the right 
cephalic tentacle, its inner surface relatively broad, with 
some inner, longitudinal folds. Vas deferens opens in 
tentacle tip as broad papilla with rounded tip. Right 
tentacle possesses deep concavity located at right from 
genital papilla. Papilla can retract into its concavity. 

Female. Ovary grey, located in same position as 
testis. Ovary in close contact with posterior wall of 
cardiac region of stomach and along course of hepatic 
artery. Oviduct very narrow; from albumen gland it runs 
to posterior end of columella, and then turns left, making 
a loop. Oviduct increases in size and runs to capsule 
gland along brood pouch. Brood pouch located in dorsal 
of oviduct. Albumen gland tongue-shaped, slightly 
curved, capsule gland on dorsal and posterior surface and 
opens to brood pouch. Walls of brood pouch thin, semi- 
transparent, smooth. Pouch filled with about 45-84 
embryos in various stages of development. 

Ecology and habitat: it lives in slow-moving water 
such as lakes, pond, irrigation canals, ditches and slow 
moving stream (Pace, 1973). It is a benthic grazer and 
filter feeder, feeding on benthic and epiphytic diatom 
(Plinski et al, 1978), often found on sandy to muddy 
substrates. 

Distribution: original range in China, Taiwan, Korea, 
and Japan, but has also been introduced in freshwater 
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ponds and lakes in Canada, the north-eastern USA, and 
Europe (Soes et al, 2011). 

Remarks: Although Yen (1943) mentioned that C. 
leucostoma (Heude) and C. diminuta (Heude) are 
synonymous with C. chinensis (Gray), the authors 
thought these two species are synonymous with C. patris 
(Kobelt) according to apex obtuse and less than 6 whorls. 
Additionally, Cipangopaludina wingatei (Smith) could 
be synonymous with C. chinensis according to shell size, 
spiral depth longer than length of aperture, and number 
of whorls. Cipangopaludina chinensis, C. dianchiensis, 
C. fluminalis, C. haasi, C. and C. 
menglaensis share the characteristic that spiral depth is 


ussuriensis 


longer than length of the aperture. However, C. chinensis 
can be distinguished from C. dianchiensis and C. 
ussuriensis by its shell smooth rather than having a shell 
with a ridge; female’s pouch has more than 80 embryos 
vs. less than 10 embryos in C. dianchiensis; embryo with 
three primary rows of chaetae on last whorl vs. 3 ridges 
in C. dianchiensis. Cipangopaludina chinensis can be 
distinguished from C. fluminalis by the following 
characters: shell thin vs. thick, umbilicus close vs. big 
and deeply. 


Cipangopaludina dianchiensis Zhang, 1990 (Figure 1F) 

Cipangopaludina dianchiensis Zhang, 1990, Acta 
Zootaxonomica Sinica, 15(1): 25-27. 

Materials: 34 ex. (2 males and 3 females were 
anatomy). KIZ 000233, collected in November 2002 by 
Cui GH; KIZ 000001-KIZ 000004, collected in June 
2008 by Du LN, Jiang YE and Barclay H; KIZ 000005- 
KIZ 000012, collected in September 2005 by Du LN and 
Aldridge D; KIZ 000212-KIZ 000217, KIZ 000287-KIZ 
000293; KIZ 000311—KIZ 000316, collected in January 
2006 by Du LN and Yuan C; KIZ 000373, KIZ 000231 
collected in June 2006 by Du LN. All specimens were 
from Lake Dianchi. 

Diagnosis: shell thick. Body whorl with 7—12 ridges 
(3-4 ridges clearly). Apex pointed. Osphradium ridge- 
like, shorter than length from anterior tip of osphradium 
to mantle border. Gill filaments and jaw plate are the 
same as C. chinensis. The size of ventricle is the same as 
the auricle. Vas deferens has 3 main branches and the 
first anterior branch has 4 “V” shaped secondary 
branches. Vas deferens opens in last 1/3 of testis. Pouch 
filled with about 6-8 embryos in various stages of 
development. Embryos with 5 whorls, body whorl with 3 
ridges and other whorls with 1 ridge. 
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Distribution: only known to live in Lake Dianchi, 
Yunnan. 

Remark: Morphology characteristics, such as shell 
thickness, body whorl with clear ridges, female’s pouch 
with few (6-8) embryos and embryo shell with ridges 
can be distinguished clearly from other species of genus 
Cipangopaludina. Additionally, phylogenetic analysis 
indicated that C. dianchiensis clustered with species of 
Margarya in Lake Dianchi (Du et al, 2013). The morp- 
hology and molecular results indicated that C. dianch- 
iensis should be placed into the genus Margarya. 


Cipangopaludina fluminalis (Heude, 1890) (Figure 1E) 
Heude, 1890, Memoires 
concemant L’Histoire naturelle e L’empire chinois par 


Paludina fluminalis 


des peres de la compagnie de Jesus, 175, pl. 39. 

Vivipara chinensis fluminalis Heude, Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p 116, pl. 
19, fig. 4-5. 

Viviparus chinensis fluminalis (Heude), Yen, 1939, 
Diese Arbeit erscheint gleichzeitig in den Abhandlungen 
der Senckenbergischen Naturforschenden Gesellschaft 
Abhandlung, 444, p. 35, taf. 3; Yen, 1942, Proceedings of 
the Malacological Society of London, 190, 200, pl. 13, 
fig 32. 

Cipangopaludina_ chinensis fluminalis (Heude), 
Prashad, 1928, Memoirs of the Indian Museum, Calcutta, 
8: 168. 

Cipangopaludina lecythoides fluminalis (Heude): 
Yen, 1943, Nautilus, 56(4): 124—130. 

Cipangopaludina fluminalis (Heude), Liu, 1991, 
Proceeding of Tenth International Malacological 
Congress, 590. 

Materials: 21 ex (1 male and 1 female were 
anatomy). KIZ 000181 collected from Jinping County, 
Honghe prefecture, Yunnan in April, 2011 by Jiang WS; 
KIZ 000204-000211 collected from West Lake, Dali 
prefecture, Yunnan in April 2006 by Du LN, Jiang YE 
and Aldridge D; KIZ 000384 collected from Songming 
County, Kunming prefecture, Yunnan in December 2005 
by Du LN; KIZ 000442-000451 
Songming County, Kunming prefecture, Yunnan in 
December 2004 by Du LN. 

Diagnosis: shell is thick, larger, up to 60 mm, with 7 


collected from 


slightly inflated whorls. Apex pointed. Spiral whorl 
depth longer than aperture depth. Umbilicus big and deep. 
Osphradium length similar with length from anterior tip 
of osphradium to mantle border. Gill filaments same as C. 
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chinensis. Kidney triple shaped, the dorsal length of 
kidney left plane about same length as ventral right plane, 
1.5 times longer than posterior plane. Vas deferens with 3 
main branches and the first anterior branch with 5-6 
parallel main secondary branches. Vas deferens opens in 
last 1/30f testis. Pouch filled with 12 embryos in various 
stages of development. Embryos with 4 whorls, body wh- 
orl with 3 chaeta and other whorls have 2 rows of chaetae. 
Distribution: Yangtze River (Yen, 1943; Liu, 1991). 
Remarks: very 
similar with C. chinensis except the shell thick rather 


Cipangopaludina fluminalis is 


than thin; umbilicus big and deep rather than small or 
closed. The anatomy characteristics were relatively 
similar with little difference, making it is possible to treat 
C. fluminalis as subspecies of C. chinensis. 


Cipangopaludina latissima (Dautzenberg et H. Fischer, 
1905) (Figure 1I) 

Paludina lecythoides latissima Dautzenberg et H. 
Fischer, 1905, Journal de Conchylologie, p418, pl. 10, 
fig 17; Kobelt, 1909, Abbildungen Nach de Natur mit 
Beschreibngen, p 144, pl. 29, fig 1. 

Cipangopaludina lecythoides latissima (Dautzen- 
berg et H. Fischer), Pradshad, 1928, Memoirs of the 
Indian Museum, Calcutta, 8: 168. 

Cipangopaludina latissima (Dautzenber et H. Fis- 
cher): Yen, 1943, Nautilus, 56(4): 124-130; Liu et al, 
1991, Proceeding of Tenth International Malacological 
Congress, 590. 

Materials: 15 shell specimens. IOZ-FG84653-84667, 
collected from Mengzi County, Yunnan in March 1927. 

Diagnosis: shell is large, round, up to 53 mm, with 6 
whorls. Apex pointed. Body whorl very inflated, shell 
width is about 85%—95% of shell depth. Due to the body 
whorl inflated, the upper of body whorl has a flat 
shoulder. Aperture ellipse. Umbilicus very big, round. 
Anatomy characteristics unknown. 

Distribution: only known to be from Mengzi County, 
Yunnan. 

Remarks: Cipangopaludina latissima can be distin- 
guished from other species of Cipangopaludina by the 
body whorl special inflated, and the upper of body whorl 
with a flat shoulder. 


Cipangopaludina lecythoides (Benson, 1842) (Figure 1J) 

Paludina lecythoides Benson, 1842, The Annals and 
Magazine of Natural History, Zoology, Botany and Geo- 
logy 9, 488; Philippi, 1846, Abbildungen und Beschrei- 
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bungen neuer oder wenig gekannter Conchylien heraus- 
gegeben, II. p.133, pl 2, fig 1; Küster, 1852, In Abbildu- 
ngen nach der Natur mit Beschreibungen. Systematisches 
Conchylien-Cabinet 2ed, p. 23, pl. 5, fig 1, 2; Reeve, 
1862, Conchologia iconica, or, illustrations of the shell 
of molluscous animals, pl. 4, fig 21; Dautzenberg et H. 
Fischer, 1905, Journal de Conchylologie, p417-418. 

Paludina (Vivipara) lecythoides Benson: Mabille et 
Le Mesle, 1866, Journal de Conchyliologie, p. 134; 
Fischer, 1891, Catalogue et Distribution géographique 
des Mollusques terrestres, fluviatiles et marins d’une 
partie de I’Indo-Chine (Siam, Laos, Cambodge, Conchi- 
nchine, Annam, Tonkin), p. 177. 

Vivipara (chinensis var.) lecythoides (Benson), 
Kobelt, 1909, Abbildungen Nach de Natur mit Beschrei- 
bngen, p 119, pl. 23, fig 8. 

Cipangopaludina chinensis lecythoides (Benson), 
Prashad, 1928, Memoirs of the Indian Museum, Calcutta, 
8: 168. 

Cipangopaludina lecythoides (Benson), Yen, 1943, 
Nautilus, 56(4): 124-130; Liu et al, 1991, Proceeding of 
Tenth International Malacological Congress, 590. 

Materials: 1 ex. KIZ 000507, collected from Hang- 
zhou City, Zhejiang Province in April 2013 by Lu HF; 

Diagnosis: shell depth 55 mm, with 6 inflated whorls. 
Shell dark brown, with clearly growth lines. Apex pointed. 
Spiral whorl depth shorter than aperture depth. Umbilicus 
small, nearly closed. Vas deferens has 2 main branches and 
the first anterior branch with 3 parallel secondary branches. 
Vas deferens opens in last 1/3 of testis. 

Distribution: Zhejiang Province. 

Remarks: Although Wang & Xie (2005) mentioned 
that this species is distributed in Yunnan, it was not 
collected during field excursions undertaken between 
2003-2008 by authors Du LN and Chen XY. Generally, 
the validity of this species is unclear, as it is difficult to 
distinguish C. lecythoides and C. cathayensis by shell 
morphology alone. Further anatomic characteristics are 
required to determine the validity of this species. 


Cipangopaludina lecythis (Benson, 1836) (Figure 1B) 

Paludina lecythis Benson, 1836, Journal of Asiatic 
Society of Bengal, 5: 745. 

Paludina chinensis lecythis Benson: Nevill, 1885, 
Hand list of Mollusca in the Indian Museum: Gastropoda. 
Prosobranchia—Neurobranchia, p. 20. 

Paludina ampullacea Reeve (not Charpentier), 1862, 
Conchologia iconica, or, illustrations of the shell of 
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molluscous animals, 12. 
Vivipara lecythis Benson, Kobelt, 1909, Abbildungen 
Nach de Natur mit Beschreibngen, p 148, pl. 30, fig 1-2. 
Cipangopaludina chinensis lecythis (Benson): Prashad, 
1928, Memoirs of the Indian Museum, Calcutta, 8: 168. 
Cipangopaludina lecythis (Benson): Yen, 1943, 
Nautilus, 56(4): 124-130; Liu et al, 1991, Proceeding of 
Tenth International Malacological Congress, 590. 
Materials: 25 ex. (1 female and 2 males). KIZ 
000018-000030 collected from Songming County, Kun- 
ming City, Yunnan in April 2006 by Du LN and Yuan C; 
KIZ 000044 collected from Lake Dianchi, Yunnan in 
April 2006 by Du LN and Yang J; KIZ 000144 collected 
from Lake Erhai in April 2006 by Du LN, Yang YE and 
Aldridge D; KIZ 000151-156 collected from Lake 
Xingyun, Yunnan in September by Du LN, Jiang YE and 
Cui GH; KIZ 000224 collected from Lake Xingyun, 
Yunnan in April by Du LN; KIZ 000284—286 collected 
from Lake Dianchi, Yunnan in March 2006 by Du LN. 
Diagnosis: shell larger, up to 60 mm, brown gree- 
nish, with 7 whorls. Apex pointed. Spiral whorl depth sh- 
orter than aperture depth. Whorls inflated, with a shoul- 
der. Aperture ellipse, the upper of aperture form a strai- 
ght shoulder. Umbilicus small or closed. Osphradium 
ridge-like, shorter than length from anterior tip of 
osphradium to mantle border. Gill filaments same with C. 
chinensis. Vas deferens with 2 main branches and the 
first anterior branch with 4 parallel secondary branches. 
Vas deferens opens in last 1/4 of testis. Pouch filled with 
about 37 embryos in various stages of development. 
Embryos with 3—4 whorls, body whorl with 2 chaetae 
and other whorls with 1 chaetae. 
Distribution: Yunnan, China. Burma (Kobelt, 1909). 
Remarks: Cipangopaludina lecythis can be distin- 
guished from other species of this genus by its special 
aperture shape, the upper of aperture straight and form a 
flat shoulder. Kobelt (1909) list C. lecythis ampulliformis 
(Souleyet) as a subspecies of C. lecythis, which are 
distributed in both Burma and Yunnan, China. According 
to Kobelt (1909) it is difficult to distinguish between 
these two species, so we opted to treat C. lecythis ampu- 
lliformis as synonymous with C. lecythis. 


Cipangopaludina patris (Kobelt, 1906) (Figure 1C) 
Paludina patris Kobelt, 1909, Abbildungen Nach de 

Natur mit Beschreibngen, p 118, pl. 23, fig 4-5. 
Cipangopaludina chinensis patris (Kobelt), Pras- 

had, 1928, Memoirs of the Indian Museum, Calcutta, 8: 
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168. 

Cipangopaludina lecythoides (Heude, non Benson): 
Yen, 1943, Nautilus, 56(4): 124-130. 

Cipangopaludina compacta (Kobelt, non Nevill): 
Yen, 1943, Nautilus, 56(4): 124—130. 

Cipangopaludina lecythis crassior Annandale: Yen, 
1943, Nautilus, 56(4): 124-130. 

Materials: 28 ex (2 females and 2 males were 
anatomy). KIZ 000268—000279, collected from Yuanji- 
ang County, Yunnan in January 2010 by Du LN; KIZ 
000369-000372, collected from Jinning County, Yunnan 
in August 2006 by Du LN; KIZ 000391-000402, 
collected from Yuanjiang County, Yunnan in January 
2013 by Du LN. 

Diagnosis: shell medium, up to 50 mm, brownish, 
with 6 whorls. Apex obtuse. Whorls inflated, especially 
last two whorls, upper three whorls slightly inflated and 
always eroded. Spiral whorl depth shorter or equal with 
aperture depth. Aperture ellipse. Umbilicus closed. The 
shape of the kidney is similar with C. ampullacea. Osph- 
radium ridge-like, longer than length from anterior tip of 
osphradium to mantle border. Gill filaments same with C. 
chinensis. Vas deferens with 3 main branches and the 
first anterior branch with 1 “V” shape and 3 parallel sec- 
ondary branches. Vas deferens opens in last 1/3 of testis. 
The loop of oviduct exceeds the posterior margin of col- 
umellar muscle. Pouch filled with 96-109 embryos in 
various stages of development. Embryos with 3 whorls, 
body whorl with 3 chaetae and other whorls with 2 chaetae. 

Distribution: Yunnan, Guangdong and Zhejiang pro- 
vinces, China. 

Remarks: Cipangopaludina patris can be disting- 
uished from C. chinensis, C. haasi, C. dianchiensis, C. 
ussuriensis and C. menglaensis by a spiral depth shorter 
than aperture length and a medium shell. Cipang- 
opaludina patris can be distinguished from other species 
by the following characteristics: larger shell, up to 50 
mm (vs. medium, up to 26 mm in C. yunnanensis); the 
kidney arch-shape is similar with C. ampullacea, (vs. 
triple-shape in C. cathayensis, C. aubryana and C. 
lecythoides); medium shell and body whorl are not very 
inflated (vs. shell larger and body whorl very inflated in 
C. latissima); shell with 6 whorls (vs. only with 4 whorls 
in C. hainanensis); the upper of aperture without straight 
shoulder (vs. with straight shoulder in C. lecythis). 


Cipangopaludina haasi (Prashad, 1928) (Figure 1F) 
Cipangopaludina chinensis haasi Prashad, 1928, 
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Memoirs of the Indian Museum, Calcutta, 168. 
Paludina 1890, Memoires 
concemant L’Histoire naturelle e L’empire chinois par 


longispira Heude, 


des peres de la compagnie de Jesus, 175, pl. 39. 

Vivipara chinensis longispira Heude, Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p 115, pl. 
19, fig 3. 

Viviparus chinensis longispira (Heude), Yen, 1939, 
Diese Arbeit erscheint gleichzeitig in den Abhandlungen 
der Senckenbergischen Naturforschenden Gesellschaft 
Abhandlung, 444, p. 35, taf. 3. 

Viviparus chinensis longispira (Heude), Yen, 1942, 
Proceedings of the Malacological Society of London, 190, 
200, pl. 13, fig 32. 

Cipangopaludina haasi (Prashad): Yen, 1943, Naut- 
ilus, 56(4): 124-130. 

Cipangopaludina longispira (Heude), Liu, 1991, Proc- 
eeding Tenth International Malacological Congress, 590. 

Materials: 44 ex (2 females and 2 males were 
anatomy). KIZ 000197-000198, KIZ 000280-000282, 
collected from Black dragon spring, Songming Town, 
Kunming, Yunnan in December 2004 by Cui GH and 
Chen XY; KIZ 000122-131, collected from Xiayu 
Village, Zhangshu City, Jiangxi Province in February 
2006 by Yu GH; KIZ 00101-00108, collected from Lake 
Yilong, Yunnan in February 2008 by Du LN and Jiang 
YE; KIZ 000477-000497, collected from Hangzhou City, 
Zhejiang Province in April 2013 by Lu HF. 

Diagnosis: shell larger, up to 65 mm, greenish, with 
7 whorls. Apex pointed. Suture shallow. Spiral whorl 
depth higher than aperture depth. Shell width less than 
body whorl depth. Whorls slightly inflated, do not form 
shoulder. Umbilicus small. Aperture ellipse. Osphradium 
ridge-like, longer than the length from anterior tip of 
osphradium to mantle border. Gill filaments same as C. 
chinensis. Esophagus has pair of dorsal, longitudinal 
folds, about 2 times of buccal mass length. Stomach 
elliptical, spoon-shape. Heart located anterior ventral of 
kidney. Ventricle is about 2 or 3 times larger than auricle. 
Kidney triple-shape. Nephrostome round, small pore, 
located in posterior 1/5 of ventral surface of kidney. Vas 
deferens with 2 main branches and the first anterior 
branch has 8 to 9 parallel secondary branches. Vas 
deferens opens in last 1/4 of testis. Female’s pouch 
length is about 4 times longer than egg gland. Pouch 
filled with 128-188 embryos in various stages of 
development. Embryos with 4 whorls, body whorl with 3 
chaetae and other whorls with 2 chaetae. 
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Distribution: Hunan, Jiangxi, Sichuan, Zhejiang and 
Yunnan Provinces, China. 

Remarks: Cipangopaludina haasi, C. dianchiensis, 
C. chinensis, C. ussuriensis and C. menglaensis share the 
characteristic that spiral depth is longer than length of the 
aperture. However, C. haasi could be distinguished from C. 
dianchiensis and C. ussuriensis by shell smooth vs. shell 
with ridge; female’s pouch with more than 128-188 
embryos vs. less than 10 embryos in C. dianchiensis; 
embryo with three primary rows of chaetae on last whorl 
vs. 3 ridges in C. dianchiensis. Cipangopaludina haasi 
could be distinguished from C. chinensis by the following 
characteristics: body whorl not inflated vs. inflated; body 
whorl depth same or shorter than body width vs. obviously 
longer. Cipangopaludina haasi could be distinguished 
from C. menglaensis by a larger shell, up to 65 mm vs. 
smaller, up to 30 mm; body whorl not inflated vs. inflated. 


Cipangopaludina hainanensis (Mollendorff, 1909) 
(Figure 1H) 

Vivipara (chinensis var.) hainanensis Méllendorff, 
in Kobelt, 1909, Abbildungen Nach de Natur mit 
Beschreibngen, p 117, pl. 19, fig. 6-7. 

Cipangopaludina chinensis hainanensis (Möllend- 
orff): Prashad, 1928, Memoirs of the Indian Museum, 
Calcutta, 168. 

Cipangopaludina hainanensis (MOllendorff): Yen, 
1943, Nautilus, 56(4): 124-130. 

Materials: 5 ex. (without specimens were anatomy). 
IOZ-FG 168660 collected from Ledong County, Hainan 
Island in December 1989; IOZ-FG 168160-168163 
collected from Chang Jiang County, Hainan Island in 
December 1989. 

Diagnosis: shell medium, up to 40mm. Apex eroded, 
left 4-5 inflated whorls. Suture deeply. Umbilicus big 
and deep. 

Distribution: only known from Hainan Island. 

Remarks: this species was evaluated as being 
extinct by Wang and Xie (2005). Existing specimens can 
be distinguished from other species by medium shell, 
apex eroded, with 4—5 inflated whorls and umbilicus big 
and deep. 


Cipangopaludina menglaensis Zhang, Liu et Wang, 
1981 (Figure 1L) 

Cipangopaludina menglaensis Zhang, Liu et Wang, 
1981, Acta Zootaxonomica Sinica, 6(1): 40. 

Materials: 10 ex. IOZ-FG 0142 (holotype), IOZ— 
FG 00143-151 (Paratypes) from Manzhazai Village, 
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Pu’er City, Yunnan in March 1965 by Zhang WZ and Liu 
YY. 

Diagnosis: shell dark brown, small, up to 30 mm, 
with 6 whorls. Apex pointed. Suture deeply. Central 
tooth with wide rounded major denticles and 3 or 4 
smaller triangular denticles on each side. Lateral tooth 
with tongue-shaped major denticles and 4 or 5 minor 
denticles on each side. Inner marginal tooth with tongue- 
shape major denticles and 3 or 4 minor denticles and 
marginal teeth with 20 almost equal sized denticles. 

Distribution: only known from type locality. 

Remarks: Without anatomized specimens, the anat- 
omy characteristics remain unknown. However, this 
species can be distinguished from other species of 
Cipangopaludina by a smaller shell, up to 30 mm and 
marginal teeth with 20 almost equal sized denticles. 


Cipangopaludina ussuriensis (Gerstfeldt, 1859) (Figure 
10) 

Paludina ussuriensis Gerstfeldt, 1859, Memoires 
des savants Strangers, 9: 507, fig 1-4; Reeve, 1862, 
Conchologia iconica, or, illustrations of the shell of 
molluscous animals, pl. 8, fig 1—4; 

Vivipara ussuriensis (Gerstfeldt), Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p108, pl. 
18, fig 1—4. 

Cipangopaludina_ ussuriensis (Gerstfeldt), Yen, 
1943, Nautilus, 56(4): 124-130; Liu et al, 1991, Proce- 
eding Tenth International Malacological Congress, 590. 

Material: 1 ex. IOZ-FG84819, collected from 
Khanka Lake, Heilongjiang Province in 1960. 

Diagnosis: shell larger, up to 54 mm, brown-gree- 
nish, with 6 whorls and two or three brown band. Apex 
obtuse. Spiral depth longer than aperture depth. Shell 
width is about 78.5% shell depth. Body whorl with 4 
clear ridges, other whorls with 2 ridges. Umbilious small. 

Distribution: Lower Amur River. 

Remarks: Cipangopaludina ussuriensis can be 
distinguished from other species of this genus by the 
following character: apex obtuse, body whorl with 4 
ridges and other whorl with 2 ridges. 


Cipangopaludina ventricosa (Heude, 1890) (Figure 1P) 
Paludina ventricosa Heude, 1890, Memoires conc- 
emant L’Histoire naturelle e L’empire chinois par des 
peres de la compagnie de Jesus, 175, pl. 39. 
Vivpara chinensis ventricosa (Heude): Kobelt, 1909, 
Abbildungen Nach de Natur mit Beschreibngen, p 114, pl. 
19, fig 1-2. 
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Cipangopaludina chinensis ventricosa (Heude), 
Pradshad, 1928, Memoirs of the Indian Museum, Calc- 
utta, 8: 168. 

Cipangopaludina ventricosa (Heude), Yen, 1943, 
Nautilus, 56(4): 124-130; Liu, 1991, Proceeding Tenth 
International Malacological Congress, 590. 

Materials: 24 ex. (1 female and 2 males were 
anatomy). KIZ 000031 collected from Yangguang Cou- 
nty, Yunnan in February 2006 by Du LN; KIZ 000036- 
38, KIZ 000174, KIZ 000227-228, KIZ 000230 collec- 
ted from Lake Dianchi, Yunnan in March 2006 by Du 
LN; KIZ 000143 collected from Lake Erhai, Yunnan in 
April 2006 by Du LN, Jing YE and Aldridge A; KIZ 
000175-180 collected from Jiuzai dragon spring, Yunnan 
in March 2006 by Du LN; KIZ 000226 collected from 
Lake Xingyun, Yunnan in April 2006 by Du LN; KIZ 
000283 collected from Black dragon spring, Songming 
County, Yunnan in December 2004 by Cui GH and Chen 
XY; KIZ 000361-362 collected from Black dragon 
spring, Songming County, Yunnan in December 2005 by 
Du LN; KIZ 000437-441 collected from Black dragon 
spring, Songming County, Yunnan in April 2006 by Du 
LN. 

Diagnosis: shell larger, medium thick, up to 65 mm, 
with 7 whorls. Body whorl inflated. Aperture inflated 
and round, but the upper angle does not forms a straight 
shoulder. Apex pointed. Umbilicus big and deep. 
Osphradium ridge-like, longer than length from anterior 
tip of osphradium to mantle border. Gill filaments same 
with C. chinensis. The dorsal length of kidney left plane 
about same length with ventral right plane, 1.5 times than 
posterior plane. Vas deferens with 3 main branches and 
the first anterior branch with 1 “V” shape and 3 parallel 
secondary branches. Vas deferens opens in last 1/4 of 
testis. Pouch filled with 33 embryos in various stages of 
development. Embryos with 3 whorls, body whorl with 3 
chaetae and other whorls with 2 chaetae. 

Distribution: Yunnan, Sichuan and Guizhou Provi- 
nces, China (Liu, 1991). 

Remarks: Cipangopaludina ventricosa is similar to 
C. cathayensis and C. ampullacea with a slightly larger 
shell, apex pointed and whorls inflated. However, C. 
ventricosa can be distinguished from C. cathayensis by 
following characteristics: shell with 7 whorls vs. 6; 
aperture round vs. ellipse; Osphradium ridge-like, longer 
than length from anterior tip of osphradium to mantle 
border vs, anterior 2/3 part of osphradium ridge-like and 
posterior 1/3 part inflated, osphradium length similar 
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with length from anterior tip of osphradium to mantle 
border; Vas deferens opens in last 1/4of testis vs. 1/3. 
Cipangopaludina ventricosa can be distinguished from C. 
ampullacea by the following characteristics: upper 
aperture does not form straight shoulder vs. having a 
straight shoulder; umbilicus big and deep vs. closed; 
Osphradium longer than length from anterior tip of 
osphradium to mantle border vs. shorter; Vas deferens 
opens in last 1/4 of testis vs. 1/3. 


Cipangopaludina yunnanensis Zhang, Liu et Wang, 
1981 (Figure 1M) 

Cipangopaludina menglaensis Zhang, Liu et Wang, 
1981, Acta Zootaxonomica Sinica, 6(1): 40-41. 

Materials: 48 ex. IOZ-FG 0152 (holotype), IOZ— 
FG164354—164400 (Paratypes) from Maliping village, 
Pu’er City, Yunnan in May 1957 by Zhang WZ and Liu 
YY. 

Diagnosis: shell greenish, small, up to 26 mm, with 
5 whorls. Apex obtuse. Shell surface with conspicuous 
growth lines and irregular malleation. Central tooth with 
wide rounded major denticles and 3 or 4 smaller trian- 
gular denticles on each side. Lateral tooth with tongue- 
shaped major denticles and 3 or 5 minor denticles on 
each side. Inner marginal tooth with tongue-shape major 
denticles and 3 or 4 minor denticles and marginal teeth 
with 13 almost equal-sized denticles. 

Distribution: only known from type locality. 

Remarks: Without anatomized specimens, the 
anatomy characters remain unknown. However, it can be 
distinguished from other species of Cipangopaludina by 
a smaller shell, up to 30 mm and marginal teeth with 13 
almost equal-sized denticles. 


Morphometric analysis by PCA 

The five shell characteristics are shown in Table 1. 
The first two factors PCA accounted for 85.4% of the 
total variance (Table 2). The first factor (55.9% of the 
variation of the selected variables) separated the ratios 
W/D, BW/D, BW/LA and BW/WA (factor loadings >80). 
The second factor (29.5% of the variation of the selected 
variables) separated the WA/D and WA/LA. A scatter 
map was made by averaging the first and second PCA 
factors per species (Figure 3). In PCA 1 axis (X-axis), W 
and BW are positively, D and LA negatively related to 
X-axis, respectively. In PCA 2 axis (Y-axis), WA is 
positive, D and LA negatively related to Y-axis, 
respectively. These variances are related to the shape of 
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aperture and body whorl depth. The pot of C. latissima 
located in the right of region II, it indicated that this 
species has a taller body whorl and shell width, but has a 
shorter spiral whorl. To the contrary, C. haasi looks 
slimmer due to taller spiral whorl and smaller shell width. 
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Figure 3 Scatterplots of scores on 1“ and 2” priniciple 
components of Chinese Cipangopaludina 


Table 2 Loadings of the first two principal components for 
5 shell characteristics of Chinese Cipangopaludina species 





PCA 1 PCA 2 
W/D 0.832* 0.322 
LA/D 0.571 0.513 
WA/D 0.397 0.914* 
BW/D 0.949* 0.288 
WA/LA 0.078 0.858* 
BW/LA 0.856* 0.001 
BW/WA 0.727* —0.682 
Prp. Total 55.9 29.5 


* loadings> 0.80. 


Additionally, C. haasi, located in left bottom of 
region II, has an oval aperture. Cipangopaludina 
hainanensis and C. cathayensis located in the top of 
region I and II, they have a bigger WA and relatively 
smaller LA, indicating these two species have a round 
aperture. 

Key to species of Chinese Cipangopaludina 
1 Shell with brown color band and ridge 


— Shell without brown color band or ridge eee 2 
2 Apex pointed a E A T 3 
= Apex obtuse OF eroded «os eresceeeeeteeeeeeeeceeeeeceeceneeeseeeeees 11 


3 The upper of body whorl form a flat shoulder 
eee eee eee eee eee eee eee eee eee eee eee eee eer ere eee eres C. lati ssi ma 
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The upper of body whorl without a flat shoulder ------- 4 
4 Shell with 5 whorls, spiral whorl not inflated 

PEIEE PAE TE E EA SEAE C. aubryana 
— Shell with 6-7 whorls, and spiral whorl inflated --.----- 5 
Vas deferens opens in last 1/3 of testis -+++ 6 


Nn 


— Vas deferens opens in last 1/4 of testis -+++ 10 


an 


The first anterior branch of vas deferens with more 
than 5 secondary branches ++-+++sssseseeseeeeeesneeeenee 7 
— The first anterior branch of vas deferens with less than 
3 secondary branches.---s-ssssscsssessssssssessnsessesssneeesness 9 
Shell width shorter than body whorl depth; the first 
anterior branch of vas deferens with 8 secondary 


N 


branches eenei enie o E EE EEEa C. haasi 
— Shell width longer than body whorl depth; the first ant- 
erior branch of vas deferens with 5-6 secondary bra- 


oo 


Shell thin; umbilicus small, nearly closed 
asna a Te a O aan n E e C. chinensis 
— Shell thick; umbilicus big and deeply -----C. c. fluminalis 
9 Vas deferens with 3 main branches-------- C. cathayensis 
— Vas deferens with 2 main branches .-.--..------ C. lecythoides 
10 Vas deferens with 2 main branches; kidney triple 
Shia PC ris riaren o A E OES C. lecythis 
— Vas deferens with 3 main branches; kidney arch shape 
ey ee ee ee ee ee ee C. ventricosa 


11 Shell medium, with 4—5 whorls---------- C. hainanensis 
— Shell larger, with 6 whorls -sesser 12 
12 Vas deferens with 2 main branches, opens in last 1/4 of 
A A EE EE C. l. ampullacea 
— Vas deferens with 3 main branches, opens in last 1/3 
OF LESEIS Hi cxies taane an C. leucostoma 
DISCUSSION 


Alongside numerous difficulties in determining the 
validity of numerous Cipangopaludina—lack of anatomy 
samples, discrepancies in the literature, or ambiguous 
morphological data—the genus itself is also problematic. 
Smith (2000) treated Cipangopaludina as subgenus of 
Bellamya, while Sengupta et al (2009) indicated that 
species of Bellamya from Asia could not belong to this 
genus. In China, there are currently seven recognized 
genera of Bellamyinae recorded, viz. Margarya Nevill, 
1877, 1912, Bellamya 
Jousseaume, 1886, Mekongia Crosse et Fischer, 1876, 


Cipangopaludina Hannibal, 


Filopaludina Habe, 1964, Angulyagra Benson, 1836 and 
Trochotaia Brandt, 1974 (Liu et al, 1979, 1993; Du et al, 
2011). Generally speaking, it is possible to distinguish 
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these, as the shell of Cipangopaludina, Margarya, and 
Trochotaia tend to be larger, while other four genera 
species are more mid-sized. Similarly, species of 
Filopaludina and Mekongia have a color band, and 
Angulyagra has ridges. For more detailed markers, 
Cipangopaludina can be distinguished from Trochotaia 
by shell round or oval vs. depressed pyramidally, last 
whorl with a sharp keel, and from Margarya by shell 
smooth vs. shell with ridges, nodules or spines; embryo 
with chaetae vs. ridges or nodules. Based on these 
criteria, we opted to treat Cipangopaludina as a valid 
genus. 

Generally, the diagnostic characteristics of 
Cipangopaludina are shell median to large, sculpture 
lacking except growth lines and axial undulations. 
Aperture subcircular. Young specimens (from brood 
pouch) with up to 3—4 whorls; protoconch smooth, three 
primary rows of chaetae on last whorl and other whorls 
have two rows of chaetae. However, C. menglaensis and 
C. yunnanensis are small, up to 30 mm, smooth and solid, 
suggesting that C. menglaensis and C. yunnanensis 
belong to the genus Mekongia. It is also possible that C. 
dianchiensis could belong to the genus Margarya due to 
several characteristics: shell thick with ridges, female 
pouch with few big embryos and embryos shell with 
ridges (not chaetae). Similarly, Cipangopaludina 
ampullacea and C. lecythis, C. fluminalis and C. 
chinensis can be distinguished by shell characteristics, 
but they do not appear to possess any obviously 
distinctive anatomy, placing C. ampullacea and C. 
fluminalis as a subspecies of C. lecythis and C. chinensis, 
respectively. Our analysis also placed C. leucostoma 
(Heude) and C. diminuta (Heude) as synonymous with C. 
patris (Kobelt) according to apex obtuse and less than 6 
whorls. However, according to the International Code of 
Zoological Nomenclature (ICZN, 1999), C. leucostoma 
should be valid due to its publication earlier than C. 
and C. 


ampulliformis were found to be synonyms of C. 


patris. Finally, Cipangopaludina wingatei 
chinensis and C. lecythis, respectively. 

Totally, our analysis of 18 Cipangopaludina species 
resulted in a revised taxonomy that includes 11 species 
and 2 subspecies recorded in China, viz. C. aubryana, C. 
cathayensis, C. chinensis, C. haasi, C. hainanensis, C. 
latissima, C. lecythis, C. lecythoides, C. leucostoma, C. 
ussuriensis, C. ventricosa and C. chinensis fluminalis, C. 
lecythis ampullacea. Based on whether or not the spiral 
whorl depth was longer than aperture depth, Chinese 
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Cipangopaludina species could be divided to two groups, 
the chinensis group with spiral whorl depth longer than 
aperture depth that includes C. chinensis, C. haasi, C. 
ussuriensis and C. c. fluminalis, and the other the 
cathayensis group with spiral whorl depth being shorter 
than aperture depth, including C. aubryana, C. catha- 
yensis, C. hainanensis, C. latissima, C. lecythis, C. lecyt- 
hoides, C. leucostoma, C. ventricosa and C. l. ampul- 
lacea. Several anatomical characteristics can also help 
distinguish between the two groups. For example, the 
chinensis group can be distinguished from cathayensis 
group by vas deferens: chinensis group species with 
more than 5 (C. haasi with 8) parallel secondary bran- 
ches in anterior branch, while cathayensis group species 
always with 4 secondary branches. The number of 
secondary branches of vas deferens are in connection 
with size of body whorl. The shell of chinensis group 
species was taller, with body whorl depth being the same 
or longer than body width, making the vas deferens 
longer and with more secondary branches, especially in 
C. haasi. However, the body whorl of cathayensis group 
species is inflated, with body whorl depth shorter than 
body width, resulting in vas deferens being short and fat, 
and the secondary branch with many tiny branches, 
looking rather similar to a lush tree. 

Within each group, there are further defining 
characteristics of individual species. In chinensis group, 
C. ussuriensis can be easily distinguished from other 
species by body whorl with 4 ridges and colour band. 
Cipangopaludina haasi can be distinguished from other 
species by body whorl not inflated, vas deferens with 8 
secondary branches in anterior branch and vas deferens 
opens in last 1/4 of testis. Cipangopaludina chinensis 
fluminalis, however, is quite similar with C. chinensis, 
except for the shell thick vs. thin; umbilicus big and 
deeply vs. small or closed. In cathayensis group, C. 
latissima could be distinguished from other species by its 
markedly enlarged aperture while C. hainanensis can be 
distinguished from other species by apex eroded, left 4 
whorls. According to kidney shape, the other species 
could be divided to two subgroups. Triple shape, 
including C. cathayensis, C. aubryana and C. lecythis, 
and arch shape, including C. l ampullacea, C. 
leucostoma and C. ventricosa. Cipangopaludina lecythis 
can be distinguished from C. cathayensis and C. 
aubryana by vas deferens with 2 branches and opens in 
last 1/4 of testis vs. 3 branches and 1/3, and C. catha- 
yensis can be distinguished from C. aubryana by the 
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anterior 2/3 part of osphradium ridge-like and posterior 
1/3 part inflated. Cipangopaludina leucostoma can be 
distinguished from C. l. ampullacea and C. ventricosa by 
vas deferens opening in last 1/3 vs. 1/4, and C. ventricosa 
can be distinguished from C. l. ampullacea by vas 
deferens with 3 branches vs. 2 branches. 

Similar to our study, most viviparids are primarily 
delineated by their shells, though a few observations on 
the anatomy have also been made. Vail (1977) for 
example mentioned that the reproductive system was a 
useful characteristics to delineate viviparid subfamilies 
and genera by anatomical studies of mature male and 
female Campeloma geniculum (Conrad, 1834), Lioplax 
pilsbryi Walker, 1905 and Viviparus georgianus (Lea, 
1834). Simone (2004) and Rao (1925) also reported 
intergeneric differences in the central nervous system, 
stomach lining and form of the ctenidial filaments in the 
five Asian genera. Du et al (2011) later reported a 
difference in kidney shape and the number of embryos 
present in the uterus of four Chinese viviparids, 
Angulygra, Cipangopaludina, Margarya and Trochotaia. 
In this study, the size of osphradium, location of vas 
deferens opens and the number of secondary branches, 
kidney shape were useful anatomical characteristics to 
identify interspecies of Cipangopaludina. The allometric 
environmental effects could result in 
shell 
operculum (Chiu et al, 2002), but these anatomical 


growth or 
intraspecific variations in morphology and 
characteristics could be used to identify young specimen. 

Given the difficulties in identifying species solely 
by shell characters, further description of such species 
anatomical characteristics of both young and mature 
individuals may prove useful in further refining the 
taxonomy and identifying species and determining their 
further 
more accurately characterize the 


validity. Unfortunately, until studies are 
undertaken that 
anatomy or even underlying genetic differences, there 
will remain a variety of unanswered or unanswerable 
questions. To date, the monophyly of Cipangopaludina is 
still unresolved. A recent study by Du et al (2013) on the 
phylogeny of Margarya based on combined CO I and 
16S rRNA sequences indicated that Cipangopaludina is 
not a monophyly, since species cluster with Margarya; 
however, the specimens of Cipangopaludina and 
Margarya are from Lake Dianchi. Likely further studies 
into the phylogeny and ecology of the Cipangopaludina 
in different environments and distributions will further 


enhance our understanding of these organisms. 
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Abstract: Cervus sichuanicus is a species of sika deer (Cervus nippon Group). To date, research has mainly focused on quantity 
surveying and behavior studies, with genetic information on this species currently deficient. To provide scientific evidence to assist in 
the protection of this species, we collected Sichuan sika deer fecal samples from the Sichuan Tiebu Nature Reserve (TNR) and 
extracted DNA from those samples. Microsatellite loci of bovine were used for PCR amplification. After GeneScan, the genotype 
data were used to analyze the genetic diversity and population structure of the Sichuan sika deer in TNR. Results showed that the 
average expected heterozygosity of the Sichuan sika deer population in TNR was 0.562, equivalent to the average expected 
heterozygosity of endangered animals, such as Procapra przewalskii. Furthermore, 8 of 9 microsatellite loci significantly deviated 
from the Hardy-Weinberg equilibrium and two groups existed within the Sichuan sika deer TNR population. This genetic structure 


may be caused by a group of Manchurian sika deer (Cervus hortulorum) released in TNR. 


Keywords: Sichuan sika deer; Microsatellite; Genetic diversity; Population structure 


The Sichuan sika deer (Cervus sichuanicus) (Guo et 
al, 1978) is a species of the Cervus nippon group (Groves 
and Grubb, 2011). It is an endemic species in China and 
is listed as a class ] endangered species, as well as an 
endangered species by the International Union for 
Conservation of Nature and Natural Resources (IUCN). 
Currently, the Sichuan sika deer is distributed in three 
unconnected areas located in northwest Sichuan Province 
and southwest Gansu Province. Compared to other spec- 
ies in China, Sichuan has the largest sika deer population 
in the wild, consisting of about 850 individuals (Guo, 
2000). 

In recent years, most research has focused on 
macroscopic aspects such as social behavior (Guo et al, 
1991), distribution and quantity surveying (Guo, 2000), 
activity rhythms (Liu et al, 2004), and food habits (Guo, 
2001). At the genetic level, Wu et al (2008) chose 16 
microsatellite loci to analyze the genetic diversity of Sika 
deer in China, including 24 samples of Sichuan sika deer. 
Lü et al (2006) and Liu et al (2003) studied the genetic 


Science Press 


relationship of sika deer using mitochondrial DNA, 
which also had a small number (16) of Sichuan sika deer 
samples. These studies focused on the differences among 
species and genetic structure of sika deer in China. Due 
to the limited samples, however, the above data do not 
reflect genetic information in one species or population. 
Moreover, using traditional capture methods makes it 
difficult to obtain abundant samples to study genetic 
diversity and population structure within Sichuan sika 
deer populations. 

The Tiebu Nature Reserve (TNR) in the Zoige 
administrative district of northwest Sichuan has the 
largest sika deer population. According to monitoring 
data from TNR, there were about 1 000 Sichuan sika deer 
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in the reserve in 2012. TNR has been chosen as the main 
area to study Sichuan sika deer due to large population 
size and better protection situation (Li, 2010; Ning et al, 
2008; Qi et al, 2010). To date, however, no genetic 
information studies on TNR sika deer have been 
conducted and the genetic diversity and possible 
population substructures remain unknown. Previous 
studies have shown that the geographical distribution 
pattern of the population is related to geographical 
isolation (Boyce et al, 1999; Gavin et al, 1999; Scandura 
et al, 1998). However, it is unclear whether the rivers, 
roads, villages, and crop protection fences across the 
reserve affect gene flow or generate genetic structures 
within the sika deer population. Therefore, we 
hypothesized that geographic and artificial factors may 
affect the gene flow of sika deer, and thus established 
seven sampling areas in TNR to obtain sika deer fecal 
samples. Overall, we collected five hides, one tissue 
sample, and 149 fecal samples. After extracting DNA 
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from those samples, nine Bovinae microsatellite loci 
were amplified and the products were used for gene 
scanning. We then used the genotype data to analyze 
genetic diversity and population structure. The results of 
this study provide a scientific basis for the protection of 
the Sichuan sika deer. 


MATERIALS AND METHODS 


Study area and sampling collection 

Fecal samples were collected in seven areas in TNR 
(E102°46'-103°14', N33°58’-34° 16’, about 206 922 km”), 
as shown in Figure 1. To avoid cross-contamination 
when sampling, we only collected one pellet from each 
fecal heap. We collected five hides, one tissue sample, 
and 149 fecal samples. The hide samples were collected 
from dead individuals in the wild. The only tissue sample 
came from an individual that had been killed by poachers. 
All samples were stored in 95% ethanol. 


Tiebu Nature Reserve 


Figure 1 Distribution of sampling areas in Tiebu Nature Reserve 


DNA extraction and amplification 

DNA was extracted using an E.Z.N.A HP Tissue 
DNA Maxi Kit and a Stool DNA kit (OMEGA) 
according to the manufacturer’s protocols. The former 
primers of nine Bovinae microsatellite loci (IDVGA29, 
BM4107, RT1, TGLA53, RM188, RT13, BM6506, 
BL42, CSSM019, ETH225, and CERVID14) were mo- 
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dified fluorescent groups. The PCR mixture contained 
1.5 uL genomic DNA (100 ng/uL), 1 uL of each pri- 
mer (10 pmol/uL), 0.5 uL BSA (20 mg/uL), and 7.5 uL 
of Premix Ex Taq (Takara), with deionized water used to 
make the sample volume up to 15 uL. The PCR amp- 
lification was carried out using an initial denaturation at 
94 °C for 3 min, followed by 35 cycles at 94 °C for 30 s, 
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annealing for 30 s, primer extension at 72 °C for 45 s, 
and a final extension at 72 °C for 10 min. The 
products were preserved at 4 °C. One tissue sample 
was used as positive control and water as negative 


Table 1 Information on the nine microsatellite bovine loci used in this research 


Annealing 


control to ensure that PCR was performed correctly 
and to avoid contamination. The PCR products were 
sent to Shanghai MAP Biotechnology Co., Ltd 
(Shanghai, China) for GeneScan. 


Sources of 





Locus Repeat motif Primer sequence Size range (bp) temperature (C) primers 

CCC ACA AGG TTA TCT ATC TCC AG 

IDVGA29 (AC)n 137-157 60 Konfortov et al 
CCA AGA AGG TCC AAA GCA TCC AC (1996) 
AGC CCC TGC TAT TGT GTG AG ‘ 

BM4107 (AC)n (TC)n (TG)n 162-172 60 Bishop et al 
ATA GGC TTT GCA TTG TTC AGG (1994) 
TGC CTT CTT TCA TCC AAC AA i ; 

RT1 (GT)n 222—240 55 Mieg it 
CAT CTT CCC ATC CTC TTT AC (1998) 
CAG CAG ACA GCT GCA AGA GTT AGC 

TGLA53 (AC)n 176-216 50 Hoffmann et al 
CTT TCA GAA ATA GTT TGC ATT CAT GCA G (2004) 
GGG TTC ACA AAG AGC TGG AC 

RM188 (AC)n 141-151 50-52 Barendse et al 
GCA CTA TTG GGC TGG TGA TT (1994) 
GCC CAG TGT TAG GAA AGA AG . 

RT13 (GT)n 293-324 60 ree al 
CAT CCC AGA ACA GGA GTG AG (1994) 
GCA CGT GGT AAA GAG ATG GC 

BM6506 (AC)n 196-212 59-63 Moore et al 
AGC AAC TTG AGC ATG GCA C (1994) 
CAA GGT CAA GTC CAA ATG CC i 

BL42 (AC)n 246-258 58-60 Bishop et al 
GCA TTT TTG TGT TAA TTT CAT GC (1994) 
TTG TCA GCA ACT TCT TGT ATC TTT ; ; 

CSSMO019 (AC)n 138-156 52-58 Wilson & White 
TGT TTT AAG CCA CCC AAT TAT TTG (1998) 
GAT CAC CTT GCC ACT ATT TCC T 

ETH225 (AC)n 140-193 60 Steffen et al 
ACA TGA CAG CCA GCT GCT ACT (1993) 
TCT CTT GCG TCT CCT GCA TTG AC 

CERVID14 (AC)n 214-231 61-65 De Woody stal 
AAT GGC ACC CAC TCC AGT ATT CTT C (1995) 


Species identification and population assessment 

We amplified all samples using two pairs of primers 
of the mitochondrial control region and compared the 
sequences in GenBank to ensure the samples in our 
research were from sika deer. We used GeneMapper v. 4 
(Applied Biosystems) to obtain the microsatellite 
fragment sizes and identified the heterozygotes and 
homozygotes. We did two repeats if the genotypes of the 
two repetitive PCRs were different; for example, if there 
was a heterozygote in one repeat but a homozygote in the 
other repeat. If there was a heterozygote in three of four 
repeats, we considered this site to be a heterozygous 
locus; the same for a homozygote locus. The probability 
of identity (PI and the PI between siblings were 
calculated in Gimlet 1.3.3 (Valière, 2002; Liu & Yao, 
2013) using genotype data of the nine microsatellite loci. 
We used the following principle when using Mstools 
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(Park, 2001) to assess the population number: we 
allowed one genetic mismatch at one allele for one locus 
to be considered as two samples belonging to the same 
individual (Bellemain et al, 2005). To detect the validity 
of the above principle, we randomly selected 50 fecal 
samples (accounting for a third of the total number) and 
amplified the nine loci twice. At the same time, we 
amplified another two loci (ETH225 and CERVID14) 
of the selected 50 samples. We tested whether the 
influence of population number may be caused by the 
repeat amplification and increasing loci. 


Genetic diversity and population structure 

We computed the expected heterozygosity (He), 
observed heterozygosity (Ho), and polymorphism infor- 
mation content (P/C) using Mstools. Arlequin 3.1 (Exc- 
offier et al, 2005) and Fstat 2.9.3.2 (Goudet, 2001) were 
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used to test the Hardy-Weinberg equilibrium and Fis 
separately. The fixation index (Fsr) (Weir & Cockerham, 
1984) among sampling areas was tested in Arlequin 3.1, 
using 1000 random iterations to calculate confidence 
intervals (CI). Software structure 2.3 (Pritchard et al, 
2000) was used to estimate the potential populations (K). 
The admixture model was used and run three times using 
10 replications of K ranging from 1 to 7. For those runs 
the following conditions were adopted: a burn-in period 
at 10 000 following 10 000 000 replicates of the MCMC. 
The AK calculated as per Evanno et al (2005) was used 
to estimate optimal K. The formulas were: 
AK=m|L(K+1)—2L(K)-L(K-1)|/s[L(K)]._ We referred to 
the structure output InP(D) as L(K). m was the mean 
value and s was the standard deviation. Software 
CLUMPP Windows.1.1.2 (Jakobsson & Rosenberg, 2007) 
and destruct (Rosenberg, 2004) were used to visualize 








the individual coefficients of membership in the 


subpopulation calculated by structure 2.3. 


RESULTS 


Species identification and population assessment 

Comparing the obtained sequences of mitochondrial 
hypervariable region I from our research to the 
corresponding sequence in GenBank, all samples were 
obtained from sika deer. Thereinto, 143 of 155 samples 
belonged to the haplotypes of Sichuan sika deer and the 
remainder belonged to the haplotypes of Manchurian 
sika deer (Cervus hortulorum). 

Most of the Sichuan sika deer samples could be 
amplified at all the nine microsatellite loci, few samples 
(11) had a loss at one locus. Analyzing the united PI 
values for nine loci, we found that the discrimination 
power of nine loci was close to 100% (Figure 2). 
According to the rule of Bellemain, we assessed 76 sika 


deer individuals from 143 face samples. After the second 
amplification of the randomly selected samples, the 
population number did not change. The additional two 
loci (ETH225 and CERVID14) increased the population 
by one individual, accounting for 3% of the total number. 
The distribution of the 76 individuals in the samplings 
areas is shown in Table 3. 
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Figure 2 Multi-loci PI (biased and unbiased) and PI (sibs) for 
sika deer in increasing order of single-locus values 

Genetic diversity and population structure 

In the Sichuan sika deer population in TNR, the 
mean number of alleles / locus was 6.56 and the Hg was 
0.562. The PIC values were higher than or close to 0.5, 
except for RT13 (0.129). All detected loci deviated from 
the Hardy-Weinberg equilibrium significantly (Figure 2) 
and the P value was 0. In addition to the locus RT13, the 
Fis values of the remaining eight loci were negative. The 
values of Hg ranged from 0.4292 to 0.6865 in the seven 
sampling areas and all Ho values were higher than those of 
Hy (Table 3). The Fsr among the seven sampling areas 
ranged from —0.04781 to 0.05243. The highest Fr 
(0.05243) appeared between Area IV and VI but the lowest 
(—0.04781) appeared between Area I and III (Table 4). 


Table 2 Microsatellite diversity indices of Sichuan sika deer in Tiebu Nature Reserve 





Locus No. of alleles Hg Ho PIC Fis 
IDVGA29 5 0.587 0.722 0.516 —0.230 
BM4107 7 0.560 0.654 0.495 —0.169 
RT1 6 0.633 0.951 0.558 —0.506 
CSSM019 6 0.635 0.926 0.560 —0.462 
TGLA53 8 0.643 0.899 0.571 —0.401 
BL42 9 0.607 0.913 0.523 —0.509 
RM188 6 0.656 0.738 0.607 —0.125 
BM6506 8 0.604 0.974 0.519 —0.620 
RT13 4 0.134 0.026 0.129 0.810 
Mean 6.56 0.562 0.756 0.496 —0.347 
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Table 3 Indices of genetic diversity of Sichuan sika deer in different sampling areas 





Sample area Sample size Loci typed Ag Ho 
I 10 9 0.4292 0.7160 
I 13 9 0.5499 0.8034 
MI 7 9 0.4831 0.7778 
IV 8 9 0.6804 0.7284 
V 20 9 0.6865 0.7444 
VI 11 9 0.4350 0.7063 
VI T 9 0.5368 0.8571 


Table 4 F-statistics (Fsr, below diagonal) and significant differences among different areas 











Sample area I I M IV V VI VII 
+ = + + + + 
II —0.00493 — — — + — 
M —0.04781 —0.03753 — — — = 
IV 0.03553 —0.0052 —0.00225 = F = 
V 0.05177 0.00789 0.01646 —0.01778 F = 
VI —0.02823 —0.00142 —0.0334 0.03287 0.05243 + 
VII 0.01389 —0.02825 —0.03837 —0.01869 0.01487 —0.00278 
+ denotes significant difference, — in contrast (95% confidence interval). 
Population structure 
Structure analysis showed that when K=2, the 4K DISCUSSION 


exhibited an obvious apex (Figure 3). The Fs; between 
subpopulations was 0.535. Individual coefficients of 
membership in subpopulations (Figure 4) showed that 
when K=2, the genetic amount of cluster 1 (Yellow) was 
larger and distributed in all sampling areas. In contrast, the 
genetic amount of cluster 2 (Blue) was less than cluster 1, 
and had a relatively smaller distribution area. Most of the 
blue genotype was distributed in Areas IV and V. 
80 
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Figure 3 Value of 4K as a function of K based on 10 runs 
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Figure 4 Individual coefficients of membership in each 
sampling area 
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Population assessment 

In recent years, non-invasive sampling methods 
have been widely used in genetic studies of endangered 
animals (He et al, 2010; Wang et al, 2012). Because of 
the endangered status and fugacious feature of sika deer, 
non-invasive sampling methods (obtaining fecal samples) 
have become an important means of population 
assessment and individual identification. In this study, 
the stability of repeat amplification of nine microsatellite 
loci had no effect on assessing the population number, 
that is, our genotype data were reliable. After adding two 
loci, there was a 3% error between previous and post 
assessment results, which is much higher than 2% 
(Bellemain et al, 2005). Bellemain et al (2005) attributed 
such errors in population assessment to pseudogenes, 


pollution and missing data. 


Genetic diversity and differentiation 
Genetic diversity helps populations adapt to 
changing environments. With more variation, it is more 
likely that some individuals in a population will possess 
variations of alleles that are suited for the environment. 
The population will continue for more generations 
because of the success of these individuals (Khater et al, 


2011). In our research, the Hg of Sichuan sika deer 
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population was 0.562, close to the mean Hg (0.559) of 
sika deer in China, but much higher than the mean Hg for 
Sichuan sika deer (0.477) and slightly lower than the 
mean Hg for Manchurian sika deer (0.584) and south 
China sika deer (Jiangxi, 0.585; Zhejiang, 0.589) (Wu et 
al, 2008). Compared to other Ungulata, such as red deer 
(Cervus elaphus) (H,=0.78, Hajji et al, 2007; Hp=0.804, 
Pérez-Espona et al, 2008; Hp=0.62—0.85, Zachos et al, 
2007), forest musk deer (Moschus berezovskii) 
(Hp=0.8-0.9, Zhou et al, 2005), and Mongolian wild ass 
(Equus hemionus) (Hy=0.83, Kaczensky et al, 2011), 
Sichuan sika deer had a relatively lower Hg and was 
equivalent to Przewalski's Gazelle (Procapra przewalskii) 
(H:=0.552, Yang & Jiang, 2011) and milu (Elphurus 
davidianus) (Hg=0.46—0.54, Zeng et al, 2007). The Hgs 
values were different in each sampling area. Area | had 
the lowest Hg (0.4292), while Area V had the highest 
(0.6865). Polymorphism information content (PIC) is a 
measure of variation between microsatellite loci. It is 
generally considered that a PIC value greater than 0.5 is 
for highly polymorphic loci, a PIC value of less than 0.5 
and greater than 0.25 is for moderately polymorphic loci, 
and a PIC of less than 0.25 is for low polymorphic loci 
(Niu et al, 2001). In this study, the PIC values of six of 
the nine microsatellite loci were higher than 0.5, 
indicating highly polymorphic loci, while two were 
moderately polymorphic loci and one was a low 
polymorphic locus. Heterozygosity and PIC information 
indicated Sichuan sika deer in TNR had abundant genetic 
diversity. 

Hartl & Clark (1997) considered that population 
deviation from the Hardy Weinberg equilibrium was 
mainly due to small populations, not random mating, 
gene mutation, migration or other factors. The present 
study showed that nine microsatellite loci of the Sichuan 
sika deer population in TNR deviated significantly from 
the Hardy Weinberg equilibrium. The Sichuan sika deer 
population numbers in TNR grew steadily from 650 in 
2000 (Guo, 2000) to 1000 in 2012 (monitoring data from 
TNR). As a result, the deviation from the Hardy 
Weinberg equilibrium is unlikely attributed to a 
reduction in the population. The coefficient of genetic 
differentiation (Fsr), on the other hand, is an important 
indicator of genetic differentiation among subpopul- 
ations. According to Wright (1978), Fsrs values bet- 
ween 0-0.5 suggest no differentiation between sub- 
groups. The Fsrs values in our seven sampling areas were 
all less than 0.5, ranging from 0.04781 to 0.05243, and 
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showing negative value in some regions. This indicated 
that the genetic variation among the seven sampling 
areas mainly came from inside the sampling area, and the 
geographical barriers and human factors do not affect 
gene flow in Sichuan sika deer in TNR. In addition, due 
to the relatively closed geographical position of TNR, 
very few little sika deer were distributed in the periphery 
of the protected areas (Guo, 2000). Gene flow through 
immigration and emigration could not significantly affect 
the Sichuan sika deer populations to deviate from the 
Hardy Weinberg equilibrium. 

Microsatellite polymorphism studies of sika deer in 
China found that 16 microsatellite loci significantly 
deviated from the Hardy Weinberg equilibrium, but Hgs 
were higher than Hos (Wu et al, 2008). In the present 
study, Hos were significantly higher than the Hgs. The 
Fis values of eight microsatellite loci were negative, 
showing heterozygosis excess. When a population 
experiences a reduction in its effective size, it generally 
develops heterozygosity excess at selectively neutral loci 
(Cornuet & Luikart, 1996). Xu et al (2009) used 
microsatellites to detect the population structure of 
Himalayan marmots and found heterozygosity excess 
due to historic population decline. Wu et al (2008) tested 
bottlenecks in four Chinese sika deer species and found 
that Sichuan sika deer had not experienced a population 
Therefore, the 
heterozygosity excess of Sichuan sika deer in TNR could 


bottleneck in recent history. 


not be caused by population reduction. 


Population structure 

A population composed of individuals from two 
different source populations will tend to have an excess 
of heterozygotes (Milkman, 1975). In testing population 
structure in the present study, we found that when K=2, 
AK exhibited a significant apex, indicating that two 
subpopulations existed within the Sichuan sika deer 
population. Considering that the Sichuan sika deer 
population had heterozygote excess and two subpo- 
pulations, we speculated that the sika deer in TNR might 
have two different sources populations. According to Wu 
(2004) and TNR staff, a captive-bred Manchurian sika 
deer group was bred and released in TNR in the late 
1960s. The mitochondrial data showed that haplotypes of 
Manchurian sika deer (unpublished data) do exist. 
Previous research has shown that introduced species 
often hybridize with local sibling species (Avise et al, 
1974; Gutiérrez, 1993; Perry & Goudet, 2001) and 
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panmictic admixtures can develop rapidly over a short 
period of time (Echelle & Connor, 1989; Echelle, 1987; 
1997). 
introduced sika deer have been shown to hybridize with 
native red deer (Senn & Pemberton, 2009). Therefore, 
there is a strong possibility that hybridization occurred 


Hybridization abounds in (Cervidae, and 


between the two sika deer species in our research. The 
released Manchurian sika deer explained the existence of 
subpopulations in the Sichuan sika deer population in 
TNR. We speculate that the Manchurian sika deer 
genotype is the minority — the blue parts in the individual 
coefficients of membership, which were mainly 
distributed in Areas IV and V- because the number of 
introduced Manchurian sika deer is not clear and of the 
155 sika deer fecal samples examined, only 12 belonged 
to the haplotypes of Manchurian sika deer. The Sichuan 
sika deer population were in larger quantities than that of 


Manchurian sika deer. So, the yellow parts in Figure 4 
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Abstract: The yellow meal worm (Tenebrio molitor L.) is an important resource insect typically used as animal feed additive. It is 
also widely used for biological research. The first complete mitochondrial genome of T. molitor was determined for the first time by 
long PCR and conserved primer walking approaches. The results showed that the entire mitogenome of T. molitor was 15785 bp long, 
with 72.35% A+T content [deposited in GenBank with accession number KF418153]. The gene order and orientation were the 
same as the most common type suggested as ancestral for insects. Two protein-coding genes used atypical start codons (CTA in ND2 
and AAT in COXI), and the remaining 11 protein-coding genes started with a typical insect initiation codon ATN. All tRNAs showed 
standard clover-leaf structure, except for tRNA*” (AGN), which lacked a dihydrouridine (DHU) arm. The newly added T. molitor 


mitogenome could provide information for future studies on yellow meal worm. 


Keywords: Tenebrionidae; Mitogenome; Evolution; Resource insect; Yellow meal worm 


Mitochondrial genes have been widely used in 
studies of phylogenetics, phylogeography and molecular 
diagnostics (Kocher et al, 1989; Simmons & Weller, 
2001; Wolstenholme, 1992). The full-length mitochon- 
drial genome (mitogenome) has become an important 
tool for studies of genome architecture, population gene- 
tics, primer design, and molecular evolution (Kim et al, 
2009; Lee et al, 2009; Liao et al, 2010; Yukuhiro et al, 
2002). Furthermore, mitogenome sequences provide im- 
portant information for pest control. For example, insec- 
ticide resistance in arthropod pest Tetranychus urticae is 
controlled by mtDNA (Van Leeuwen et al, 2008). 

Tenebrionidae is a large family of Coleoptera and 
more than 20000 species have been described worldwide 
(Bouchard et al, 2011). Despite this large taxonomic 
diversity, information about the Tenebrionidae mito- 
genome is still limited. To date, only two complete mit- 
ogenomes are available in GenBank (NC_003081: Trib- 
olium castaneum and NC_013554: Adelium sp. NCS- 
2009). Due to the importance of mitogenomes in popul- 
ation studies and pest control, studies on the mitogenome 
of tenebrionids are required. 

The yellow meal worm (Tenebrio molitor L.) is an 
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important resource insect, which has great value as an 
animal feed additive because of its high protein content. 
It has also been used in numerous experimental studies 
(Gomez et al, 2013; Simon et al, 2013). Its relatively 
large body-size and ease of rearing make it attractive for 
biological research. An increasing number of researchers 
have used this beetle as a model system for studies in 
biology, biochemistry, evolution, immunology and 
physiology (Dobson et al, 2012; Goptar et al, 2013; Ohde 
et al, 2013; Prokkola et al, 2013; Yan et al, 2013; Zhu et 
al, 2013). However, the genetic background of this beetle, 
in both nuclei and organelles, has not yet been determ- 
ined, and information is urgently needed. 

In this study, we sequenced the whole mitogenome 
of T. molitor for the first time and analyzed its nucleotide 
composition and major characteristics. The entire 
mitogenome was found to be 15 785 bp long and harbo- 
red 13 protein-coding genes (PCGs), 22 tRNA genes, and 
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2 rRNA genes. The gene arrangement in the T. molitor 
mitogenome was identical to the ancestral insect mitog- 
enome arrangement (Boore et al, 1998). The mitogenome 
of T. molitor is an important addition to the continuing 
efforts in developing T. molitor as a model system for 
biological studies. 


MATERIALS AND METHODS 


Sample and DNA extraction 

An adult T. molitor was collected in Kunming, 
Yunnan Province, China, on July 29, 2012. The adult 
insect was one of the descendants from our inbred line, 
which has been established and maintained in our 
laboratory since May 2008. All members of this inbred 
line have one common ancestral mitochondrial geno- 
me background. After collection, the fresh materials 
were immediately preserved in 100% ethanol and 
stored in a —20 °C refrigerator before genomic DNA 
extraction. Total genomic DNA was extracted with the 
WizardTM Genomic DNA Purification Kit, in accor- 
dance with the manufacturer's instructions (Promega, 
USA). 


PCR amplification and sequencing 

To sequence the complete mitogenome, long PCR 
primers and some short PCR primers were designed by 
multiple sequence alignments of all available complete 
Coleoptera mitogenomes using ClustalX1.8 (Thompson 
et al, 1997) and Primer Premier 5.0 software. Primer 
sequence information can be obtained from the authors 
upon request. Long PCRs were performed using TaK- 
aRa LA Taq polymerase with the following cycling 
parameters: initial denaturation for 5 min at 95 °C, 
followed by 30 cycles at 95 °C for 50 s, 50 °C for 50 s, 
68 °C for 2 min and 30 s; and a final extension step of 
68 °C for 10 min. Short fragments were amplified 
with TaKaRa Taq polymerase: initial denaturation for 
5 min at 94°C, followed by 35 cycles at 94 °C for 1 
min, 45—53 °C for 1 min, 72 °C for 2 min, and a final 
extension step of 72 °C for 10 min. The PCR 
products were detected via electrophoresis in 1.5% 
agarose gel, and purified using the QIAquick PCR 
(QIAGEN, USA). 
reaction was performed using a BigDye Terminator 


Purification Kit Sequencing 
v3.1 Cycle Sequencing Kit (Applied Biosystems, 
City, CA). Electrophoresis 
sequencing products was performed on an ABI-3730 
DNA Analyzer (Applied Biosystems). 


Foster of purified 
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Gene identification and tRNA structures 

Sequences from overlapping fragments were assem- 
bled with the neighboring fragments using the SeqMan 
program included in the Lasergene software package 
(DNAStar Inc., Madison, Wisc.). Comparison of the 
DNA or amino acid sequences with homologous regions 
of known full-length insect mitogenome sequences 
indicated 13 protein-coding genes (PCGs) and two rRNA 
genes based on MEGA 5.0 software (Tamura et al, 2011). 
The nucleotide sequences of the PCGs were translated on 
the basis of the invertebrate mtDNA genetic code. 
Transfer RNA gene analysis was conducted based on 
tRNAscan-SE software v.1.21 (Lowe & Eddy, 1997), and 
the folding of the predicted tRNA sequences were further 
confirmed by visual inspection. Finally, the T. molitor 
mitogenome sequence was deposited in GenBank with 
accession number KF418153. 


Sequence analysis 

The A+T-content of the whole genome was calcul- 
ated based on the EditSeq program included in the 
Lasergene software package. The nucleotide composition 
at each codon position of the PCGs and codon usage 
were calculated based on MEGA 5.0. Gene overlaps and 
intergenic-space sequences were hand-counted. 

Nucleotide composition skew was calculated for the 
PCGs and the whole genome base on the EditSeq prog- 
ram included in the Lasergene software package (DN- 
AStar Inc., Madison, Wisc.), using the following formula 
proposed by Perna & Kocher (1995): GC-skew=(G— 
C)/(G+C) and AT-skew=(A-T)/(A+T), where C, G, A and 
T are the frequencies of the four bases. 


RESULTS AND DISCUSSION 


General features of Tenebrio molitor mitogenome 

The complete mitogenome of T. molitor was 15 785 
bp in length, and consisted of 2 rRNAs (srRNA and 
IrRNA), 13 PCGs (ATP6, ATP8, COI-III, ND1-6, ND4L, 
CytB), and 22 tRNAs (Figure 1). Like many other insect 
mitogenomes, the major strand (J strand) carried most of 
the genes (9 PCGs and 14 tRNAs), whereas the 
remaining genes were on the minor strand (N strand) (4 
PCGs, 8 tRNAs and 2 rRNA genes) (Figure 1). Gene 
order and orientation were the same as the most common 
type suggested as ancestral for insects [18, 23] (Boore et 
al, 1998; Taanman, 1999). 

All genes were closely assembled and only nine 
intergenic spacers were observed. These intergenic 
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Tenebrio molitor L. 


Mitochondrial Genome 


15,785 bp 





Figure 1 Circular map of the Tenebrio molitor mitogenome 
CO I , COII and COIII refer to the cytochrome C oxidase subunits; CytB refers to cytochrome B; ATP6 and ATP8 refer to subunits 6 and 8 of FO ATPase; 


ND1-6 refer to components of NADH dehydrogenase. tRNAs are denoted as one-letter symbols consistent with the [UPAC-IUB single letter amino acid codes. 


Gene names not underlined indicate a clockwise transcriptional direction, whereas underlines indicate a counter-clockwise transcriptional direction. L* and S* 


denote tRNA'"(UUR) and tRNA*"(UCN), respectively. 


spacers were 53 bp in total, with individual sizes ranged 
from 1 to 22 bp. In addition, there were 36 bp overl- 
apping sequences totally in 15 overlaps between different 
genes, ranging from 1-7 bp (Table 1). 

Comparative analysis with two other published 
Tenebrionidae mitogenomes of Tribolium castaneum and 
Adelium sp. NCS-2009 exhibited highly conserved geno- 
me architectures including gene order, nucleotide compo- 
sition, codon usage, and amino acid composition (Friedr- 
ich & Muqim, 2003; Sheffield et al, 2009). 


Protein-coding genes (PCGs) 

Eleven of the 13 PCGs of T. molitor used standard 
ATN (Met) start codon, while the remaining two PCGs 
used atypical start codons (CTA in ND2 and AAT in 
COXI) (Table 1). ATP6, ATP8, ND4L, ND6, and CytB 
genes had the common stop codon TAA; ND1 and ND3 
had the stop codon TAG; while ND2, COX1, COX2, 
COX3, NDS, and ND4 terminated with incomplete stop 
codon T or TA (Table 1). Similar cases have been found 
in other insect mitogenomes (Yin et al, 2012). For 
example, a single T residue was deemed as the stop 
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codon for CO I , COI, ND5 and CytB, and TA was 
deemed as the stop codon for ATP6, ND4, ND4L and 
ND6 in Coreana raphaelis (Kim et al, 2006); similarly, a 
single T was considered as the stop codon for CO I , CO 
II and ND4, and TA residue as the stop codon for ATP6 
in Hyphantria cunea (Liao et al, 2010). Incomplete stop 
codons also exist in vertebrate mitochondrial genomes such 
as in human mitochondria (Andrews et al, 1999). It has 
been demonstrated that incomplete stop codons can produce 
functional stop codons in polycistronic transcription 
cleavage and polyadenylation processes (Ojala et al, 1981). 


Transfer RNA and ribosomal RNA genes 

Twenty-two tRNA genes (one specific for each 
amino acid and two for leucine and serine) were identi- 
fied within the 7 molitor mitogenome. The 22 tRNA 
genes were interspersed throughout the coding region 
and ranged from 60 to 70 bp in length (Table 1). All 
tRNAs, except tRNA*” (AGN), were folded into typical 
cloverleaf secondary structures (Figure 2). The unusual 
tRNA” (AGN) lacked a dihydrouridine (DHU) arm, 
which was replaced by a simple loop (Figure 2). Another 
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unusual feature was the use of TCT as the tRNA*” (AGN) et al, 2011; Liao et al, 2010; Salvato et al, 2008; 
anticodon instead of GCT, which is used in most other Wolstenholme, 1992; Yang et al, 2011). 

arthropods. This incomplete tRNA*” structure has also Twenty-two non-Watson-Crick base pairs were 
been detected in other insect groups (Kim et al, 2006; Li observed in the T. molitor mt tRNAs, including 17 G-U 


Table 1 Summary of Tenebrio molitor mitogenome 





Gene Direction Nucleotide number Size (bp) Anticodon Start codon Stop codon 

tRNA" F 1-64 64 29-31 GAT = m 
tRNA" R 62-130 69 98-100 TTG z g 
tRNA”“ F 131-198 68 160-162 CAT = = 
ND2 F 204-1206 1003 = CTA(M) T(AA)* 
tRNA”? F 1207-1272 66 1240-1242 TCA = = 
tRNAS” R 1295-1355 61 1324-1326 GCA = = 
tRNA” R 1356-1421 66 1388-1390 GTA £: & 
COXI F 1423-2956 1534 2 AAT(M) T(AA)* 
tRNA" (UUR) F 2957-3021 65 2986-2988 TAA = = 
COX2 F 3022-3709 688 = ~ = 
tRNA” F 3710-3779 70 3740-3742 CTT = a 
tRNA“? F 3779-3843 65 3809-3811 GTC = = 
ATP8 F 3844-3999 156 a ATA(M) TAA 
ATP6 F 3993-4664 672 = ATG(M) TAA 
COX3 F 4664-5448 785 -= ATG(M) TA(A)* 
tRNA” F 5448-5509 62 5478-5480 TCC - - 
ND3 F 5510-5863 354 = ATT(M) TAG 
tRNAM F 5862-5927 66 5891-5893 TGC = = 
tRNA“? F 5927-5990 64 5957-5959 TCG = = 
tRNA" F 5990-6054 65 6021-6023 GTT = = 
tRNA*"(AGN) F 6055-6114 60 6077-6079 TCT 5 = 
tRNA" F 6116-6177 62 6146-6148 TTC = = 
tRNA’ R 6176-6239 64 6206-6208 GAA = z 
ND5 R 6240-7953 1714 4 ATA(M) T(AA)* 
tRNA‘ R 7951-8013 63 7981-7983 GTG 7 z 
ND4 R 8014-9346 1333 - ATG(M) TAA)" 
ND4L R 9340-9627 288 = ATG(M) TAA 
tRNA" F 9630-9692 63 9660-9662 TGT z = 
tRNA?” R 9693—9758 66 9726-9728 TGG z = 
ND6 F 9761-10258 498 = ATT(M) TAA 
CytB F 10258-11394 1137 £ ATG(M) TAA 
tRNA” (UCN) F 11393-11458 66 11423-11425TGA - - 
NDI R 11476-12429 954 - ATA(M) TAG 
tRNA‘“"(CUN) R 12427-12491 65 12460-12462TAG = = 
IrRNA R 12491-13771 1281 - z g 
tRNA” R 13774-13842 69 13810-13812 TAC Š = 
srRNA R 13844-14615 772 A = = 
A+T-rich region R 14616-15785 1170 = = = 


tRNA abbreviations follow the IUPAC-IUB three letter code. 
#: TAA stop codon is completed by the addition of 3' A residues to the mRNA. 
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Figure 2 Predicted secondary clover-leaf structure for the 22 tRNA genes of the Tenebrio molitor mitogenome 


The tRNAs are labeled with abbreviations of their corresponding amino acids. Point (e) indicate Watson-Crick base-pairing. Arms of tRNAs (clockwise from 


top) are amino acid acceptor (AA) arm, TyC (T) arm, anticodon (AC) arm and dihydrouridine (DHU) arm. 
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pairs, which form a weak bond in the tRNAs. The 
remaining five were atypical pairings: one mismatch in 
tRNA‘ (CUN) (U-U), one in tRNA‘ (UUR) (U-U), two 
in “RNA” (2 U-U), and one in tRNA”? (A-G) (Figure 2). 
These mismatches found in tRNAs can be corrected 
through RNA-editing mechanisms (Lavrov et al, 2000). 

The /rRNA and srRNA genes of the T. molitor 
mitogenome were | 281 and 772 bp in length, respec- 
tively. As observed in other insects (Boore et al, 1998; 
Kim et al, 2009), these two genes were located between 
tRNA‘ (CUN) and tRNA™ as well as between tRNA” 
and the A+T-rich region, respectively (Figure 1). The 
A+T contents of the rRNA and srRNA genes were 
78.92% and 77.33%, respectively. 


A+T-rich region 

The 1 170 bp A+T-rich region of the T. molitor 
mitogenome was located between the srRNA and tRNA” 
genes (Table 1). Like the two other Tenebrionidae 
mitogenomes, the A+T-rich region also exhibited the 
highest A+T content (85.38%) in the 
mitogenome. Unlike some insect species (Kim et al, 
2006; Yin et al, 2012), the region did not have large 
tandem repetitive sequences; however, it did have some 


T. molitor 


microsatellite-like repeats (for example, (T)ıs, (AT), 
(A)v, and (G)8). The poly-T stretch (15 bp) has been 
suggested to function as a possible recognition site for 
the initiation of replication of the mtDNA minor strand 
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(Andrews et al, 1999; Kim et al, 2009). The A+T-rich 
region in insect mitogenomes is equivalent to the control 
region of vertebrate mitogenomes, which harbor the 
origin sites for transcription and replication (Andrews et 
al, 1999; Yukuhiro et al, 2002; Zhu et al, 2013). 


Nucleotide composition and codon usage 

The A+T content of the whole T. molitor mitogen- 
ome coding region was 72.35%, showing an obvious AT 
mutation bias (Eyre-Walker, 1997), as observed in two 
other published Tenebrionidae species (Friedrich & 
Muqim, 2003; Sheffield et al, 2009). The AT-skew value 
of the major strand was 0.204, indicating the occurrences 
of As in the major strands were greater than those of Ts. 
The mean value of the A+T content of the 13 PCGs was 
69.27%, with a strong A+T bias. The A+T content at the 
third codon position (76.30%) was higher than the first 
(64.75%) and second position (66.77%), indicating that 
the third codon position was most susceptible to AT 
mutation bias (Eyre-Walker, 1997). 

The relative synonymous codon usage (RSCU) in 
the T. molitor mitochondrial PCGs was investigated and 
the results are summarized in Table 2. The four most 
frequently used codons were ATT (Ile), ATA (Met), TTT 
(Phe), and TTA (Leu), accounting for 25.7% of all 
codons in the T. molitor mitogenome. These four codons 
were composed of A or T nucleotides, indicating that 
their biased usage resulted from strong AT mutation bias 


Table 2. Codon number and relative synonymous codon usage in Tenebrio molitor mitochondrial protein coding genes” 























Codon Count RSCU Codon Count RSCU Codon Count RSCU Codon Count RSCU 
UUU(F) 227 1.34 UCU(S) 99 2.28 UAU(Y) 113 1.47 UGU(C) 25 1.61 
UUC(F) 113 0.66 UCC(S) 23 0.53 UAC(Y) 41 0.53 UGC(C) 6 0.39 
UUA(L) 198 2.29 UCA(S) 101 2.33 UAA(*) 0 0 UGA(W) 74 1.61 
UUG(L) 73 0.84 UCG(S) 11 0.25 UAG(*) 0 0 UGG(W) 18 0.39 
CUU(L) 90 1.04 CCU(P) 46 1.39 CAU(H) 31 0.87 CGU(R) 26 1.82 
CUC(L) 9 0.1 CCC(P) 18 0.55 CAC(H) 40 1.13 CGC(R) 2 0.14 
CUA(L) 131 1.51 CCA(P) 61 1.85 CAA(Q) 50 1.54 CGA(R) 26 1.82 
CUG(L) 18 0.21 CCG(P) 7 0.21 CAG(Q) 15 0.46 CGG(R) 3 0.21 
AUU(D 274 1.52 ACU(T) 53 1.07 AAU(N) 116 1.21 AGU(S) 26 0.6 
AUC(D) 86 0.48 ACC(T) 40 0.8 AAC(N) 75 0.79 AGC(S) 6 0.14 
AUA(M) 252 1.67 ACA(T) 98 1.97 AAA(K) 78 1.46 AGA(S) 66 1.52 
AUG(M) 49 0.33 ACG(T) 8 0.16 AAG(K) 29 0.54 AGG(S) 15 0.35 
GUU(V) 91 1.69 GCU(A) 48 1.22 GAU(D) 42 1.2 GGU(G) 50 0.94 
GUC(V) 14 0.26 GCC(A) 24 0.61 GAC(D) 28 0.8 GGC(G) 20 0.38 
GUA(V) 92 1.7 GCA(A) 81 2.05 GAA(E) 53 1.43 GGA(G) 102 1.92 
GUG(V) 19 0.35 GCG(A) 5 0.13 GAG(E) 21 0.57 GGG(G) 40 0.75 





a) RSCU, relative synonymous codon usage. Letters in brackets are the single-letter amino acid codes. Total number of codons, excluding the stop codons 


(*)=3696. 
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(Powell & Moriyama, 1997; Rao et al, 2011). 

The total number of non-stop codons (CDs) of 
PCGs was 3 696. Among all amino acids encoded by the 
13 PCGs, Leu (14.04%), Ile (9.74%), Ser (9.38%), and 
Phe (9.20%) were the four most abundant amino acids, 
three of which were encoded by the AT-rich codons (see 
above), suggesting that the AT bias affected amino acid 
composition (Foster et al, 1997; Min & Hickey, 2007). 


SUMMARY 


Recently, T. molitor has been used as a model to 
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